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CXIV. X-ray Transition Probabilities with Special Reference to 
K Absorption in Lithium 


By J. FRIEDEL 
H. H. Wills Physical Laboratory, University of Bristol* 


[Received August 8, 1952] 


SUMMARY 


To interpret the position of x-ray absorption edges in monovalent 
metals, the author assumed in a previous paper that the positive hole 
created in an inner shell was screened by a bound electron. It is 
shown here that this assumption gives also a correct order for the 
magnitude of the absorption coefficient near the edge. 

In an introductory treatment we take the absorbing system to be 
the ejected electron alone. We use the hydrogenic wave-function for 
the bound state and the plane wave for the continuum. We thus 
obtain approximate formulae for the total jump in the absorption 
coefficients at the K, L and M edges, both for x-rays and in the optical 
region. 

In a further treatment, taking lithium as an example, we treat the 
absorbing atom and the Fermi electrons as a whole, and with this model 
analyse the fine structure of the edge. We show how the screening may 
be obtained from the conduction band. We analyse various possible 
types of transitions producing edges or lines in the absorption, and 
relate them to the structure observed near the edge. To avoid too much 
complexity in treating the Fermi electrons, we enclose the metal in a 
large but finite sphere. Quantitative agreement with experiment is 
not altogether satisfactory, possibly because of inaccuracies in the 
wave-functions used. 

Finally, we describe briefly the reverse process of K-emission in 
lithium, and discuss the possible presence of lines as well as bands in 
the x-ray spectra of monovalent metals. 


§1. INTRODUCTION 


A THEORETICAL treatment of photoelectric absorption by free atoms 
in the non-relativistic region has been given by Stobbe (1930). Stobbe’s 
very simplified model consists of one electron in the field of a bare 
nucleus of charge Z. The wave-functions of both the initial bound state 
and the final free state are then hydrogenic ; and an exact computation 
gives for the K-absorption coefficient (in cm’): 

Z® exp (—4 €arccot €) 


* Communicated by Professor N. F. Mott, F.R.S. 
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Here N is the number of atoms per cm*; ¢)=87e4/3m*c? is the classical 
Thomson scattering factor; the energies hy and H are expressed in 
atomic units* ; and hy=}Z°-+E is the energy of the photon, where LH is 
the kinetic energy of the expelled electron ; finally =}Z°H™. 

To apply formula (1) and a similar one for Ly yy to actual atoms, 
Stobbe makes two corrections. He replaces Z by an ‘ effective’ 
charge Z, which takes into account the screening of the charge acting 
on the K-electron by the other electrons (2,—Z—0-3). He also uses 
for hv the experimental energy of absorption instead of 32Z°+-L. 
Good agreement with experiment is then obtained for solid absorbers 
(Stobbe, loc. cit.; cf. also Bethe 1933 a, and Heitler 1936 a). Some of 
Stobbe’s examples show, however, large deviations near the edges, : 
and agreement is especially poor for light elements. For (solid) lithium; 
for instance, hv=54:5 ev (Skinner and Johnston 1937) and at the edge. 
E=~4-7,ev (Mott and Jones 1936). With these values formula (1) 
predicts a jump of 86-5x104cm~1; the experimental measurement 
(Skinner and Johnston, loc. cit.t) gives about 3104. The theoretical 
value is not very sensitive to £. 

It is therefore of interest to analyse the absorption process more » 
closely. Two criticisms can be made of Stobbe’s treatment: the field 
acting on the expelled atom is certainly less attractive than Z,r—!; 
also the wave functions of the remaining electrons are altered by the 
absorption process. These two effects are most important for light 
elements near the absorption edge. We shall study them in turn, 
restricting ourselves in the second case to the simple case of lithium - 
K absorption. 


§2. SCREENING OF THE PosiTIvE HOLE 


An improvement on Stobbe’s treatment for a gas would be to include 
the effect of all the other electrons by treating the absorbing electron 
as moving in a field of which the potential behaves like Zr! near the 
nucleus and 7~! at large distances. In a metal, however, the effect of all 
the other electrons is more complicated. The conduction electrons after a 
K-absorption process must be in equilibrium in the field of the positive 
hole which replaces the missing K electron. The motion through the 
lattice of this hole is slow and can be neglected; and for transitions 
near the absorption edge, its coulomb field must be screened by the | 
mobile conduction electrons. The screening, after absorption at the - 
edge, is of the type most favourable energetically. In lithium, this is 
provided by a 2s bound electron ; and the energy of the edge agrees 
satisfactorily with the energy of excitation of the Lit ion from the 1s2_ 
to the 1s2s state (Friedel 1952 a). an 

In a metal, therefore, near the absorption edge, it seems a better 
approximation than Stobbe’s is to assume that the expelled electron 


* Atomic units are used throughout this paper. 
} There is an obvious error of sign in the value reported on the figure. 
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moves in the unperturbed lattice. Neglecting the lattice structure, 
we therefore write for its wave function a plane wave exp (vk. r), and 
for the Is electron a hydrogenic wave function : 

RG mmuie Zee eexp (ar ye eer. at (8) 

Let Op=exp (—ix.r)k, be the operator of the transition (Mott and 
Sneddon 1948). In this expression, k, is the projection along the 
direction of polarization e of the momentum k=—iV of the expelled 
electron, considered as an operator. x is the momentum of the absorbed 
photon. We have, with obvious notations : 

(1s |Op |exp {ik . r})=k,(1s, exp {i(k—x)r}). 
With the use of this relation, one obtains easily the K-absorption 
coefficient near the edge (cf. Heitler 1936 b) : 

Li ele 
GZ 2-H) hy’ 

The definitions are the same as in (1), and momentum x has been 
neglected in comparison with k.* 

To take account of the fact that the bottom of the conduction band 
has a negative potential energy —H,, we replace in K the energy change 
(42+) of the expelled electron by the more exact formula : 

hy'=32,/7+E, 


K=4,/(2)N¢,(137)? 


where 
LZ =Z 7 — 2K 
and thus write : 
. BWIA : 


A similar computation gives for the L and M absorptions : 
,2.°E "(32 e—B) 


V2 (hv Ph , 
8Nb.9 ance Za? EO? 
Ly, m= De Tw’ hv’ 
No Z PH? ![3(3Z .)*— 2032, B+ 12K) / 
LINES Serra) as Ta hv , a) 
124/(2)N'¢ Z.,1 B55 Z,2—By? 
“igus Tene ee a 
Z.7H5P(3(4Z,,)'—50(4Z)2H-+ 16 B22 
M se No (137)3 x [ (3 iE : aut: ) ; 
VV 405 56274/2 (hv’) hv 
with hy'=4Z,2-+-E=32+E—E, 


for L absorption, and 
Ww’ =42,2+ H=7,2/+H-E, 
for M absorption. 


* Heitler’s final formula is an approximation for large #; this exact form is, 
however, easily deduced from his computations. 
AG2 
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The absorption jumps computed from (3) and (4) are compared in 
table 1 with experimental values (Skinner and Johnston, loc. cit. ; 
Johnston 1939, Meier 1910). We use Sommerfeld’s value for the energy 
E at the top of the conduction band, and the ‘ effective charges Zi 
according to Slater’s empirical rule (Slater 1930) : Z, is Z—0-3 for K, 
Z—4-15 for L and Z—11-25 for M. Finally, 2) is deduced from the 
cohesive energy E,, the first ionization potential H; and the number 
of valency electrons by the formula (Mott and Jones 1936) : 


E,x~nE,+E,+0-6E. 


Table 1. The Discontinuities t of the Absorption Coefficient at the 
K, L, M edges 


3-0 
0-07 or 1-8 2 


1:3 
2:5 
1-4 
3:6 
2-7 
1:3, 
1-7 
1-5 
0:7 


A correct order of magnitude is obtained down to the optical region 
(Cu Myy, y), except perhaps for L; and M,; absorption. The interpretation 
of the experimental results for those edges is however not certain 
(cf. Johnston, loc cit.). 

We show in Appendix (1) how in a better approximation the perturbing 
potential acting on the expelled electron could be taken into account. 
The correction introduced is small. But to interpret the fine structure 
of the absorption edge it is necessary to take into account the changes 
in the wave functions of the remaining electrons, both in the excited 
atom and in the conduction band. This is a more complex problem, 


which we shall now study in §§3 to 6. Only. the simplest case, that 
of lithium, is analysed. 


§3. Transitions In Many ELEcrron-SystEMs 


The final state of a metal after the absorption of a quantum of radiation 
is by hypothesis a stationary state, and the wave functions of all the 
electrons are of course slightly different from what they were in the 
initial state. We study in this section a simple formula which takes 
this change into account. It is based on determinantal wave functions, 
an approximation valid in the case of light elements where a definite 
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spin may be attributed to every electron. We shall apply this formula 
to Fermi electrons in §4, before using it in §$5 and 6 to study the fine 
structure. 

The initial and final states of the system are thus represented by 
two orthonormal wave functions of the type : 


P(N NMS —1)P half al2) - ++ + Pin(N), 
He Saleh PZ 1)PSP pip (Libga(2) «a N), 
af 


where p;=1 or 2 according to the parity of permutation P;. The 
individual initial wave functions %;, are assumed to be orthogonal to 
each other ; and similarly for the #,,° But the ¢,,, are not assumed to be 
orthogonal to the y,,. 

The matrix element defined in §2 for the transition may be given a 
simple form (Pauling and Wilson 1935 a). We have 


J=JPOp¥, dr 


N - 
=(N 4S (1 MPP, D | dp *(iba(l) dry 
Py, P; Hei 


| bo* DOP A) drs... | Hy" axl) dry. 


We call now P the permutation of the f’s with respect. to the 7’s : 
ize PP, 
J may be written : 


N f 
J=(N !)-12(—1) PP ZP; Z | er*(Ubia(L) dry. - 
P P; Yh 


| ® DOP HI) dry.» | dys*N aN) dry. 


P, only permutes the terms in the last summation over j without 
altering its value. As there are NV! such permutations, we have simply : 


N 
J=3(—1"PD | dsb Pe | Opt, Cae { bateigdes (5) 


To obtain J, we form therefore the determinant | f,;,*,;d7| of the 
projections between all the initial and final individual wave functions ; 
we then introduce operator Op in the terms of the jth row, and add the 
determinants obtained for j=1 to N. 

Formula (5) may be interpreted as follows: the total transition 
probability is the sum of all the possible transition elements fi} >,,*Ope yj dr 
between two individual wave functions, each one combined with all the 
possible projection elements fib>,,*;,, dt of the remaining electrons. Hach 
system of projection elements expresses a one to one correspondence 
between the initial and final wave functions, which may be pictured 
by arrows if the y,,, and %,, are placed on two parallel rows (cf. for 
instance, figs. 5, 6). The sign (—1)’ of a projection is then given by the 
number of crossings of the arrows. 
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We stressed that the initial and final electronic wave functions 
, and x,, are not orthogonal. This is because in the Hartree-Fock 
approximation the central fields V and V’ acting on them are different. 
We have for the two functions : 

Vb = (P—2 Vb 5= 0, 
Vabpy t+ (hk? —2V "iby =0. 
The use of Green’s theorem gives readily : | 
Us e= J tpn ty dT= 2k? —K)™ J by*(V' — Vb y dr. 
The projections may thus be considered as secondary transitions induced 
by the change of potential V—V’—V due to the main transition. 

Hylleraas (1937) has treated the simplest application of formula (5), 
K absorption in helium: the remaining K electron must then be 
projected from a state of the full K shell to a hydrogenic wave function 
round a bare nucleus of charge 2. In metallic lithium, however, we must 
take into account the complex projections occurring in the conduction 
band. We devote the next section to this problem. 


§4. PRogECTIONS IN A FERMI GAS 


In this section, the initial and final states ¥; and Y; are relative to an 
X-ray absorption process in a metal. Our purpose is to calculate some of 
the elements [¥,* dz which we cali projection elements, when Y; and Y, 
are both entirely made up of conduction electrons and do not contain 
any bound electron. We obtain first an expression for the individual 
projection elements w,,,, for one electron (§4.1, formula (9)). We then 
study in some detail the projection of a gas of many electrons (§§4.2 
and 4.3). These elements will be used in §§ 5 and 6 to compute the total 
transition elements. 


4.1. Simple Form for the Individual Projections Un 

We give here an approximate form for the projections U,,= fob, * im U7; 
where 7,,, and yw, are one electron wave functions. 

As this is an absorption process, the initial potential V is that of the 
lattice. We take it as a constant. The final potential V’'=V+8V 
differs from V by a certain amount 8V, the difference between the 
potential of a normal ion and that of an ion excited by the absorption. 
In Li K absorption, for instance, 5V is roughly the difference between 
the potentials of a Lit ion in the states 1s? and 1s2s (ef. §5). 

We assume 8V to be spherically symmetrical and take the metal to 
extend over a sphere of large radius R. The initial and final wave 
functions are thus of the type : 


1 /21+1\12 ¢,,(k, 1) 
Bb im= 2 ( mre ) ane care P,(cos Cole, 


A /21+1\12 6,,(k’, r) 
Brn = Fee 


P,(cos@,)x; 


me 7R 
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where the spin wave functions y and quantum numbers / must evidently 
be the same in the two functions. P, is the Legendre coefficient ; 
6, the angle of r with momentum k. The radial parts satisfy : 


Foe og adi ten 
a =e (ie ) bm=9, 
a4, Ul-+1) 


me + (1 5 = 23) b,=0. 


They are normalized so as to give for large r: 4,,~sin (kr—27m) and 
$,~sin {k'r--n(k’)—3l7}, where 7 is the phase shift produced by 8V. 

* A is a normalizing factor which takes into account the displacement 
of charge in the state of momentum k’ due to the perturbing potential 5V. 
If 6V attracts locally a charge q, this reduces the density of charge 
outside the perturbed region by a factor 1—q. One shows easily that 
g=(1/R) [dn(k’)/dk’] (Friedel 1952 a). Thus: 

, 1 dy (k’) 

Fai grime ce eee aes gs WO) 


There are 2/1 independent functions ¢,,, of given J, and similarly 
for %;,. Only the functions with parallel momenta k and k’ have 
non-vanishing projections, so w,,,, reduces to: 


thy =4AR-U hk’? — Ke)! [bn (k'r)8 V4 ( ker) Fe Ge 
J0 


Formula (7) is exact. We replace now the upper limit R of the 
integral by infinity. This adds a term, the magnitude of which oscillates 
as R increases, with an average value zero. This term may be neglected. 
Then developing ¢,,(kr) and ¢,(k’r) successively in powers of (k’—k), 
and using the exact relation (Mott and Massey 1949 a) : 


J Saldr)BV by (br) dr—— Ble sin {nh}, 


we have: 
2[ J,(k'r)BV4(kr) dr—=—3[k' sin{y,(k’)} +h sinfy(k)}+OH'—K). (8) 
0 


A more detailed discussion (cf. Appendix (2)) shows that this is a good 
approximation if k’ is not too different from & and 5V not too large. 
. The boundary condition #,,,(k’R)= (kR)=0 gives : 

. kR~3lr+mr and k’R~4ka+nr—7(k’), 
where m and n are integers which can be used to define the initial and final 
momenta, ‘Thus: 
= _ Afke sin n,(k)-+’ sin 9 (k')] 
mn (hk) (k')+(m—n)r] ° 
where A is given by eqn. (6). 

Formula (9) shows some general properties of the projection 
‘elements wt, First, if 5V vanishes, w,,, becomes a 6 function 5,,,, for 
free electrons have orthogonal wave functions. If now 6V +0, u,,,, has, 


(9) 


U 
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for every m and n, a finite value, which is function only of the phase 
shifts 7, and does not depend otherwise on 5V—this at least if JV | 
and |k’—k]| are not too large. The function |%nn| has a single maximum, 
for that value of n—m which is nearest to 7,/7 ; it tends to zero as m—n 
tends to infinity. 

When the metal becomes infinite (R infinite), u,,,, will be negligibly 
small except for infinitely small values of k’—k. For such values : 

Sin (hk) 
mn— n(k)+(m—n)r : 
|v,,.»| is then a 8 function if y,is a multiple of 7, But it has its broadest 
maximum when d|w,,,,|/dk=0. This occurs when either dy,/dk=0 or 
4,2 (A+4)z (A integral).* 

Finally, w,,,, changes its sign, as m—n varies, for 7,—A7z and immediately 
after the maxima of |1,,,,|. 

These remarks will be used in the next section. We study now the 
projection elements for a system of N Fermi electrons. We start with 
a case where the final state is the most stable possible. 


U 


4.2. Fundamental Projections I(p/N) 
4.2.1. General Description 
We restrict ourselves to a definite spherical component / and a definite 
spin, and study projections of a simple kind pictured in fig. 1. In the 


Fig. 1 
‘ n 
N —@—— 
Mf —e— Nty-! 
l : 
—@o—— 
Stk 2 =e 
5 SO 


ig peeks —e—_ jt2 
le eee Ne —@o—_ jl 


(a) (0) 
Scheme for projection I(p/N). (a) Initial, and (0) final state. 


initial state, N—1 electrons and a hole occupy the N states of lowest 
energy up to momentum ky;. In the final state, the N—1 electrons occupy 
the states of lowest energy in the conduction band. In fig. 1 (b) 8V was 
assumed to subtract 7 bound states from the bottom of the band, so 
that the first state in the band is n=)+1. The total projection is 
denoted by J(p/N) when the hole is in the pth state. Projections of this 
type will be used in §5, the hole being replaced by an electron which, 
in the final state, drops into a bound level. 

The quantity (—1)? J (p/N) is deduced from the matrix (%,,) 
(m=1 to N, n=j+1 to N-+j—1) by suppressing the m=pth row (§3). 


* More exactly tglmi+(m—n)r)=[n1+(m—n)a]?. 
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Its large terms are near to the main diagonal, and its absolute magnitude 
is evidently smaller than one. Some general characteristics are easily 
deduced from the behaviour of w,,,,. 

Let k=ky(p-+-31)/(N-+41) be the momentum of the hole p. If first 
9 (k)=Ar (A integer), the projection I,(p/N) tends to zero as R tends to 
infinity. For, when Rk becomes very large, the projections of the 
n=(p-+d)th column tend towards a 8 function (cf. §4.1) : 


Um, pa = Oni ps 


Therefore in the determinant I,(p/N) where mp, all the terms of the 
(p-+A)th column tend to zero, and hence I (p/N)=~0. As an application, 
I,(0)~0, for at the bottom of the eondustion band 7 (0)=jz7 if there 
are j bound states (Mott and Massey 1949 b). 

On the other hand, dJ,/dk=0 if d|u,,,,|/dk=0 near the main diagonal. 
Hence (§ 4.1) either dy (k)/dk=0 or 7 (hk) & (A+4)z (A integral). 

In conclusion, the projection is maximum if the momentum k of the 
hole is in a zone of maximum ‘ misfit ’ with the quantized momenta k’. 
This happens if the momentum & falls just in the middle between two 
values of k’ (y,=(A+4)z) or if the misfit, hence the phase shift 7,, are 
stationary (dy,(k)/dk=0). The projection on the contrary vanishes if 
the momentum k of the hole is just equal to one of the momenta 
k’ (n (k)=2n). 

4.2.2. Simple Example 


Keeping the band width ky constant and allowing F to increase, 
each term in the development of I(p/N) decreases, ‘and tends to zero for 
an infinite volume of metal. But we believe that the number of terms 
in Z(p/N) must increase in such a way that the total tends to a finite 
value. This point is illustrated, if not proved, by the study of a simple 
example. 

When 8V has the form of a spherical well of radius a and constant 
depth —D, an exact and straightforward computation gives : 


Ad NEY DG sin (K’—k)a_ sin (K’+k)a 
Van B—t| Kb OKTHE 
with 
X= R-"[a—(2k’)-1 sin 2K'a]+sin? K'a sin~(k'a+79) 
[I —R-1fa—(2h')-1 sin 2(k/a-+n9)} 
k’ cotg (no +k’a)=K' cotg K'a 
and 
K?=k?+2D. 
We took as an example (in atomic units) : a=3 and D=2/9. There is 


then one bound s-state and ,(k) varies as pictured in fig. 2. 
I(p|N) was computed for ky);=9-63 and k=5, 10, 15, 20 and 25, hence 


N=1, 2, 3, 4, and 5 (fig. 3). 


1124 J. Friedel on the X-ray Transition Probabilities with 


For N>3, the curves seem to converge rapidly towards a Ae aoe 
represented by the full curve. This convergence has been chec i a 
other examples and seems reasonable, for there are already abou 
100 Li atoms in the sphere of radius R=15, and more than 450 for N == on 
The small scatter of the points may be due to slight errors in the 
computation. As pointed out in §4.2.1, L(0)=0 and dI/dx=0 for 


Ho=n7/2. 


Fig. 2 
1 
1 
to 
3% 
ie) k 
k—> or 


no(*) in a simple example. 


Fig. 3 


a 


! 
Oo X=P/N —> 


I(p/N) for various N’s. 


Points marked A for N=1; © for N=2; «for N=3; x for N=4; 
®@ for N=5; O for N=3 by approximation (9). 


Failing a rigorous proof, we feel therefore justified in evaluating the 
projections [(p/N) for small values of R(N=3 or 4). The values obtained 
using approximation (9) are in fair agreement with the exact ones 
(fig. 3 for N=3). We shall therefore use this approximation in all 
further computations. 


4.3. Projections over Excited States 


Projections of the initial state fig. 1 (a) over excited states will be 
used in §§5 and 6. : 
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In the final state pictured in fig. 4, the electron with highest energy 
is separated from the rest by q holes. We denote by [%(p/N) such a 
projection. 41) is deduced from I(1) by replacing the terms Up», y4j—1 
by Um, N+j+q-1 Which decrease rapidly as q increases. So only projections 
with small ¢ are of importance. As for J(1), one can show that, for given 
q, £%(1) tends rapidly towards a limit as the volume of the metal increases ; 
and it is sufficient to evaluate it for small R. J %(p/N) is then deduced 
by the approximate relation : 


I%p/N) = {14(1)/T()}1(plN). . se... (10) 


More complex excited states may be studied in the same way. They 
are of little importance in the transitions studied below. 


Fig. 4 
——_@—— N+} +q-1 


—_@——-_ +2 
Be j+t 


Final state for projection [4(p/N). 


§5. Martn Apsorprion EpGe in LiTHIumM 


In a metal the bound states of the excited ion and the Fermi electrons 
are seriously affected by an x-ray absorption. They must all be taken 
into account in a detailed study of the transition; and, according to 
formula (5), we must add the various matrix elements obtained when 
one electron makes a transition to any one of the new electronic states, 
all the other electrons being ‘ projected ’ over the remaining states. 

In lithium, three of these elements are preponderant: the transition 
may affect the inner bound electron ; this electron is then thrown into 
the continuum (fig. 5) or goes over to the new bound state which screens 
the hole created in the inner shell (fig. 16) ; the transition may also occur 
between a Fermi electron and the screening bound state (fig. 6). These 
transitions are appreciably large only when the Fermi electrons remain 
in their most stable state, thus when there is no hole left in the conduction 
band. The transition pictured fig. 16, where no electron is thrown into 
the continuum must correspond to an absorption line, while the two 
others give an absorption edge. The transition of fig. 16 is actually just 
an inner excitation of the monovalent ion, and only occurs when allowed 


by the usual selection rules. 
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5.1. General Description 


In lithium, the energetically most favourable screening of the hole 
in the K shell is by a 2s electron with spin parallel to that of the remaining 
K electron. For this configuration gives the most stable form of the 
lithium ion in the state 1s2s (Bacher and Goudsmit 1932). 

We picture therefore the transition near the main absorption edge 
(hv=54-5 ev) in the following way. The initial state contains the 
K shell (1s?) and » Fermi electrons (f;). The final state contains a 
ls and a 2s electrons with parallel spins (1s’, 2s’), one electron in the 
continuum (c’) and n—1 Fermi electrons (F,’). The dashes emphasize 
the fact that the potential provided by the excited ion is different from 
that of a normal ion; hence the quantized states after absorption 
(1s’, F,’...) are different from those before absorption (1s, F’;). 

A possible scheme of transition is pictured in fig. 5. It may be written 
with an obvious notation : 


T= (c,'|Op|ts)(1s’, 1s)(28’, F(x) )To(w)Z,3(1)159(1) . . . To(1)L,3(1)72°(1). . -« 


Fig. 5 


energies ——+> 


(a) Cb) 


First schema of absorption with 1s’2s’. (a) Initial and (0) final states : 
—> transition ; — — > projection. 


The terms barred are those with spin antiparallel to that of 2s’ and 
F(x) denotes a Fermi state of s symmetry with momentum k=zky. 
The factor in J, is obtained assuming that at the end of the transition the 
Fermi electrons are in their most stable state (cf. §4.2), i.e. there is 
no hole left in the Fermi gas at the end of the absorption. 

Such transitions, whatever the value of 2, are between states of the 
same difference of energy. The total transition probability is obtained 
by integrating 7*7 over x after multiplication by the density of 
states p, of the Fermi electrons F,(x) : the Fermi electron which happens 
to be nearest to the atom during the absorption drops into the 2s’ state ; 


its energy may have any value in the range of the conduction band.* 
ene ap Le 
* As [(0)=0, the electron cannot come from the bottom of the conduction 


band. We shall see that the i ibuti 
ie Gara Caeser e main contribution comes usually from the top of 
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The total transition probability is therefore finite, although each 
individual 7’ tends to zero if the volume of the metal becomes infinite. 
Transitions of the type pictured in fig. 5 therefore correspond to a finite 
absorption edge at the critical frequency. 

There are many possible schemes of absorption other than 7. We 
shall however neglect those where the electronic transitions are between 
Fermi electrons, which are certainly unimportant. Conservation of 
spin and usual selection rules for electronic transitions then restrict the 
possibilities. The transitions leading to the most stable states, i.e. those 
which leave no hole in the Fermi gas, add up to two terms corresponding 
to two different configurations : 


A(x)=(c,'|Op|1s)[(1s’, 1s)(28’, F,(2)) 
—(1s’, F,(x))(28’, 1s)|Zg(x)1,3(1) .. . Zo(1) . 
Blew; y)= {(28'|Op|F,(y) Ls’, 18)(¢,, 1s)Z9(1)—(c,’, 18)(1s’, F(x) I o(x)] 
—(1s’|Op|F ,(y))(c,’, 18)[(28’, 1s)Zg(1) 
(250 (or). 1), Is)Z g(x) 144) P20) Pa8(1). « . Zo(1).. - 
The various signs are easily obtained by counting the number of crossings 


of the transition arrows (§3). The total transition probability is obtained 
by integration over x and y: 


v1 Ly eal) 
| A*4p, det | | B*Boapy de dy. Bey a1 1 
0 0/0 


Transitions leading to holes in the Fermi gas will be considered in §5.4. 
But we shall first study the main term (11). 


5.2. Reduction of Formula (11) 

We use spherical harmonics for the states in the conduction band 
and the continuum, and hydrogenic wave functions for the bound 
states (cf. (2) for 1s). 

Thus : 
ela) ek iret sn Ky, 
Cp’ = V-123tk’—*r- (sin k'r—k'r cos k'r)P, (cos 0), (12) 
29’ = (44/(27)) 12’ ,2/2(2—Z,'r) exp (—3Z’,r). 
We have then : 
(1s’, Is)=8(Z,+Z)*°(Z,Z)°”, 
(28’, Is) =24/(2)(Z,+32'.) (Z.2'.)?"(Z.—Z'e), (13) 
[29’, F,(a)]== —V-2224/(277)(4Z' 2+?) 90" 5 ?(4Z' P —k?). 

As in §2, we take into account the negative energy — EH, of the bottom 
of the conduction band by replacing in this last projection Z’, by 2’, 
chosen so that $Z',?= 2’ ?—E). 
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We may attribute to the screening electron 28’ the difference of energy 
between the absorption edge hv) and the second ionization potential 
of lithium Lill (cf. Friedel 1952 a) : 

E2,™ Lill—hyp. 
We have hvy>=54-5 ev (Skinner and Johnston loc. cit.) ; Li II=75-26 ev 
(Bacher and Goudsmit loc. cit.). Hence Hy.= $Z’ 2=20-8ev and Z ow 2-4. 

(13) then gives (1s’, 1s)=0-9975, (28’, 1s)=0-0445 and in (11) only 
the terms containing (1s’, Is) are of importance. The absorption 
coefficient thus reduces to : 


os fA , 
ry KI,A(1)Zg (1)... Tp2(1) -» (18", 18)? [es F,(«))? 


[T,2(x)E42(1)lo2(1)2a2() on G01 ee) 2 
with 
Pp, de= (20?) Vie dk. 
K is given by (83) and the various projections by (13). 
The two terms in 7, are described in figs. 5 and 6. 


(a) (b) 


Second scheme of absorption with 1s’2s’ : (a) initial and (0) final states. 
— transition ; —--> projections. 


5.3. Jump in Absorption at the Edge 


We shall now compute (14) when the electron thrown into the 
continuum is just at the top of the conduction band : k’=hy,. 

We have to evaluate the projection factors I,, hence the phase shifts nt 
produced in the Fermi gas by the perturbing potential SV. The main 
part in 6V is due to the replacement of a normal Li*1s? ion by an excited 
Is'2s’ one. We shall neglect the other terms, namely the contribution 
from the charge displaced in the conduction band and the corrective 
term in Wannier’s theorem (cf. Friedel 1952 a). These should be considered 
to obtain a more correct self-consistent value of SV. 

We divide 8V into two parts, an electrostatic 5V, and one due to 
exchange 6V, : 


8V=8V,+18V,. 
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With (2) and (12) we have: 
OV =r *(2(1+Z,r) exp (—2Z,r)—(1+Zr) exp (—2Zr) 
pita 64-42) ,F*) expd—Z' 7)! 
For 5V,, we used Slater’s average value (1951) : 
L272 a? (r8V.) Je" 
We checked this formula a a direct evaluation of the exchange term 


between the 2s’ electron and the Fermi electrons with s symmetry 
(cf. Fock and Petrashen 1935). Figure 7 shows a fair agreement in 


6V,.2 (15) 


fF atomic 
units 


Term in 6V, due to electron 2s’. —— Fock ; - —- Slater. 


Fig. 8 Fig. 9 


r atomic 
units 


-2 
Perturbing potentials for 1s’ 2s’. Phase shifts for 1s’ 2s’. 


the interesting range of radii (r>2 for k<ky=0-:59). Thus for Fermi 
electrons with spin parallel to that of 2s’: 


8V = 1-272[274Z, exp (—3Z,r)— 
{4Z3 exp (—2Zr)+42' 2(4—Z' +42’ ,r*) exp (—Z’-r)}3*] 
and, for Fermi electrons with opposite spin : 
SV ,=1-272 x 228 Z, exp (—3Z,r). 


SV and 8V are given fig. 8. Owing to exchange with 2s’, the 
potential 5V is attractive enough to subtract from the conduction band 
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one s and one p state with spins parallel _to its own: 7o(0)=m(0)=7 
(fig. 9). The 7, are small. As our dV and 5V are not self-consistent, the 
7, do not exactly satisfy the screening condition (Friedel 1952). We have 


actually : 
a 3>Z (20+ 1) [9 (ke) +k) ]—4= — Lb 1 £0. 
l 


From 7, we deduce J,(x) as described in §4.2. Figure 10 gives J, and 
I,. The others are practically 1 for z=1 and 0 elsewhere. 

With these values for the J;, relation (14) gives : 

7,=(0-04-+0-04,) x K=0-11,; x 104. 

The two terms correspond to the schemes of figs. 5 and 6. They are 
of the same order of magnitude: it is here equally easy to subtract an 
electron from the s or p part of the conduction band, because the 
corresponding phase shifts 7) and 7, hence the projection factors 
I, and J,, are similar. We see also that the major contribution of the 
Fermi electrons comes from the top of the conduction band. 


Fig. 10 Fig. 11 


Ff atomic 


re) ° units 
} 
1, : + 
- x —> ! ae 
Projections J; (x) for 1s’ 2s’. Perturbing potentials for 1s’ 2s’. 


5.4. Variation of r with the Energy of Absorption 


_ 7, represents the two major terms in (11). The others contribute 

to only a few per cent of 7,, as is easily computed, and may be neglected. 

But to every ‘ fundamental ’ transition so far studied there corresponds 
a series of transitions of the type described in $4.3 leading to an excited 
state of the Fermi gas. Here, the only ones of importance are the [4%(x) 
with g small (<q =3 or 4). So on the short wave lengths side of the edge, 
the absorption coefficient increases rapidly from 7, to a saturation value 7. _ 
The change occurs within a finite number cf Fermi states, and hence 
in an infinitesimal range of frequencies ; and the coefficient experimentally 
measured is +. Using relation (10), we must therefore replace every 
projection [?(a) by 


P(a)+ 2 [eee = Peo (42 Eran ) 


This is greater by 1-12 for J,2 and 1-12, for J,2; thus 7 is increased by 
a factor 1-6: 


71-67, =0-2 x 104, 
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Finally, + varies with the energy of absorption proportionally to K. 
Thus from (3) 


dr 3 4 i 


7 should therefore increase slightly on the short wave lengths side of the 
edge, up to a maximum at a distance from the edge found to be about 
25 ev. However, before this energy is reached, other features occur, 
probably due to screening by excited electrons. 


§6. SCREENING BY ExcirEep Bounp States in Liratum 


In a previous paper (Friedel 1952 a) we interpreted tentatively the fine 
structure on the short wave lengths side of the main absorption edge 
as due to screening by more excited electrons. We show in this section 
that one should expect an absorption edge for every configuration of 


the inner ion (1s’2s’, 1s’2s’, 1s’2p’, 1s’2p’...) and an absorption line 
when a direct excitation of the ion is possible (Isls+Is’2p’....). 
Fig. 12 Fig. 13 
is ! 
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k-— x— 
Phase shifts for 1s’ 2s’. Projections for 1s’ 2s’. 


6.1. Secondary Absorption Edges 

The same analysis as for the main edge gives the following results. 
Belcl. Te 2s' 

The energy difference between the Is’2s’ and 1s’2s’ states of 
Lit ion is about 1-5 ev (from optical data, cf. Bacher and Goudsmit 
loc. cit.). Attributing to electron 2s’ this difference, we have 
17’ 2=E y= 20-8—1-5=19-3 ev and Z’,=2-38. The coulomb potential 
8V, is therefore slightly more attractive. The exchange term due to 1s’ 
is now added to 5V and subtracted from 6V (fig. 11). So the 7 are 
positive and 7, smaller for large k (fig. 12), and J) has a maximum for 
no=in (cf. §4.2). Ig and I, are given in fig. 13. J,(1)~0-98 and the 
other J(1) may be taken as 1. 

The possible transitions are, with a change of spins, those pictured 
in figs. 5 and 6. And (14) gives: 

712 (0-11+0-04) x 104=0-15 x 104, 
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Projection over excited states (§5.4) multiply 7, by a factor 1-7. Thus: 
71-77, =0-255 x 104. 
6.1.2. Is’2p’ 
We take for 2p’ a hydrogenic wave function : 

2p’ =[44/(27)]-1Z"",2r exp (—3Z” r) cos? . . . {16) 
where the value of Z’’, is deduced from the position of the absorption 
line L, (fig. 17) which we interpret as corresponding to the same 
screening Is’2p’. We have hy=58-5 ev (Skinner and Johnston loc. cit.), 
hence Z’",~2-22. 8V,7,, I) and I, are very similar to those for 1s’2s’, 
and-J4(1)==0:9. 

Thus practically only Fermi electrons with s or p symmetry can 
drop into the bound state. And the absorption is the sum of two 
transitions similar to those pictured in figs. 5 and 6, with a permutation 
of s and p components in the Fermi gas : 


Tee OKT ACL) 15201 ee ee el 
[ 12D’, F(x) PUPP) +LA°ADL2@)) pe ae. (17) 
0 


where factor 3 is due to degenerescence of 2p’ states. With (12) and (16) 
one finds : 
(2p’, F(x) P= 8a(38V) IZ" Phe (gZ" P+ k?)-® 

where Z”’, has been replaced by Z’’,, such that 12"? =41Z" 2H, (ef. §2). 

Finally : 

7, (0-03-+ 0-05) 104=0-08 x 104. 
Projections over excited states give a factor 1:5. Thus 
71-57, =0-12X 104, 

6.1.3. 1s’2p’ 

We neglect the difference of energy between this case and the 
preceding one (less than 1 ev from optical data), and thus use the same 
value for Z’’,. SV is then similar to that for 1s’2s’ but stronger. Hence 
no(0)=27 (fig. 14); ZI presents a maximum for small k and is nearly 
zero at the top of the conduction band (fig. 15). 

The same formula (17) as before is valid for 7, and gives : 

7,~(0-000;-+-0-01,)104=0-01, x 104. 
With a factor 1-5 for excited states, the total absorption coefficient 7 
is small : 
71-57, =0-02 x 104, 


6.2. 1s’2p’ Absorption Line 
With 1s’2p’ a direct transition of the inner Lit ion is possible (fig. 16) : 
(Y,*] Op |¥,)=(2p'| Op |1s)(Is’, 1s)Z(1)Z,9(1) ....Zp(1).... 


There is no electron c’ thrown into the continuum ; so this gives an 
absorption line. 
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The maximum strength of the line is given by (Heitler 1936 c) : 
NihwWyp%y" 
(hv)Phy * 
N is the number of atoms per cm*; hv the resonance energy and hy 
the line width (both expressed in atomic units) ; and w,, the transition 
probability, proportional to hv| (Y%*| Op |¥;) [?. 
With (2) and (16), we have 


Shp (2,2) 


ty= 2a (137)? (18) 


hw »=6 — SS OOo Fr (I8’, 18)? Le? eC IN eau oC 1) tees 
Wab °3x(i37) (Z,+42"7,)8 48 ? Is) I) (1) I, (1) Ty°( ) 
Fig. 14 Fig. 15 
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Phase shifts for 1s’2p’. Projections J; for 1s’2p’. 
Fig. 16 Rigo 
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@) (b) hv 
Schema for absorption line with 1s’2p’. Li K absorption, after Skinner 
(a) Initial and (6) final state. and Johnston. 


The factor 6 comes from the two K electrons and (2x3) 2p’ states. 
An experimental line width y~ 1-5 ev (Skinner and Johnston loc. cit.) 
‘and the values obtained in §6.1.2 give » 
Wap= 0-45 X 101° CGS. 
Hence 
Ty, 12x 104, 


4H2 
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Projections over more excited states in the Fermi gas (cf. § 4-3) introduce 
a factor 1:5. Thus finally 
inn 1: 8 x 1 04, 


The large broadening y may be due in part to Auger effect, the 2p’ 
bound electron making a transition to the neighbouring 2s’ state when 
colliding with a conduction electron. The same broadening should 
affect the 1s’2p’ and 1s’2p’ edges. Some of y may also be due to lack 
of resolving power in the instruments ; the form of the line is actually 
different in the two measurements by Skinner and Johnston. 


§7. COMPARISON WITH EXPERIMENT 
Skinner and Johnston’s measurements (loc. cit.) give three successive 
absorption edges E,, E, and E;, E, and E; starting with somewhat broad 
and weak lines L, and L, (fig. 17). 
E, probably corresponds to the formation of 1s’2s’ Lit ion, and the 
computed value of 7 turns out somewhat too small (table 2). 


Table 2. Comparison of the Fine Structure with Experiment 


AtX 10-4 E, E, L, 
exp. 1-0 1-2 0°5 
comp. 0:2 0-4 1:8 


The same happens with E,, if we attribute it to the 1s’2s’ and 
broadened Is’2p’ and 1s’2p’ edges. The computed value for L, on 
the other hand is too large. We believe this is because the 2p’ wave 
function is actually somewhat swollen and deformed by the surrounding 
lattice. For E, and Ly, the approximations concerning both bound 
and free electrons would certainly be unsatisfactory. Finally the 
absorption seems more or less constant after each edge, as predicted 
in §5.4. 

The computation gives therefore an acceptable order of magnitude, 
but no very good agreement, the calculated relative value of the L, line 
with respect to the edges being about ten times too strong. (cf. table 2). 

The fine structure E,, L,, E3;, L,...could on the other hand be a 
Kronig structure (Jones and Mott 1937). The first peak L, would 
then correspond to the peak in the first zone arising from reflection on 
the 110 planes. As pointed out by these authors, this peak would be 
observed in a K absorption only if the first zone is ‘ inverted’ so as to 
have s character at the bottom but p character at the top. In that 
case the K emission should have a sharp maximum at the edge. But 
actually no such maximum is observed (Skinner 1940). Also a fine 
structure similar to that of lithium, with same position of the secondary 
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maxima, occurs for LiOH which has of course another crystal structure 
(Skinner and Johnston, loc. cit.). These two facts seem to suggest that 
the structure observed is not a Kronig structure. 


§8. QuaLiTatrvE Description or Li K Emisston 


The transitions occurring in Li K emission are more or less the reverse 
of those for absorption. We point out some features of the process. 

The initial state in the emission contains mainly the most stable 
Lit ion, 1s’2s’, and the main transitions for the emission band correspond 
to those of figs. 5 and 6 for the absorption. They are pictured in figs. 18 
and 19. 


Fig. 18 Fig. 19 
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FT 
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ee Tels a My 
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© i) ’ @ (b) 
First type of transition for the emission Second type of transition for the 
band. (a) Initial and (6) final state. emission band. (a) Initial and 


(5) final state. 


An electron drops from the part of the conduction band with spin 
antiparallel to that of 2s’. This part is little perturbed (8V and 7, 
small in figs. 8 and 9) and the projections J; (#) are zero for x ~ 1 (§5.3). 
Conversely the perturbed Fermi gas with a p (or 8) hole of momentum 
k<ky must have a negligible projection over the stable unperturbed 
gas. A/p(orS) hole must therefore be created, with momentum practically 
equal to that k of the electron dropping into the Is level. As for 
absorption, the process leaving a p hole (fig. 18) is probably the more 
important. s holes can nevertheless be created (fig. 19). 


§9. CONCLUSION 
PossiIBILITy OF LINES IN THE FINE STRUCTURE 


In conclusion, experimental results seem compatible with a screening 
of the hole created in the K shell by absorption. In lithium, an 
absorption edge should correspond to every possible configuration of 
the screening and a line to every possible direct excitation of the inner 
ion Lit. The computed strength of the absorption jump at the edges 
agrees in order of magnitude with experiment. It is smaller for more 
excited screenings, because the Fermi electrons are more perturbed. 
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Thus the edge computed for a 2p’ screening is smaller than for a 2s’ 
one; for 1s’2p’, smaller than for Is’2p’. In our interpretation of 
Johnston and Skinner’s measurements, the 2p’ line is weaker than as 
computed, probably because the 2p’ state is somewhat swollen and 
deformed by the surrounding lattice. = 
Similar absorption lines should appear in various transitions of 
monovalent metals. At the edge of Ly yz, absorption, a weak line 
may actually be seen in Na (O’Bryan 1940), but not in Cu (Coster and 
Veldkamp 1932) nor in Au (Veldkamp 1935). We think that the line 
is not observed with heavy elements for the following reason: as the 
atomic number increases, the radius of the outer screening electron 
increases. The overlap on it by the inner L electron decreases therefore 
more quickly than the overlap on the states in the conduction band, 
which keep a nearly constant density near the origin. This decreases 
the importance of the lines (transitions of the type pictured in fig. 16) 
relative to the edges (cf. figs. 5 and 6). The rapidity of this decrease 
may be deduced from a remark by Bethe (1933 b) that, for a given atom 


Fig. 20 
s Pp s p 
2p 
©) Cb) 


Li, m line when screening is by Fermi electrons. (a) Initial and (b) fual states. 


or ion, the relative intensity of say transitions 2p—ns with various n 
are nearly the same as for a hydrogen atom. Thus the ratio of the 
oscillator strengths of transitions 2p—3s, 4s, 5s, 6s to that for transitions 
towards the continuum (a constant in Bethe’s case) are proportional 
to 1, 0-2, 0-08 and 0-04. 

As in lithium, the intensity of these absorption lines is also probably 
reduced by perturbation of the ns bound state by the surrounding 
lattice, especially in noble metals with their large d shells. Furthermore, 
the screening electron is probably replaced in noble metals by a heaping 
up of charge in the conduction band (Friedel 1952 a). In other words, 
the perturbing potential 5V is perhaps not strong enough to subtract 
-a bound state from the bottom of the conduction band, but increases 
only the density of Fermi states in k space, mainly in the s component. 
Transitions of the type pictured in fig. 16 are then impossible, and one 
may wonder whether any absorption line may be observed in those 
conditions. Actually a transition of the type pictured in fig. 20 may now 
take place ; owing to the crowding of the Fermi s states, the projection 
-of the Fermi electrons pictured in fig. 20 may have a finite value. It will 
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correspond to an absorption line, for projections over more excited 
states of the Fermi gas will be vanishingly small. By integration over 
all possible final states, one obtains a transition probability with a finite 
value, which we believe to be of the same order of magnitude as if 
screening was by a bound state. If this is so, x-ray evidence cannot 
help in differentiating the two types of screening. Thus for instance 
we interpreted the secondary peak observed at about 25004 in optical 
absorption of copper as a direct 3d—4p excitation (Friedel 1952 b), but 
the 4p screening is certainly obtained by a heaping up of Fermi electrons 
and not by a bound state. We think therefore that it is only because 
of the weakness of the transition 2p—4s that no line is observed at the 
edge of Cu Lay ar- 

Finally, in Na Ly yy absorption, the line should correspond to the 
less stable singlet state Ly, (configuration 1P, of Nat), as is indeed 
observed. In Na Ly yy emission, the main contribution to the band 
is thus the more stable triplet state Ly (configurations 3P,, 3P,, 3Pp 
of Na*), for which no direct transition of the Na+ ion is possible. The 
contribution by Ly; should contain a line at the emission edge. But 
this term is very weak in Mg (Skinner 1940) ; and actually no line has 
been observed either for Na Ly, yy (Skinner loc. cit.) or for the analogous 
K My, m emissions (Kingston 1951). 


The author wishes to thank Professor N. F. Mott for enlightening 
discussion and advice and the Direction des Mines for allowing him 
to stay in Bristol. 


APPENDIX 


(1) APPROXIMATE CORRECTION FOR THE EXPELLED ELECTRON 


In computing K (§2), we treated the ejected electron c’ as a plane 
wave. We indicate here how one could take somewhat into account 
the perturbing potential 6V acting on it. 

We use a simple model where 5V is a spherical well of potential of 
radius a and constant depth, at the centre of a large sphere of radius R 
representing the metal. For an s state with momentum k and if ka>z, 
the average radial density of charge may be considered as constant 
both inside and outside the well. When 5V=0, the charge of the state 
inside the well is thus 

Ga ah. 

Let g+4q be the charge of the state inside the well when 5V #0. 
The radial density of charge of the state outside the well is now multiplied 
by (l1—q—4g)(1—q)"!. The density of states per energy unit is 


multiplied by 4 
1 dy\-1 
(.- Rit) 
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(Friedel 1952 a). And the product of the two, that is the total density 
of charge per energy unit outside the well, must remain constant. 


(Laue 1914). Thus 
1—q— 4q = (1—q)(1—R dy/dk) 
or (q+ 4g) 7? = 14+-(¢'1-1) Rt did. 


The total density of charge per energy unit inside the well is then 
multiplied by a factor 


f= (1+(@1—B*) dy[dk\[1—R~ dna] 


or, for large R 


If 8V has a more general form, we take for a the radius giving, in the 
WKB approximation, the major contribution to the phase shift 7. 
For the 5V obtained in this paper for lithium, f is then always slightly 
smaller than unity, and therefore reduces + somewhat. 


(2) Form or THE PrRosection ELEMENTS ¥,,,, when |k’—k| 1S LARGE 


The development (8), carried one step further, gives 


Olea oe ae eV E- at) PE CEE) 


Ob n (Kk, a 
=, (hyr) aS | ar OCR" by, 


A linear combination of the Schroedinger equations for ¢,, and ¢, 
(cf. § 4.1) gives easily 


k'—k ca sf , ’ r) t ’ 
; | r3V(r) ar | SV (r' \bin( yt’ )bn (Ky 0’) dr’ +-O(e'—k)?, 
0 0 


In the cases considered in this paper, SV is large only in a region near 
the origin where ¢,, and ¢, have no zero. An upper limit for the 
absolute magnitude of O(k’—k) is then obtained by replacing the limit 
of integration r by infinity. The second integral is then equal to 
—3k sin {ny (k)} (ef. §4.1), and 


O(k'—k) ~ —(k'—H) sin {n (4)} IL. r3V dr. 
0 


The maximum of |k’—k|, given by the maximum momentum ky, 
is a fraction of atomic unit in monovalent metals. This development 


is therefore justified, and O(k'—k) may be neglected, provided le roV dr 
is not too large, 0 
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For very large values of | k’—k | on the other hand, one of the functions 
dm OV $, is rapidly oscillating and may be replaced by a sinus function. 
Application of a lemma by Riemann and Lebesgue (Whittaker and Watson 
1950) shows that w,,,, is of the order of 


tah SV dr, eas) SV dr or (kk!) |" 8V dr, 
0 0 0 


according respectively to whether & tends to infinity, k’ remaining finite ; 
k’ tends to infinity, k remaining finite or kh, k’ and |k—k’| tend all 
three to infinity. 
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SUMMARY 


The abrasion properties of several diamonds have been investigated 
using a multiple beam interferometric technique. The cube, dodecahedron 
and octahedron planes were abraded in various directions using a special 
micro abrasion tester. Abrasion only occurred in directions parallel to 
crystallographic axes ; that is in four directions on cube planes and in 
two on dodecahedron planes, whereas no abrasion at all was produced on 
either natural or cleavage octahedron planes. Not all these directions were 
equally favourable; although it was found that an exact orientation reduced 
the vector effect. Taking the most favourable direction in each case the 
dodecahedron surfaces were easier to abrade than the cube surfaces. 

In addition the interferograms show that the abrading wheel moved in 
‘both lateral and longitudinal directions whilst abrading and at the same 
time probably suffered elastic deformation. It can be seen that the sides 
of an abrasion were highly polished and that there was no large scale 
deformation of the surrounding surface. 


$1. IyrRopDUCTION 


Diamonp polishers have known for several centuries that only certain 
faces of diamond can be polished, and that not all directions on these 
faces are equally favourable for grinding. However, no systematic study 
was undertaken until 1920 when Tolkowsky investigated the rate of 
grinding in different directions, and compared the rate of grinding for 
different crystallographic faces. Further investigations and interpre-’ 
tations have been advanced by Stott (1931), Kraus and Slawson (1939), 
Whittaker and Slawson (1946), Slawson and Kohn (1951). These workers 
concluded that polishing occurs most easily in directions parallel to 
crystallographic axes. Moreover, on a given crystallographic face it is 
equally easy to produce abrasion in each of these directions. For example, 
a cube face, parallel to two crystallographic axes, has four equally 
favourable directions for grinding. They found that the optimum 
grinding directions on a dodecahedron face were easier to abrade than the 
corresponding directions on a cube face: there was also a gradual change 


in hardness in moving from one crystallographic face to another, except 
a eee 
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that the octahedron face, equally inclined to all three crystallographic 
axes, could not be abraded in any direction. However, according to 
some diamond polishers (Grodzinski and Stern 1949 a, Slawson 1946) 
there is a pronounced vector effect of abrasion hardness; on a cube face 
for example, abrasion will only occur in one direction. 

The above results were obtained from the polishing of entire crystallo- 
grephic faces; the weight of diamond removed in a given time being 
determined. The disadvantage of this method is that any variations in 
abrasion hardness due to the presence of submicroscopic twins (Sutton 
1928) will not be revealed. By using a micro-abrasion tester (Grodzinski 
and Stern 1949 b) for the present investigation, only a small part of a 
surface was abraded at a time, and since the length of the abrasion made 
was only of the order of 10-?.cm several abrasions could be made on each 
face (the average area of a face being 4sq.mm). Further, any submicro- 
scopic twins would be revealed and it would therefore be possible to decide 
whether there were in fact fundamental variations in hardness properties 
between diamonds from different sources. 


§2. EXPERIMENTAL TECHNIQUE 


The Grodzinski micro-abrasion tester consists primarily of a double 
conical grinding wheel of 1 inch diameter with included angle of 110°. 
It was charged with diamond powder of 0—2 micron size, and rotated at 
approximately 10 000r.p.m. The diamond to be tested was so mounted 
on a balance arm that it was held against the abrading wheel by a force 
of 75g wt. The time of each test run was 15sec. The abrasions were 
made by W. Stern of the Diamond Research Department of the Industrial 
Distributors (Sales) Limited, and the abraded surfaces were studied using 
multiple beam interference techniques (Tolansky 1948). The depth of 
each abrasion was determined from the number of Fizeau type fringes 
contained within each outline, but when the depths were less than a light 
wave, fringes of equal chromatic order were used (Tolansky 1948). 
Examples of the two forms of fringes are shown in figs. | and 2 (figs. 1-8, 
Plate LX XII). Figure 1 (magnification x 44) shows Fizeau type fringes on 
the cube face of a glassy diamond, with arrows indicating the direction 
of abrasion, and fig. 2 (x 220) shows the fringes of equal chromatic order 
across the most shallow abrasion on this face. The length and width of 
each abrasion were determined from photomicrographs, but for reasons to 
be discussed later, the depth measurement was considered to give the more 
accurate estimate of the abrasion hardness; the latter being finally 
expressed as the depth per revolution of the abrading wheel. 

Two series of tests were made : in the first, two diamonds were tested, a 
glassy and an alluvial, both originating from South Africa. Cube and 
dodecahedron faces of both stones were abraded, several directions on each 
face being tested. The natural and cleavage octahedron faces of the 
glassy were also tested. These faces were first orientated visually from an 
inspection of the natural features, but in a second series the surfaces were 
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orientated by x-ray diffraction. In this series cube faces of four diamonds 
were abraded, including the glassy and alluvial stones of the first series, 
together with two stones from the Belgian Congo (ref. No. BCK No. 27a 
and BCK No. 71a). 

The hardness of the diamond had a most noticeable effect on the cast- 
iron wheel, which became very blunted with use. Also the amount of 
diamond powder on the wheel decreased with use, but although the wheel 
was recharged with powder after an average of 20 test runs, the wheel was 
not reground until it had become excessively blunt. This had a progressive 
effect on succeeding abrasions, for example, on the cube face of the diamond 
BCK No. 71a, four runs were made in the same direction, and the depth 
decreased steadily with successive runs. However, the variations in 
abrasion hardness between different faces, and between different directions 
on the same face, were always greater than the variation due to this lack 
of uniformity in the abrading conditions. Therefore, although no 
strictly quantitative assessment can be drawn from the experiments, a 
good qualitative estimate of abrasion hardness variations cen be obtained. 


§3. EXPERIMENTAL RESULTS 


With respect to the cube and dodecahedron faces of the glassy and 
alluvial stones, it was found that abrasion occurred only when the direction 
of grinding was parallel to a crystallographic axis, and when the directions 
of least resistance to abrasion were compared, it was found that the 
dodecahedron faces were some 25°% easier to abrade than the cube faces. 
These results are in qualitative accord with the results of Tolkowsky and 
Slawson (loc. cit.). However, an interesting feature emerged when the 
abrasion hardness values for the most favourable directions on the two 
stones were compared. The dodecahedron face of the glassy was easier 
to abrade than that of the alluvial yet on the other hand, the cube face of 
the glassy was more resistant than that of the alluvial. This fact makes 
suspect any general statement about the variation of abrasion hardness 
between different stones unless similar crystallographic faces are being 
compared. No abrasion mark was produced at all on either the natural or 
the cleavage octahedron faces of the glassy stone, in agreement with the 
traditional view that octahedron faces cannot be polished in any direction. 

With respect to the cube faces, it was found that there was a marked 
difference in the apparent abrasion hardness values between the faces that 
had been set by visual examination, and those that had been set by x-ray 
diffraction methods, the effect being more pronounced on the alluvial stone. 
Table 1 below gives the average hardness values expressed as depth of 
abrasion per revolution for the alluvial stone. 

It is seen that although abrasion occurred in all directions parallel to the 
two crystallographic axes, on the visually orientated face there is one 
markedly favourable direction of grinding together with a pronounced 
vector effect, whereas on the face orientated by x-ray diffraction there are 
two equally favourable directions of grinding and the vector effect has 


of the Abrasion Hardness Properties of Diamond 1143 


been reduced. These results may partly resolve the discrepancy between 
theoretical prediction and familiar experiment, for it would appear that 
diamond polishers may seldom be abrading the true cube faces. The 
abrasions produced on the glassy and on the two Belgian Congo stones which 
were orientated by x-ray diffraction showed similar results, in that there 
were two directions at right angles to each other along which the abrasion 
hardness values were nearly equal, whereas in directions at 180° to them 
the resistance to abrasion was high. 


Table 1. Abrasion Hardness Expressed as Depth rev-1 


Orientated 
by X-ray 
diffraction 
em x 10? 


Visually 
Alluvial Abrasion direction orientated 


em x 10? 


Parallel to a crystallographic axis 
90° to (1) 
180° to (2) 
180° to (1) 


Table 2 shows the average abrasion hardness values in the most 
favourable directions for the cube faces of the four different stones. 


Table 2 


Diamond Hardness value 


BCK No. 27a 0:80 


Glassy 1-13 
Alluvial 1-34 
BCK No. 7la 1:39 


(The hardness value is expressed as the depth of abrasion per revolution, 
the unit being cm x 10~* rev—1.) 

The magnitude of the differences between these values, as much as 40% 
between the two stones from the Belgian Congo, indicate a fundamental 
variation in the hardness properties of diamond as distinct from variations 
due to the presence of sub-microscopic twins. 

The dodecahedron faces of the glassy and alluvial stones were orientated 
visually, and fig. 3 (x36) shows the Fizeau type fringes over the surface 
of the alluvial stone, with the arrows indicating the direction of abrasion. 
It was found that there was a distinct difference between the results 
obtained from the two stones. On the alluvial stone, abrasion occurred in 
both directions parallel to the crystallographic axis whereas on the glassy it 
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only occurred along one of these directions, although a small abrasion was 
effected at right angles to the axis. Moreover this favourable direction 
on the glassy was so easy to abrade that the time of abrasion had to be 
reduced to 5 sec. However, these differences may have been caused by 
an inexact orientation of the two faces. 

In addition to the hardness determinations, the interferograms of the 
abrasions yielded some information concerning the effect on the abrading 
wheel of the hard diamond surface, for the geometrical outlines of the 
abrasions revealed several anomalies. Firstly, the included angle of the 
minor axis was found to be 176°, the angle being approximately the same 
for all the abrasions, yet the included angle of the double conical wheel 
measured before and after the test runs, was only 110°. Then, when the 
depths of the abrasions were related to the corresponding lengths and widths, 
it was found that although there was a general correlation between the 
length and depth, there was no such correlation between the width and the 
depth, occasionally the widest abrasion was also the shallowest. Finally 
when the measured depth was compared with the depth computed from 
the length of the abrasion and the geometry of the cast iron wheel, it was 
found that the experimental depths were always less than the computed 
depths, by between 10 and 25% ; the greater factor applying to the shorter 
abrasions. 

Three causes were considered for these anomalies, 

(a) that elastic deformation occurred during the abrading test, 

(b) that the wheel moved laterally in its bearings during the test, and 

(c) that the wheel moved longitudinally during the test. 

Subsidiary experiments were therefore carried out, using a wheel with a 
controlled lateral movement. The resulting width of the abrasion was 
considerably increased, and the form of the abrasion showed that the wheel 
had slid over the surface, making several slight abrasion marks, until a 
groove was formed; once the wheel was in the groove, the lateral move- 
ment ceased. It was therefore thought that the excessive width of the 
abrasions was caused by a mechanical flexibility in the supports of the 
wheel which had resulted in a lateral movement over the surface of the 
diamond; the effect being more pronounced the more resistant the surface. 

To investigate the difference between the computed and experimental 
depths, abrasions were made on a glass surface, and it was found that the 
computed depth was twice as large as the experimental depth. Since the 
wheel could not have been elastically deformed while abrading the glass, 
it was concluded that there was a longitudinal motion of the wheel which 
resulted in an increased abrasion length. It is probable that the uneven 
outline at one end of the abrasions which is shown by the interferograms is 
due to this longitudinal motion of the wheel. 

Thus the mechanical movements of the wheel could have caused the 
dimensional anomalies of the abrasions ; the longitudinal motion affecting 
the length, and the lateral motion affecting the width ; the latter motion 
also tending to increase the included angle of the minor axis. Yet, owing 
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_ to the large discrepancy between the included angle of the cast iron wheel 
and that of the minor axis of the abrasion, the possibility of elastic 
deformation of the wheel must not be completely disregarded. 

The interferograms of the abrasions also enable some conclusions to be 
drawn concerning the mechanism of abrasion. The figures show that the 
fringes within all the abrasions were extremely sharp and smooth, 
indicating that the surfaces within the abrasions were all highly polished, 
even though there is an anisotropy in the hardness property. When cube 
and dodecahedron faces were ground in directions at 45° and 90° respectively 
to the axis, the resulting surfaces were so uneven that no fringes were 
obtained, thus confirming earlier suggestions (Finch 1937, Bowden and 
Tabor 1950) that diamond is ground solely by the removal of. its 
constituent particles. 
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ABSTRACT 


A thin magnetic lens beta-ray spectrometer has been modified to 
enable the beta-gamma angular correlation of 17°Tm to be measured as 
a function of beta-ray energy. Assuming that a single matrix element is 
operative in the low energy beta-ray transition of 1°Tm, the results give 
a good agreement with the theoretical predictions for the matrix element 
fox r occurring in the first forbidden tensor interaction but appear to 
rule out the corresponding foxr of the first forbidden axial vector . 
interaction. It is, however, possible that a suitable mixture of matrix 
elements might also fit the experimental data. 


$1. [INTRODUCTION 


A BETA-GAMMA angular correlation in the decay of 17Tm has been 
discovered by Novey (1950). Assuming the correlation function 
W(0)=1-+a cos? 6, he obtained for a the value —0-259-+-0-053. However, 
as he points out, this result is somewhat uncertain on an absolute basis 
because his experiment was not carried out in a vacuum, and scattering 
in air would be expected to reduce the correlation. 

The existence of the correlation implies that the low energy beta-ray 
spectrum of 1”°Tm (see fig. 1) cannot be an allowed transition, and a 
recent investigation into the shape of this spectrum by Richmond and 
Rose (1952) has indicated that the transition is first forbidden with axial 
vector or tensor interaction, the latter interaction being slightly favoured. 
By investigating the beta-gamma angular correlation as a function of 
beta energy it has been possible to draw more definite conclusions 
regarding the type of interaction obeyed in the low energy beta-ray 
transition. 

§2. EXPERIMENTAL DETAILS 


A thin magnetic lens beta-ray spectrometer was used to focus beta- 
particles of specific energies while gamma-rays from the same source 
were detected with a NaI crystal which could be rotated about the source. 
The source itself was mounted at the centre of a source-holder which was 
essentially a hemisphere of aluminium, so that the gamma-rays penetrated 
equal thicknesses of aluminium at all angles of detection. 


SS Ea Ee eee 
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The spectrometer was set at a momentum resolution of 9% to give 
a relatively large transmission. At this setting the spectrometer baffles 
limited the detection of beta-particles to those leaving the source at 
angles between 3° and 8° to the spectrometer axis. 

The beta-particles were counted by an anthracene crystal which was 
mounted on the vacuum side of a glass disc sealing the exit hole in the 
end-plate, and the light pulses were piped off by a Perspex rod 18 cm long 
to a photomultiplier (E.M.I. type 6262). The thallium-activated Nal 
crystal which detected the gamma-rays was about 4:5 cm? in cross-section 
and 1 cm thick and subtended a solid angle of 0-08 steradian at the source. 
The crystal was sheathed in a lead cylinder 1-5 cm thick to screen it from 
all radiation other than that coming directly from the source. A 25 em 
Perspex rod acted as a light-guide between the Nal crystal and a second 
photomultiplier, the whole being capable of rotation through 90° on either 
side of the spectrometer axis. 


Fig. 1 


7° Tin ENS 


ed 0-970 Mev. 


0:084 MN 26 
0 XG>) 


Decay scheme of 1”°T'm. 


The Perspex rods acted not only as light-guides but served also to 
remove the multipliers to some distance from the spectrometer lens coil. 
It was then possible to provide sufficient magnetic shielding to reduce 
the effect of the magnetic field of the coil on the multipliers to an 
inappreciable amount over the range of spectrometer currents used in the 
experiment without resorting to steel end-plates and a steel tank enclosure 
for the coil as was necessary in the beta-gamma angular correlation 
_ investigations of Stevenson and Deutsch (1951). 

Coincidences were recorded between focused beta-particles and gamma- 
rays with the gamma detector at 90°, 180° or 270° to the spectrometer 
axis. To minimize errors due to changing efficiencies the angle of the 
gamma, counter was changed and the observations repeated a considerable 
number of times for each beta-ray energy selected. The coincidence rates 
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were normalized with respect to the single counting rates in the gamma 
counter, The resolving time of the coincidence circuit used was 
8 x 10-8 sec. 


§3. EXPERIMENTAL RESULTS 


To ensure that the experimental arrangement did not suggest a 
correlation where none existed, tests were first carried out on Cs 
and Au, both of which are known to have isotropic correlations. 
(Stevenson and Deutsch 1951, Frankel 1950). 

134(\5 was deposited as an aqueous solution on to 0-25 mg cem™ 
aluminium foil and evaporated down to form a source of thickness about 
0-5mgem-2. The ™8Au source was electro-deposited on to a similar 
source backing and was only about 30yngem™ thick. Both sources 
covered about 1 cm? in area. 

No correlation was observed in the decay of either source. The results 
obtained are given in table 1. 


Table 1 
| 
Energy of a (assuming 
Isotope beta-rays W (0)=1-+a cos? 6) 
(Mev) 
1341, 0-300 —0-006 0-024 
. 0-489 +0-021 +0-032 
198Ay 0-614 +0-013+0-035 


From some ‘Specpure’ Tm,O, of high specific activity (4 me mg!) a 
source of thickness 1 mg em~? and area 1 cm? was prepared on 0-25 mg em-? 
aluminium foil. A 1 gem brass absorber was placed in front of the Nal 
crystal. This transmitted about half of the gamma-rays but absorbed 
almost all the K x-radiation of 17°Yb, which was expected to be isotropic 
with respect to the beta-rays (Novey 1950). 

The 84 kev gamma-rays are to a large extent internally converted, 
leaving only a few per cent for detection. Because of this relatively few 
genuine coincidences were observed, and in order to obtain the statistics 
shown in fig. 2 it was necessary to count for periods of 40 hours or more 
at each point on the beta-ray spectrum. 

The angular correlation function W(@) is of the form 1-+-a cos? @ for all 
first forbidden beta transitions followed by a gamma-ray (Falkoff and 
Uhlenbeck 1950). The values obtained for the differential correlation 
coefficient a at various points on the beta-ray spectrum, after correction 
both for the finite solid angles of acceptance of the detectors and for the 
spectrometer resolution, are given in fig. 2. The errors indicated are the 
probable errors based on the number of counts. 

Assuming an allowed spectrum shape (Richmond and Rose 1952) for 
the low energy beta transition an integral correlation coefficient 
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a——0-19--0-04 was calculated on the basis of the experimentally 
observed differential coefficients of fig. 2. This appears at first glance to 
be rather lower than Novey’s reported value, but it must be borne in 
mind that in his experiment he used a 3 mg cm~? mica window over the 
beta detector to prevent counting of conversion electrons. Consequently 
his value is not a true integral correlation coefficient but represents a 
correlation integrated over beta energies between about 150 kev and the 
end point. Such a procedure, when applied to the results of fig. 2, yields 
a value for a of —0-24+-0-04, which is in good agreement with that of 
Novey. 


a 
-- 
=- 
-- 


0-7 0-8 0-9 1-0 


Beta-ray energy (MeV) 


0-6 


Beta-gamma angular correlation of 1Tm. The solid line denotes the 
theoretical correlation (4=0) for the first forbidden matrix element fo xr 
occurring in the axial vector interaction, the dotted line that for the first 
forbidden matrix element {Bo xr of the tensor interaction. 


$4, Discussion 


Since 17Yb is an even-even nucleus its ground state should have 
spin zero and even parity. Measurements of conversion electron 
intensities (Grant and Richmond 1949) and of the lifetime of the 
84 kev excited level (Bell and Graham 1950) have shown that the 
84 kev transition is electric quadrupole and, therefore, the excited level 
has spin 2 and even parity. The shape of the low energy beta-ray 
spectrum (Richmond and Rose 1952) indicates that the transition is 
first forbidden and requires 17°T'm to have spin | and odd parity. 

With this assignment of angular momenta the only first forbidden 
nuclear matrix elements which can give the negative coefficient a are 
foxr and {foxr occurring in the axial vector and tensor interactions 
respectively. The theoretical correlations (Falkoff and Uhlenbeck 1950) 
to be expected for these two matrix elements are plotted in fig. 2. It 
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can be seen that whereas they yield the same correlation for the high 
energy beta-particles, they differ considerably at intermediate energies. 
The experimental points agree well with the theory for the tensor inter- 
action and seem to eliminate the axial vector as a possibility. 

The theory of Falkoff and Uhlenbeck (1950) is strictly correct only 
for the case Z=0, and little is known about the precise influence of Z. 
Although the effect of Z is generally quite small (Fuchs and Lennox 
1950) it can substantially lower the correlation in absolute value from 
that computed in the Z=0 approximation for cases where matrix elements 
like foxr are involved (Stevenson and Deutsch 1951). This seems to 
rule out completely the possibility of the axial vector being the true 
interaction in the low energy beta transition of 1%°T'm but might explain 
why the experimental points fall somewhat below the theoretical curve 
for the tensor interaction. There is also the possibility that a suitable 
mixture of first forbidden matrix elements might fit the experimental 
results, but no theoretical calculations of explicit correlation functions 
assuming such mixtures are at present available for comparison. 
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SUMMARY 


A theory of work-hardening is put forward ; the salient points are as 
follows. There are two types of hardening ; the normal rapid hardening 
of cubic metals and the slower hardening characteristic of hexagonal 
crystals and cubic metals under conditions of easy glide. A quantitative 
theory of the rapid hardening is given. This does not occur unless dis- 
locations are retained in the slip planes. The retained dislocations are 
supposed to occur in groups of about 1000 at the ends of the slip-bands. 
The strains round these groups are responsible for hardening as in Taylor’s 
(1934) theory. The obstacles which stop slip are sessile dislocations, 
and the formation of sessile dislocations stabilizes the groups of dis- 
locations. It is shown that the formation of slip lines must be a dynamic 
process, and on this assumption a parabolic stress-strain curve is deduced. 
Some discussion is given of the formation of vacancies by moving dis- 
locations, and it is shown that this can account for the clustering of slip 
lines observed under the electron microscope. The theory agrees 
qualitatively with observations on the amount of retained energy in a 
cold-worked metal. 


§1. InTRoDUCTION 


THE aim of this paper is to give a theory of work-hardening of metal 
crystals, based on the properties of dislocations. The basic assumption. 
made is that a crystal contains a network of dislocations, and that a 
dislocation in an otherwise perfect crystal moves when the crystal is 
subjected to a small stress. We shall attempt from these assumptions 
to build up a quantitative theory for close-packed cubic crystals, and shall 
make a few qualitative remarks about hexagonal crystals. 

The first theory of work-hardening based on dislocations is that due 
to Taylor (1934). Qualitatively this may be expressed as follows. The 
stress round an edge dislocation is of the form 


Pee (CDIEVEL 0 See eat cy eke a (L) 


* Communicated by the Author. 

+ It is not suggested that it is inherently impossible to prepare crystals 
without dislocations ; the dislocations in real crystals are a result of the method 
of growth. Some of the factors involved have been summarized by Frank 
(1952). A recent note by Herring and Galt (1952) shows that very thin fibres 
formed on tin crystals have yield strengths approaching the values expected. 
for a material without dislocations. 
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where (is the shear modulus, b the Burgers vector, 7 the distance from the 
dislocation line and f(@) a function of angle. — If, then, a cold-worked 
solid contains a random array of dislocations of both signs at a mean — 
distance J from each other, there will be in the material a stress o,, varying 
with position in some roughly sinusoidal way, of which the magnitude 
is of order 


o,=6b/L +. 4 on 


If each dislocation of the array is in some way ‘ stuck ’, any additional — 
dislocation introduced into the lattice cannot move through distances 
large compared with / unless the applied stress o exceeds o;. ; is then 
identified with the flow stress of the crystal at any point of the stress— 
strain curve. It is supposed in Taylor’s model that dislocations are 
constantly being created while the substance is being deformed, and 
that they move a distance L before they are stuck. The plastic strain « 
is thus 

e=Db/P. . .. - ee 


Eliminating / between (2) and (3), we find 


a= | zee vO aaa 


To obtain agreement with experiment for the stress-strain curve of 
aluminium, Taylor had to assume that Z is about 10-4cm, and he 
deduced that, after a strain of order unity, 1 ~ 3x 10-& cm. The concept 
of an oscillating internal field o; has been adapted by Mott and Nabarro 
(1948) and Nabarro (1952) to the explanation of age-hardening and the 
high yield-stress of alloys ; while an explanation of yield points in iron 
and similar materials has been given by Cottrell (1948 ; see also Cottrell 
and Bilby 1949) in terms of the segregation of impurities actually in or 
very near the dislocations. 

The explanation of strain hardening which we shall give is an extension 
of Taylor’s ; we suppose that the hardness is due to an internal field 
due to dislocations generated in the process of plastic flow. The actual 
distribution, however, is very different from that of Taylor, dislocations 
of like sign being grouped together, some hundreds being in each group. 
It shares with Taylor’s theory, however, the feature that the hardening 
is, to a first approximation,* independent of the direction in which slip 
has occurred ; in a cubic crystal dislocations generated by slip along 
one set of (111) planes harden the crystal for slip along another set of 
(111) planes. 

To complete Taylor’s theory, we have to ask : 

(a) How are the dislocations generated during slip? For this a 


satisfactory mechanism has been proposed by Frank and Read (1950), and 
will be discussed in the next section. 
eee ee 
* This point has been examined further by Stroh (1952) who i 
. ; é Bass : t 
be published shortly finds a small deveeienee ee ale 
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(0) How do the dislocations become ‘ stuck’ after moving a certain 
distance, what limits slip on a given plane, and why do not the dislocations 
move back when the stress is reversed? Here we propose different 
mechanisms for the edge and screw components of a dislocation loop. 
The edge components in a cubic close-packed crystal become stuck by 
combining with edge dislocations on other octahedral planes, as proposed 
by Lomer (1951) and Cottrell (1952). A join of two such dislocations 
produces what we call a sessile dislocation. This cannot happen in 
hexagonal crystals, which is why they harden so much more slowly. 
The screw dislocations are probably stopped by a process connected with 
the generation of vacancies by moving dislocations (Seitz 1952). These 
vacancies also play a large role in the temperature-dependence of strain 
hardening and’ in the increase in electrical resistance due to strain 
hardening. 
Some of the experimental facts which have to be explained by any 
theory are: 
(a) The occurrence of slip lamellae. ; 
(6) The fact that at low temperatures, in aluminium at any rate, a 
slip lamella consists of slip on one plane through distances of the order 
1200-1800 A. As the temperature, however, is raised, the visible slip 
bands are found to consist of clusters of these elementary steps, about 
200 A apart and with roughly the same slip on each (Heidenreich and 
Shockley 1948, Brown 1949, 1951, 1952). 
(c) The very considerable temperature-dependence of hardening. 
(d) The form of the stress-strain curve. 
(e) The formation of deformation bands and asterisms in cold-worked 
cubic metals, and their absence in hexagonal crystals. 
(f) The phenomenon of easy glide (Andrade and Henderson 1951, 
Liicke and Lange 1952). 


§2. THe FRANK-READ SOURCE AND THE ORIGIN OF SLIP BANDS 


Several accounts have already appeared of the explanation given by 
Frank and Read of the origin of slip bands (Frank and Read 1950, 
Frank 1951, Mott 1951). The model is illustrated in fig. 1. It is based 
on the concept of a dislocation line EF lying in a slip plane, which at the 
_ points E, F moves on to other planes lying at an angle to it. We consider 
movement in which EF, rather than its continuation in the other planes, 
is displaced. 

Consider then what will happen to the line EF, when the material is 
stressed in such a way as to cause slip. If the line moves from the 
position shown in fig. 1 (a) to that shown in fig. 1 (6), the material over 
the area shown shaded in fig. 1 (b) slips relative to that underneath it 
The line will move in such a way as to extend this area. On the other 
hand the points E, F are locked. Successive forms of the line are shown 
in figs. 1 (a) to (e). It will be seen that ultimately a complete loop 
is generated, and that at the same time the line EF is formed again. 
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A line such as EF is thus capable of generating an infinite series of loops, 
when a large enough stress is applied. Moreover, it will be shown later 
in this section that the stress is of the order of that for which yield of 


actual metal crystals occurs. 


Fig. 1 
E E = 
, a 2. 
(a) (b) (Cc) 
E 
F 


! (d) (e) 


Successive stages in the formation of a dislocation ring from a Frank—Read 
source. 


For close packed cubic crystals—to which much of this article will be 
devoted—the model becomes more plausible when we consider the Burgers 
vectors possible for slip in the (111) planes. Such a plane is shown in — 
fig. 2(a). The possible values of the Burgers vector are OX,, OX, and — 
OX,;. It is clear that the displacement when two dislocations with 
vectors OX,, OX; move across the plane is the same as if a dislocation 
with vector OX, moves across it. From this it follows that it is geo- 
metrically possible for three dislocations to meet in a point, as in fig. 2 (5), 
the Burgers vectors of the three being respectively OX,, OX,, OXs. 
Further, it is obviously not necessary that the three lines should all lie 
in the same slip plane. 


(b) (c) 


(a) A close-packed plane, showing the Burgers vector for complete and half 
dislocations. (b) Join of three dislocations. (c) Network of dislocations 
(from Mott 1951). 
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A possible arrangement of dislocations in a cubic close-packed. solid, 
and probably the most stable in view of the line energy of dislocations, 
is therefore a network such as that shown in fig. 2 (c). This is the model 
that we shall adopt throughout this report for the arrangement of 
dislocations in an annealed metal of close-packed structure; it will 
probably replace the old model of a ‘mosaic crystal’ so far as the 
interpretation of intensity measurements in the reflection of x-rays is 
concerned. With this model, any element such as AB can act as a 
Frank—Read source. 

There are certain further points that we may note about this model. 
These are : 

Extended Dislocations. As first pointed out by Heidenreich and 
Shockley (1948), a dislocation in a close-packed plane with Burgers 
vector say OX, in fig. 2 (a), can and will split into ‘ half-dislocations ’ 
with Burgers vectors OY, YX,. These two half dislocations repel each 
other; it is easily seen that the elastic energy, proportional to the 
square of the Burgers vector, is diminished when they separate in the 
ratio 

20 Y7/0 Xi" 
and hence in the ratio #:1. On the other hand, when they separate 
they create a stacking fault between them, the energy of which increases 
with the distance between them. According to Heidenreich and 
Shockley, they will be in equilibrium when separated by about 2 x 1077 cm. 

As a result of this dissociation, the dislocations in close-packed planes 
will have low energy, and the elements such as AB (fig. 2 (c)) in cubic 
close-packed structures will mainly be in close-packed planes. This is 
believed to be the reason why glide in close-packed structures shows a 
marked preference for close-packed planes, while for materials with 
some other structures (e.g. body-centred metals, silver halides), various 
planes show glide. 

Yield Points. Ina pure single crystal flow will begin when the applied 
stress is great enough to initiate the process illustrated in fig. 1, namely 
the generation of loops at the most favourably situated Frank—Read 
source. The force resisting the generation of loops is due to the energy 
per unit length, or tension of a dislocation. If o is the resolved part of 
the stress about a line lying in the slip plane parallel to the Burgers 
vector, the force exerted on unit length of the dislocation is ob. If W 
is the energy of the line per unit length, the resultant force in the opposite 
direction is W/p, where p is the radius of curvature. The maximum value 
which 1/p can reach as the loop expands is 1/2/, where / is the length EF 
of the source. Thus the source will generate loops if 


o> W/2lb. 
Putting in estimated numerical values for W for isotropic materials, 
Frank (1950) finds that in practical cases W ~ 2G6?, so that the condition 


for yield is 
COU mee ta iee ge 5, tee te cre = (9) 
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For the best single crystal of, for example, silver (Andrade and Henderson 
1951), yield points o of the order 50 g/mm? are observed, giving values 
of 1 of about 3X10-3cm; for aluminium Rosi.and Mathewson (1951) 
observed yield points of the order 200 g/mm?* give 1~ 3x 10-* em: 
We emphasize that the values of / deduced are for the largest values of 
the length AB in the dislocation network discussed ; the average value 
may be smaller. No theory has, as yet, been given which enables us to 
estimate the mean value of / that will result from a given method of 
preparing a crystal. . 

There are two further effects to be considered. Firstly, there will 
exist a random internal stress due to other dislocations or impurities 
which changes sign in a distance A smaller than the estimated values of 
1; this effectively introduces a frictional force, say ob per unit length 
resisting the motion of a dislocation ; a dislocation initially at rest will 
not move unless the applied stress exceeds op, and a Frank—Read source 
will not generate dislocations unless 


Oy} Goll ao. 24 re 


Secondly we have to consider the locking effect of impurities which 
have by diffusion moved into the immediate neighbourhood of a dis- 
location. This effect and its relation to yield points have been considered 
in a number of papers by Cottrell and his co-workers (Cottrell 1948, 
Cottrell and Bilby 1949). For our purpose we need to emphasize the 
following effects of such locking : 

(a) The stress (og) required to initiate a slip band is now simply 
that required to make a dislocation break away from its impurities, 
provided that og is greater than o,+Gb/l. Once the dislocation has 
broken away, however, the resistance to further slip is still given by (6). 

(6) Since ogo, will be the same for all loops and is thus independent of J, 
the conditions are favourable for sharp yield points. 

(C) Goo, is very dependent on temperature (Cottrell and Bilby 1949), 
while og is very little dependent on temperature (Mott and Nabarro 1948). 

It is of some interest to consider the effect of a network such as that 
considered here on the apparent elastic modulus. If the dislocations 
are locked in the Cottrell sense, there should be no movement for stresses 
less than ogy. Up to this stress the material should be perfectly elastic 
with the theoretical shear modulus. In the absence of Cottrell locking, 
however, a stress greater than o, but less than o +Gb/l will lead to a 
reversible movement of any dislocations present in the material. If the 
material is highly purified so that o) is very small, this will lead to an 
apparent change in the elastic modulus. We shall now estimate the 
magnitude of this for the network illustrated in fig. 2 (c). If the stress 


applied is o in the plane considered, the dislocation takes up the form 
of a circle of radius p given by 


ob=W]/p. 
Substituting for W we find 
p=2bGc. 
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The area swept out by the dislocation (the shaded area in fig. 1 (b)) is 
thus p°(@—4 sin 20), where sin@=11/p. For small 6 this reduces to 


3/12p=[o/24b@. 


If there are N Frank—Read sources per unit volume, it follows that the 
shear produced by a stress o is 


Nic 2a ge piae. Tor es AS ear 


Since WV/3 is probably slightly less than unity, it follows that the correction 
to the elastic modulus is of the order of at most a few per cent, unless o 
approaches the stress Gb/l. 

It follows that the stress-strain curve of a pure substance having a 
number of Frank—Read sources shows a (reversible) strain slightly in excess 
of that given by the ideal elastic modulus. Under alternating stresses 
there will be a slight dissipation of energy, of amount proportional to the 
temperature, due to a mechanism first put forward by Eshelby (1949). 
If, however, there is any frictional resistance to the motion of dislocations, 
one will obtain a hysteresis loop (cf. Nabarro 1951). Such a resistance 
will certainly exist in alloys for the reasons already stated. 


Fig. 3 
LAE Ae 
if 
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Glide as the result of stress. 


We may now sum up the elementary consequences which follow from 
the assumption that the crystal contains a number of Frank—Read 
sources, or a network of the type shown in fig. 2 (c). If the crystal is 
submitted to a uniform shear stress, then as the stress is increased a value 
will be reached at which one of the sources, that for which / is largest, 
begins to generate loops. The process will then continue indefinitely ; 
nothing introduced into the model so far will stop it. Glide as shown in 
fig. 3 is thus to be expected, until owing to the change of shape of the 
erystal the stress is no longer uniform, when other sources would be brought 
into play. If by some mechanism slip on the first plane was stopped, 
the source with the next largest value of / would be brought into play. 
The next task of a theory is to explain why slip ever stops on a given 
plane. In particular we have to explain the distance of c. 2000 A 
shown by the work of Heidenreich and Shockley and of Brown to be 
characteristic of slip on each active slip band. 
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§3. Tor Prtinc-up OF DISLOCATIONS AGAINST A BARRIER 


The model outlined in the last section provides for a limited number 
of sources, say 101°—10% per cubic em in an annealed crystal each of 
these sources is capable of producing an unlimited number of dislocation 
loops when a large enough stress is applied. If in the material there 
are no barriers to the movement of dislocations, these loops will spread 
to the boundaries of the crystal where they will disappear, causing glide 
as in fig. 3. We shall have, however, to consider what happens when the 
material contains some internal surface, which acts as a barrier against 
which dislocations may pile up. Such a barrier may be, for instance, a 
grain boundary ; later in this article we shall consider the influence of 
trapped or ‘ sessile’ dislocations and of deformation bands. 

To fix our ideas, we shall consider a row of dislocations, either edge or 
screw, moving in a single slip plane, piled up against a barrier under the 
action of a stress o, as illustrated in fig. 4. This important case has been 
treated in detail by Eshelby, Frank and Nabarro (1951). The following 
results are obtained : 

(a) If there are n dislocations present, the force on the leading dis- 
location is n times that ona single dislocation, namely nob per unit length. 


Fig. 4 


BARRIER 
SOURCE 


Pile-up of dislocations against a barrier. 


(b) The length L between the first dislocation and the barrier is given by 
L=kn@blag,, 1s 3. ee ee 


where k is unity for a screw, 1—v for an edge dislocation, G is the shear 
modulus and v is Poisson’s ratio. 
(c) The distance between the two leading dislocations is 


d=1-84 dy/n, 
and that between the last two of order dyn? where 
dy=Gb/2ck. 
(d) The shear stress on the glide plane at a distance x ahead of the 
leading dislocation is practically the same as that due to this dislocation 


by itself in the region where a<d. At distances large compared with d 
but small compared with L the stress is approximately 


of{1+(L/a)t2}, 
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The stress concentration is thus large, the dislocations producing the 
same stress concentration as does a freely slipping crack of the same 
length. 

(e) The total slip produced by the sequence is the same as would be 
produced if a single dislocation of strength nb moved a distance 3L/4. 

Some numerical values will be of interest here. First of a well, 
estimate the magnitude of the distance dj, which is equal to the equili- 
brium distance between two dislocations of the same sign, one of which 
is held up by a barrier and the other forced towards it by a stress o. 
For a metal such as Cd or Al, we may take G ~ 2x 10" dynes/cm? ; 
for the yield point of a pure metal, o ~ 100 g/mm? ~ 107 dynes/cm ; 
this gives dj) ~ 3108) ~5x10-5 cm. For a crystal that has been 
heavily deformed the flow stress will be from 10 to 100 times greater and 
the magnitude of d, is 10-° or 10-*cm. Any theory which seeks to ascribe 
work-hardening to the interaction between dislocations fairly uniformly 
distributed (e.g. Taylor 1934) must assume that the distance between 
them is of this order, and thus considerably less than the distance between 
observed slip bands. One has to assume, either that in a cold-worked 
metal the space between the slip bands is filled up with dislocations, or 
that in the slip bands there are clusters of piled-up dislocations, formed 
as described above. ach of these clusters will behave at large distances 
like a dislocation of large Burgers vector. In the latter sections of this 
report we shall develop particularly the second of these hypotheses, 
though we believe that there will be some formation of dislocations 
between the slip bands. 

A second point of interest is the high stress concentration that exists 
near the limit of a row of piled up dislocations. Frank (1951) has suggest- 
ed that, when this occurs at a grain boundary, the stress may be great 
enough to generate loops. This point we believe to be important, and 
it will be considered in further detail in § 6. Here we wish to emphasize 
that such a process is likely to be temperature-dependent ; temperature 
can assist the formation of a loop in a region of very high stress. 

When the surfaces of aluminium crystals which have been deformed 
at low temperatures are examined with the electron microscope, single 

_steps of the order of 1000-2000 A in height are observed. At higher 
temperatures the steps are clustered ; a hypothesis to account for this 
clustering will be given in §8. At low temperatures, however, this 
suggests that about 1000 dislocations are produced by each Frank—Read 
source. The edge dislocations of one sign producing the visible slip band 
have come out of the crystal; those of the opposite sign, at the opposite 

- side of the ring, must in general be piled up against some barrier. We 
shall now estimate the distance from the surface that this barrier must be. 
Let us assume for the moment, then, that when each ring expands the 
screw elements go to the boundary of the crystal or at any rate go a 
large distance compared with the edge components, but that the latter 
are stopped at a barrier after moving a distance L. We have to ask how 
big this distance ZL must be in order to account for a value of n of c.1000. 
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We shall take the case of aluminium single crystals subjected to a 
shear stress of 500 g/mm? or 5x 10? dynes/em’. If G~3x 1011, and 
n ~ 1000, then formula (8) gives L=0-05 cm. The distance of the first 
barrier from the surface may be as much as twice this. Thus, for slip of 
2000 A to take place for the lowest stresses at which yield occurs, a slip 
distance L of nearly 1 mm is necessary. . 

We shall now discuss the dynamics of the formation of a slip line. 
From the theoretical point of view, it is not yet certain how much a 
moving dislocation in a perfect crystal is damped by the thermal vibration 
of the lattice (Leibfried 1949, 1950, Nabarro 1951). We believe, however, 
that the facts about the formation of slip bands show conclusively that 
dislocations must acquire a velocity not less than about half the speed 
of sound. For all the observations of slip bands (Brown 1951) show that 
the elementary slip band through c. 2000 A is formed completely or not 
at all, though the formation of other members of a cluster may occur 
after some time. The evidence suggests, then, that the elementary slip 
band is formed quickly, and not gradually as the deformation proceeds. 
While a crystal is being strained, the stress rising slowly, slip from a 
given source begins then presumably when o reaches the value of Gb/l 
for that source, and the band is formed before o has risen appreciably 
above this value. If the dislocations moved slowly, say with one tenth 
of the velocity of sound, this would be impossible ; as soon as, say, 4n 
dislocations were piled up at the end of the slip distance L, they would 
exert a stress AnG/27L at the source, preventing the generation of new 
dislocations. This would stop the source from acting long before the 
final value n was reached. If the motion of disclocations is heavily 
damped, we should expect n to increase gradually with o ; it could not 
reach its full value until o reached about twice the yield stress. 

If on the other hand dislocations can move with nearly the speed of 
sound, once a Frank—Read source is set into action, it will go on generating 
dislocations under its own momentum, the momentum generated in 
moving from the configuration of fig. 1 (e) carrying the line EF to the 
configuration of fig. 1 (c). Generation of loops will then go on until the 
stress of the piled-up dislocation stops it, the whole process being over in 
a few microseconds.* 

Moreover, if the process is rapid as described here, we can, assuming 
the presence of a barrier at a distance L from the source, account for 
the fact that the dislocations do not move back when the stress is taken 
off. ‘This can be explained in terms of the formation of sessile dislocations, 
which is likely to occur as soon as a piled-up group of more than about 
a a ee re ee eee 


* A somewhat similar suggestion has been put forward by J. C. Fisher, 
EK. W. Hart and R. H. Pry (to be published), These authors consider that 
the stress due to the still moving dislocations generated by the source will 
stop the source from acting after some hundred are formed, without their 
coming up against any barrier. I am indebted to Dr. Fisher for allowing me 
to see his manuscript prior to publication. 
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six dislocations is brought to rest, but will not occur while they are in 
motion. We have already stated that when two dislocations on two 
octahedral planes join, they form a sessile dislocation. A piled up group 
of n dislocations has a stress field round it. of order n@b/2ar. There are 
other sources at a mean distance from it of order J, and a stress of order 
Gb/l is required to make any one of them operate. Thus a piled-up 
group of six is sufficient to cause plastic flow round it, and to produce 
one or more edge dislocations on a different plane which will form a 
sessile element with one of the six. We suppose, then, that this is what 
happens, and a proportion of about one in six of the dislocations piled up 
at the end of a slip plane become sessile, thus preventing any appreciable 
recovery of strain when the stress is removed. 

We emphasize that this phenomenon cannot occur for hexagonal 
crystals, where there is only one set of slip planes. 

It is clear that in ordinary plastic flow, only a small proportion of the 
available sources are used. For supposing that the slip distance for edge 
and screw dislocations is L, then the strain « produced by the operation 
of NV sources per cm? is 

e=NLL'bn. 


With L~L’~107%cm, n~ 1000 and «~ 1, this gives N ~ 108, 
considerably less than the value 101°—10! which we believe to represent 
the total number. 

Evidence in favour of this conclusion is provided by the work of Brown 
and Honeycombe (1951), who have shown that the slip markings distant 
c.1 micron from each other are characteristic only of surfaces which have 
been mechanically polished, and that they do not appear on surfaces that 
have been electrolytically polished and etched. We believe, however, 
that the former markings are characteristic of the interior, and for the 
following reason. As first pointed out by Holloman (1952), dislocations 
which terminate on the surface of an otherwise perfect crystal will behave 
as though locked between the real locking point and its reflection in the 
surface. They will have, therefore, much longer effective values of 1. If 


l1~ 10-4 cm, 


there will be per unit area of surface c. 108 sources which terminate in 
the way shown ; and if a stress is applied which will operate any of the 
sources in the interior, it will operate also most of the surface ones. 
If the dislocations generated move a distance 10° cm, we see that a 
density of slip lines 10° per cm is to be expected on the surface. Under 
such conditions it is natural that each line can only have a small dis- 
placement. 

If the surface has been cold-worked or hardened in any other way, the 
situation is quite different ; sources on the surface do not move more 
easily than those on the interior, and if the hardened layer is thin compared 
with the slip distance (say 10cm), the markings shown should be 
characteristic of the interior. 
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§4. STRAIN-HARDENING OF CUBIC CLOSE-PACKED CRYSTALS 


We give in this section a theory, allied to that of G. I. Taylor, of the 
work-hardening of cubic close-packed crystals at low temperatures. 
The theory will be limited to one rather rapid hardening characteristic of 
impure crystals; the ‘easy glide’ sometimes obtained with pure 
crystals will be discussed in the next section. The assumptions that we 
shall make are : 

(a) There exist on all active slip planes barriers at a mean distance 
2L from each other ; the nature of the barriers will be discussed below. 

(b) The screw components of the dislocation loops move a distance L’ 
comparable with or greater than L ; the screws will be discussed in § 6. 

The hardening will be ascribed to the stress round the piled-up groups 
of n edge dislocations already described. If the surrounding material 
were elastic, the stress at a distance r from such a group would be 


nGb/2nr (n ~ 1000). 


However, at a distance greater than / (J ~ 10-4cm) the material is 
plastic, as it contains Frank—Read sources. As, however, the stress is 
extremely unhomogeneous, we do not believe that the dislocations from 
any source are likely to move a distance much greater than / before meeting 
a dislocation on another plane and becoming sessile. We therefore have 
to relieve a local strain «(r) (=nb/27r) by allowing dislocations to move 
from sources distant / apart through a distance J. It is clear that, at a 
distance r from the pile-up, each source must contribute n’(r) dislocations, 
where n’(r)=l/e(r)/b, to do this. Round each primary group of n which 
have moved a large distance L, we have other secondary groups of n’ 
which have moved a small distance /. Moreover, these produce a stress 
field o; at distance r from the primary group having ‘ wave-length’ 1 
and magnitude 
o ;=Gobn' /l=Ge(r)=Gbn]2z7. 


The magnitude is thus the same as if secondary slip had not occurred, 
If, then, a material has slipped so that slip bands are at a distance 
x apart, the mean distance between the primary groups of piled-up 
dislocations is (wL)1/, 
The mean strain o, in the material is then of wave-length J, the 
dimensions of the network, and of magnitude 


G == Zar (EL) Pe ee 
The strain, on the other hand, is given by 
a e~nb/x. Me ont Saal 
Eliminating a, we find 
Op 1S RO 
a=s/T 5 ee ee eee 


This formula is to be compared with Taylor’s formula (4). 
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This stress o; cannot, as in the original theory of Taylor, be identified 
with the flow stress at a strain «. The situation is to be pictured as in 
fig. 5, in which the stress acting on a dislocation is plotted against position. 
Before an external stress is applied, the situation will be as in fig. 5 (a), 
and the loop of a Frank—Read source will be at or near a point such as D, 
where the stress vanishes. As the applied stress is increased, the stress 
within the material will change to that shown in fig. 5 (b), the dislocation 
remaining at D. Buta slip band will not be formed until the configuration 
of fig. 5 (c) is reached, the stress at the maximum reaching the value 
Ac=Gb/l required to operate the source. As soon as this is reached, the 
(dynamic) generation of a slip band occurs ; if there is a barrier at B, the 
dislocations generated will all pile up in its immediate neighbourhood, the 
mean stress acting on them beinga,+ do. But the number n of dislocations 
generated depends on Jo only ; for n dislocations near B give at the 
source a stress field 


Showing internal stress on a material ; dislocations in fig. 5 (c) are shown as 1. 


so that the source will stop generating when 


nGb/2rL= Ao ; 
when this is so, 
L=nG@b/27 Ao, 


a formula differing by a factor 2 from (8). 
Substituting Gb/l for 4c, we have 


nb/L=2nb/l, 
whence formula (11) for o; takes the form 


ge (es ene ea) 
ces 


We may take the flow stress to be 
Treble mt Panes aie >.) (13) 
except for small strains the second term is small compared with the first. 
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Several points of interest may be noted about formula (12). The first 
is that the stress-strain curve is independent of the slip distance L. This is 
a most important result, because the discussion that we shall give below 
suggests that the slip distance may depend rather sensitively on the 
physical conditions. The analysis above shows how the stress-strain 
curve may not be so dependent. 

A second point is that, if we assume that / remains constant throughout. 
the hardening process, formula (12) gives a parabolic stress-strain curve, 
and if ~10-4 cm, it is of the form 


o/G~100 €1? ; 


the constant has about the value required for agreement with experiment. 
In this case, too, it is instructive to write the formula 


gmoytomatd [Fer <> Cal Sr Genel ae 


where a, is the yield stress. It is, however, likely that / will decrease 
during work-hardening, owing to the processes already described ; in that 
case the stress will increase as «”, where m is slightly greater than }. 

One criticism may perhaps be made of this theory. We have to ask 
why, when the stress is raised, further slip does not occur on planes that 
have already been active. All the dislocations in the active slip plane, 
when o>Gb/l, will have gone to the end ; there are none left in the main 
part of the slip plane. We suggest that further slip may happen on this 
plane ; but as the random internal stress is constantly changing, there is 
no reason why the particular centre that has slipped should be more 
favoured than any other ; and there are always far more sources than are 
used. In the same way, if the stress is reversed, it will normally be 
other sources that are active, and not those that have already slipped. 

We now turn to a discussion of the nature of the barriers. We remark 
first that various investigators (Rohm and Kochendérfer 1950, Andrade 
and Henderson 1951, Liicke and Lange 1952) have observed easy glide in 
crystals of close-packed metals. The conditions seem to be 

(a) Very pure materials. 
(6) Orientation for slip on one plane only. 
(c) In some cases absence of oxide layer. 

By easy glide we mean either the absence of hardening, or slight 
hardening compared with that normally observed, of magnitude comparable _ 
with that in hexagonal metal crystals. Easy glide is not accompanied 
by the formation of asterisms. It appears that in crystals showing easy 
glide—as perhaps in hexagonal crystals, the dislocations pass out of the 
crystals. We believe that when rapid hardening takes place the barriers 
which stop the dislocations, apart from grain boundaries and oxide layers, 
are likely to be sessile dislocations, either those already existing or those 
formed by slip on two planes. If this is so, easy glide on one slip system 


will occur most easily in very carefully annealed crystals of small total 
diameter—as in Andrade’s work. 
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This assumption about the nature of the barriers enables us to estimate 
the slip distance in the absence of other boundaries. If a reasonable 
proportion of the dislocations in a material are sessile, there should be 
b/P lines, each of length J, per unit area in any plane. If any one of these 
is capable of stopping a dislocation in the plane, and if the slip distance L is 
the same for the screw as for the edge component, we have 


bL?/B~1, 
whence 
L~ly/(1/b)~10-? em. 


This gives a value of L of the right order of magnitude, perhaps greater 
than 10-? cm, since certainly the distance between sessile dislocations will 
be greater than between the normal ones. Moreover, it gives for n 


n= Ao L/bG=L|l= /(Ib), 


a value which varies with J only as the square root, and may perhaps 
explain the comparatively constant value of n. On the other hand, 
conditions which promote flow on two or more sets of planes are bound to 


Fig. 6 


\ 
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A DEFORMATION BAND 


increase the number of sessile dislocations, and so decrease the slip distance 
and hence n ; we should certainly expect some decrease in 7 as the strain 
is increased. Further experimental work on this point would be of 
interest. 

§5, DEFORMATION BANDS 


Up to the present we have assumed that all slip bands, except those 
that reach the surface, terminate in isolated groups of piled-up dislocations. 
In fact, however, there are frequently observed slip lines which terminate 
in deformation bands ; a schematic representation of a deformation band 
is shown in fig. 6. Deformation bands are an origin of asterisms, and it 
has frequently been suggested that they are the origin of hardening.* 


* e.g. Mott (1951). Asterisms are not caused by isolated groups of piled-up 
dislocations. It should be possible in principle to harden a material without 
forming asterisms. 


4K2 
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In our view, however, a group of piled-up dislocations in a deformation 
band cannot contribute to the hardening at all, because a grid of disloca- 
tions or of groups of dislocations of this type has no long-range stress 
field.* Groups of dislocations which terminate in a deformation band are, 
from the point of view of work-hardening, wasted. 

A deformation band grows in our view in the following manner. The 
regions P and Q form, owing to the stress field there, very effective traps 
for dislocations of sign such as to lengthen the band. In a material in 
which the density of barriers is otherwise low, deformation bands will 
tend to grow. But unless the density of barriers is very low indeed, there 
is no reason why a high proportion of slip lines should terminate there. 

It is known (Chen and Mathewson 1951, Honeycombe 1951, Laloeuf and 
Crussard 1951) that deformation bands do not appear when a crystal is 
stressed under conditions such that slip occurs simultaneously on two sets 
of planes. On our model, this will be because, with slip on two sets of 
planes, a large number of sessile dislocations will be formed, leading to a 
low value of Z and inhibiting the growth of deformation bands. 

Cottrell (1952) has mentioned another factor which may contribute to 
strain-hardening when slip on two or more planes occurs. Under such 
conditions, a dislocation cannot move far without cutting other dislocations 
and so forming jogs. Energy is required for this, and so an increasingly 
large stress is required to drive a dislocation through the material. We 
believe that this cannot be the main cause of hardening ; if it were, then, 
for impure materials for which there is no question of easy glide, the 
hardening would depend markedly on the orientation. This is not the 
case. As Cottrell points out, however, the activation energy required to 
form a jog is of the order of a only few electron volts even in the absence of 
astress. Thus after plastic flow such dislocations as are held up by other 
dislocations will be easily freed by heat. Cottrell (1953) has given an 
analysis of logarithmic creep in pure aluminium which suggests that it is 
due to this process. We believe that Cottrell’s mechanism introduces a 
small temperature-dependent frictional term in the resistance to the motion 
of dislocations. 


§6. MovemEent or Screw DIsLocATIONS AND THE GENERATION OF 
VACANCIES 


The theory of hardening given in the last sections is based on the 
hypothesis that groups of piled-up edge dislocations can be formed through 
the agency of the sessile dislocations which results when edge dislocations 
on intersecting planes coalesce. We have made no mention of the screws, 
except to assume that their slip distance L’ is at least as great as L. We 
shall in this section discuss the behaviour of the screws ; they are 
important, because 
UL nn ie 

* ef. Nabarro 1952. 
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(a) We believe that moving screw dislocations generate vacancies, 
and that these are important for the explanation of the dependence of 
work-hardening on temperature, clustering of slip lines, electric resistance 
due to cold-work and other effects. 

(6) In hexagonal crystals or others in which only one slip plane is opera- 
tive, the blocking of edge dislocations through the formation of sessiles is 
impossible. We believe that in hexagonal crystals the slip distance L 
for the edge components must be the same as linear dimensions of the 
erystal ; apart from oxide layers there is nothing to prevent their passing 
right out of the crystal. The (much slower) hardening which occurs in 
hexagonal crystals may be associated with piled up groups of screws. 

(c) A discussion of the behaviour of screws is necessary in order to under- 
stand cross-slip. 

All these phenomena depend in our opinion on the properties of vacant 
lattice sites or vacancies, and we shall give a few preliminary facts about 
these. According to the vacancy theory of diffusion (Johnson 1939), 
a vacancy is characterized by the energy W required to form it and the 
activation energy U tomoveit. In terms of these quantities the coefficient. 
of self-diffusion D is given, apart from a numerical factor, by 


D=va? exp {—(W+ U)/[kT}. 


According to the calculations of Huntington and Seitz (1942), W and U 
should be about equal for a close-packed metal. Thus, if vacancies are 
formed artificially, e.g. by quenching from a high temperature or by 
neutron bombardment, self-diffusion is possible at about half the temper- 
ature at which it normally occurs. Recent work by Nowick (1951) has. 
confirmed that the exchange of atoms in an alloy (AgZn) is more rapid 
after quenching, the effect being explained along these lines ; Nowick 
finds for this alloy that, if this hypothesis is correct, U must be somewhat 
greater than W. ~ 

At temperatures at which self-diffusion is possible, the absorption or 
generation of vacancies enables an edge dislocation to ‘climb’, i.e. to 
move in a direction perpendicular to its slip direction. As described by 
the present author (Mott 1951), such a process is almost certainly respons- 
ible for the polygonization of a deformation band. 

There is much evidence that moving dislocations generate vacancies at 
temperatures too low for them to be formed thermally. The evidence has 
been fully reviewed by Seitz (1952). The two main facts quoted by Seitz 
are: 

(a) Cold-work increases the conductivity of alkali-halides, a conductivity 
thought to be ionic and to be due to vacancies. After cold-work the 
resistance gradually increases at a rate consistant with this hypothesis. 

(b) The increased electrical resistance of metals due to cold work can be 
partly or mainly removed by heating to a temperature much below that at 
which the mechanical ‘ recovery ’, or drop in the flow stress, occurs. It is 
thought that a large part of the electrical resistance is due to vacancies. 
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formed by the cold-work; this will disappear if the metal is raised to a 
temperature at which they can move; while to get rid of mechanical 
hardness one has to allow edge dislocations to ‘ climb ’ out of the piled-up 
groups (Mott 1951), which will occur only if the specimen is heated to 
near to the temperature of self-diffusion. 

Recent work at Harwell (Dugdale 1952) has given striking confirmation 
of the correctness of this hypothesis about the origin of the residual 
resistance due to cold-work. Platinum wires have been exposed to fast 
neutron bombardment in the pile, and an increase in the electrical resistance 
of the order of 4% observed. It is found that this resistance begins to 
diminish on annealing at about 100° c, and the activation energy for this 
process was determined, and found to be about 1-1 ev. Specimens with a 
residual resistance due to work-hardening were annealed in the same way 
and for these specimens too the activation energy for resistance decrease 
was found to be 1-lev. Since in the first case the resistance was almost 
certainly due to vacancies and interstitial atoms and the annealing 
involved the movement of the former, it seems highly probable that this is 
true also for the cold-worked specimens. 

We turn now to the mechanism by which vacancies are formed by 
moving dislocations. Seitz (1952) gives a review of the various possibili- 
ties; of these we consider the following the most likely. 

Consider the movement of dislocations in a given slip plane. Suppose 
that a certain number of screw dislocations cut this plane ; in this context 
by a screw dislocation we mean one with Burgers vector not in the slip 
plane considered.* We shall call such dislocations ‘ crossing dislocations ” 
and the points where they cut the slip-plane under consideration * crossing 
points’. Now, whenever a moving dislocation passes a ‘ crossing point ’, 
a jog is generated in that dislocation ; and although a jog once formed is 
no impediment to the motion of an edge dislocation, a screw dislocation 
with a jog in it cannot move without leaving behind it a row of vacant 
lattice sites or interstitial atoms, according to the sign. This geometrical 
result of the dislocation concept was first put forward by Thornton Read, 
and is fully discussed by Seitz (1952) and by Mott (1951). 

We may thus accept the hypothesis that moving screw dislocations do 
produce rows of vacancies and interstitial atoms, and that this process 
means that there is a resistance to the motion of screws which is not 
shared by edges. However, it is not immediately obvious that the passage 
of a group of dislocations along the slip plane, as in the formation of a 
slip band, will produce the vacancies and interstitials that we require. 
For, as pointed out by Mott (1951), provided the crossing screws do not 
move during the formation of a slip band, the rows of atoms added or 

‘subtracted by each succeeding dislocation will be side by side ; all that 


will result from the motion of n screws a distance L’ beyond a crossing 
ee ee ee eee 
In the simple cubic structure such a dislocation must have a screw component, 


at ae is also the case for hexagonal structures, but not for the close-packed 
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point will be the removal or addition of a slab of atoms of width nb and 
length L’. After the passage of the first few dislocations, the passage of 
further ones will meet with little resistance. 

Seitz has suggested that, by a mechanism to be described below, enough 
energy is released in the formation of vacancies to disperse them. _ If this is 
true for the whole slip distance L’ (say 10-* cm), far too many vacancies 
are formed and the resistance to the motion is much too large to be 
compatible with the observed facts. In order to explain the observed 
resistance, Seitz has to assume that a moving dislocation forms one vacancy 
for every 10* atoms that it crosses. If 1/I?, the density of crossing points, 
is such that /~10~4 em, this will be the case after a dislocation has moved 
a distance of the order / ; after it has moved a distance L’ it will have jogs 
at intervals of (*/L’ along its length, and, if L’ is 10-? cm, will be creating 
vacancies at 100 times the observed rate. 

Moreover, if x is the distance between jogs, the resistance to the motion 
of a screw is equivalent to a stress of order Gb/x, and only if x is greater than 
Zis this comparable with or less than the yield stress. 


Fig. 7 


(a) (b) (c). 


Motion of a screw dislocation cutting another screw. 


We believe that the correct model is a compromise between the two 
points of view. Each screw dislocation, when it approaches a crossing 
point C, is repelled by the crossing dislocation, and curls round it, as in 
fig. 7(a). When the two edge components, now on contiguous planes, 
join together and form a row of vacancies, enough energy is released to 
_ disperse them. Thereafter, however, and thus at a distance of order / 

from ©, a continuous screw line is formed as in fig. 7 (6), which leaves 
behind it a row of vacancies which are not dispersed. The first process 
generated about J/b dispersed vacancies every time a screw passes a crossing 
point ; the second generates rows of vacancies next to each other, and thus 
removes layers of atoms, presenting no effective resistance to the movement 
of dislocations after the first few. 

Such a mechanism accounts for the generation of about the right number 
of vacancies ; we have to ask further 

(wz) Can we give any explanation of the limited slip along one plane, 
similar to that provided by sessile dislocations in the case of edges ? 

(b) Can we find a mechanism which prevents the screws from moving 


back under a reversed stress ? 
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(c) Can we account for the (slow) work-hardening of hexagonal single 
crystals ? . 

As regards (a) we suggest the following: if screw dislocations have 
moved a distance L, for each of the L/?? crossing dislocations cut per unit 
length, a plane of atoms of width nb is either added or removed. The 
crossing dislocations are pushed sideways, therefore, a distance nb, for 
which energy is required to increase their length. We may suppose that 
the resultant force pushing the n screws back is equivalent to a stress 


GnLb?|I°. 


Which increases with increasing 7. 

The model suggests, then, that the slip on each band increases with 
increasing strain. That the slip on each plane increases with increasing 
strain is in accordance with the observations of Andrade and Roscoe 
(1937) on zine. 

The model will also lead us to expect that as the strain is increased the 
steps on the surface of a cubic metal may lengthen, though they do not. 
increase in height (at low temperatures). 

With regard to the sticking or anchoring of screw dislocations, we suppose 
that they will always have a certain number of ‘ free jogs ’~i.e. jogs 
which are not anchored to a point such as X in fig. 7 (c). These, and also 
any movement of the crossing screws since the slip bands were formed, will 
hinder movement in either direction. 

Finally, we discuss the work-hardening of hexagonal crystals (or of 
cubic crystals under conditions of easy glide). This may by essentially 
of the same nature as in cubic crystals, the screws replacing the edges. 
We suppose that in hexagonal crystals the edges escape to the surface, so- 


that only the screws are stuck, the worked crystal being full of screws. — 


The slow rate of hardening of hexagonal crystals will be due primarily to- 


the much greater tendency that one would expect, there being no possib- | 


ility of sessile screws, for the groups of dislocations to form into grids, 
with the result that there is no long-range stress from any group. Some 
evidence that this occurs is presented in the next section (on cross-slip). 
We have not, however, been able to formulate any quantitative theory for 
these crystals. 

An alternative suggestion is that in hexagonal crystals and in cubic 
crystals under conditions of easy glide the screws escape as well as the edges, 
and that the hardening is simply due to the strains round vacancies formed. 


§7. CROSS-SLIP 


Up to this point we have implicitly assumed that dislocation loops remain 
in the plane in which they start to move. If this were the case, slip bands. 
observed on the surface of metals would be straight. Actually, in cubic 
metals, the slip bands observed frequently show intimate cross-slip. 


Examples of cross slip have been given by Cahn (1951). The points to 
note are: 


<-> 
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(i) Slip bands on the (111) planes on which the stress is largest are 
connected by short bands on the other sets of (111) planes. 
(ii) Cross slip is more marked at high temperature. 

(iii) The photographs at high magnification show that the step connec- 
ting the main bands is less sharp than the main steps. ; 
An explanation of cross-slip (due to Frank) was given in the author’s 
Guthrie lecture (Mott 1951). This explanation depends on the different 
kinematical properties of the edge and of the screw dislocation. The 
edge dislocation, if it is straight, can move only in the plane defined by its 
Burgers vector and itself. Moreover each element of an edge dislocation 
line which switches from one (111) plane to another has this property too. 
The same is true of dislocation lines which are partially edge, partially 
screw. A pure screw dislocation, on the other hand, has its Burgers 
vector parallel to itself, and so is not associated with any particular one 


<—— 
Stress 


(c) 
Model for cross-slip. — 


of the (111) planes in which the Burgers vector lies. There is no kinematic 
reason, therefore, why a screw dislocation should not move for a certain 
distance in one of the (111) planes and then switch to another. 

The explanation of cross-slip already referred to envisages, then, two 
loops spreading out from sources, marked §, 8’ in fig. 8, in parallel planes 
not too far from each other. PQ, P’Q’ are the screw elements. When 
PQ reaches a position above P’Q’ the two lines having opposite sign attract 
each other. They can then leave the planes in which they were moving 
originally and move towards each other until they join and annihilate 
each other. One is then left with a ring in which the two edge elements 
switch from one plane to another (fig. 8 (b)). This may happen for a 
number of sources, so that finally one is left with the two kinky edge 


1172 N. F. Mott on a 


dislocations of fig. 8 (c), connected by screws. When these edge dis- 
locations appear on the top surface of the crystal, and are followed by 
others produced in the same way, they will form kinky steps of the type 
observed. 

This model shows why the steps joining the main slip bands are less 
distinct ; when one loop after another starts out from the sources 8, 9’, 
there is no obvious reason why the screws should move from PQ to P’Q’ 
always along the same plane. 

We believe that this description is geometrically correct. It needs, 
however, some amplification. One may ask, for instance : 

(a) Why should the sources SS’ be generating loops at the same time ? 
They will not in general have the same values. of | and o; and one 
would expect one source to generate its full complement of loops before 
the next was brought into play (cf. Crussard 1951). 

(b) In close-packed cubics—as opposed to the simple cubic, a very 
considerable activation energy will be required to make a screw 
dislocation moving in one plane switch to another. This is because, as 
already stated, all dislocations, screw or edge, split into half dislocations. 


Fig. 9 


Showing activation process for the motion of an extended screw dislocation 
out of its plane, the dotted lines being in new plane. Each line represents 
a half dislocation. 


The half screw dislocation does not have its Burgers vector parallel to 
itself; it cannot move out of its plane. Ifa screw is to move out of its 
plane, a small length must be moved out first, and at the points where 
the two join, the two halves must be brought together, as in fig. 9. The 
necessary activation energy will be several electron volts, probably too 
large for the process to occur. 

We believe that the clue to both criticisms is to be found in another 
property of the model which has already been emphasized. A screw 
dislocation, after it has cut a certain number of other screws, contains 
jogs, so that the resistance to the motion of a screw will increase as it 
moves along, and other screws will pile up behind it. The stress induced 
by this pile-up will readily induce some nearby source to generate loops, 
so that one is left with the configuration shown in fig. 10. Slip has 
occurred from A to B and A’ to B’; energy will be gained if the piled-up 
screws at B, B’, can move into the plane BB’ (shown dotted) and so 
annihilate each other. The experimental evidence shows that this only 
occurs when the temperature is raised ; we suggest that the clue is the 
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high stresses due to the pile-up ; these either allow spontaneous generation 
of loops in the plane BB’, or much reduce the activation energy for the 
process of fig. 9. 

In a recent paper Smith (1952) has shown that pure copper bars 
extended by c. 10% show a drop in Young’s modulus of about 10%, but 
that this is recovered on annealing at 100°c. This may well be due to the 
motion of screw dislocations bounding the slip bands (the edges being 
locked), and the process by which annealing takes place exactly the same 
as cross slip, namely the mutual annihilation of the screws. 


§8. DEPENDENCE OF WORK-HARDENING ON TEMPERATURE 


The theory of work-hardening of cubic crystals given in §§ 2, 3 and 4 
is applicable as it stands to work-hardening at low temperatures only. 
In this section we shall suggest explanations of the following phenomena 
which depend on temperature : 

(a) The ‘ clustering’ of active slip planes at distances of the order 
of 200 A from each other observed in aluminium at room temperature and 
above and in other metals, but absent according to Brown (1951) in 
aluminium at the temperature of liquid air. 


Fig. 10 
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Showing groupings of screw dislocations. 


(b) The dependence of the whole stress-strain curve on temperature. 
(c) The behaviour observed when a crystal is strained at room 
temperature and then strained further at a lower temperature. 

We remark first that the theoretical values of all the flow stresses 
derived in this article are linear functions of G, the shear modulus, and 
that G@ itself depends on temperature ; for aluminium the change is of 
the order 7% per 100° c (Ke 1947) above the characteristic temperature. 
Thus changes of from 10 to 20% between the absolute zero and room 
temperature may be expected from this cause alone. Apart from this, 
one expects very little temperature-dependence for a yield stress 
determined by a relation of the type 

o,= GOI. 
The activation energy necessary to operate a Frank—Read source for 
stress o less than o, is (Mott and Nabarro 1948) 


ha, bE (1 olen eens or 03 ree (18) 
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where k is a numerical factor of order 0:2. The term o,6/? is equal to 
Gb3(i/b) and thus to c. 104 ev. If, then, a yield point with a dependence 
on temperature greater than G is observed, as for instance in. Rosi and 
Mathewson’s (1950) work on aluminium, the most probable conclusion 
is that the yield point is determined by Cottrell locking.* 

Another process in which a large dependence on temperature may be 
expected is the spontaneous formation of dislocation loops in regions of 
very high stress; this has already been mentioned in connection with 
cross-slip. Another is Cottrell’s mechanism of the cutting of one 
dislocation by another ; this we believe responsible for transient creep 
at low temperatures. 

Apart from these and the variation of G with temperature, it seems 
impossible to account for any temperature dependence of work-hardening 
in terms of an activated motion of dislocations ; the term / in formula (15) 
should be at least 10-4em or possibly even comparable with the slip 
distance. We consider, therefore, that the dependence on temperature 
is largely to be accounted for through the action of vacancies. 

We assume that thermal recovery is associated with the ‘ climb’ of 
dislocations out of their slip planes, a process possible at temperatures 
for which self-diffusion can occur ; a theoretical account has been given 
by the author (Mott 1951). Such a process enables the piled-up arrays 
of dislocations, to which we have ascribed work-hardening, to disperse 
and deformation bands to polygonize. We suggest, however, that at 
considerably lower temperatures, namely those at which vacancies can 
move, the vacancies created by the moving screw dislocations will soften 
the material in the immediate neighbourhood of the planes which have 
already slipped. This could occur through the diffusion of vacancies 
to the secondary groups of piled-up dislocations, the existence of which 
has already been deduced. The temperature above which this can occur 
will be c. —70° c in aluminium, room temperature in copper, c. 120° ¢ 
in nickel. f 

If this occurs we have, following Brown, an immediate explanation 
of clustering. Once slip has occurred on a given plane there will occur a 
relatively slow reduction of co, in the immediate neighbourhood of that 
plane. Slip will then occur again, either on the plane which has slipped 
already or on some neighbouring one ; but the source that has generated 
vacancies already is in no way preferred ; in most cases it will be another 
one. The nearest of these should be at a distance from the first plane 
BILL’; with L~L’~10- cm and 1/~10! em-%, this gives 10-6 cm. 
ae URE 

* Note that Neurath and Koehler (1951) find almost the same yield point of 
pure single crystals of copper at liquid air and room temperature. 

+ Brown (1951) has suggested, without proposing a mechanism, that the 


material round a slip plane anneals, and has explained clustering and certain 


other phenomena in terms of this assumption. The discussion given here 
owes much to Brown’s work. 


{ The basis for this estimate is given below. 


Theory of Work-hardening of Metal Crystals 1175 


We thus expect sources at about this distance to become active. The 
lines of a cluster, on this hypothesis, will be similar in nature to those of the 
primary slip, and will, like the primary slip band, have groups of piled-up 
edge dislocations at the extremities. We suppose, moreover, that the 
limited supply of vacancies affects the secondary rather than the primary 
groups of piled-up dislocations. 

If this hypothesis is correct, we may, following Brown, account for the 
temperature-dependence of the work-hardening curve. At temperatures 
T), T, (fig. 11), too low for vacancies to move, we shall expect a small 
change in o proportional to that for G. At a temperature 7’,, at which 
vacancies can move, we shall expect : 

(a) A rather rapid drop of the whole curve, it being assumed that most 
new slip bands start near old ones. 

(6) The disappearance of the main part of the residual resistivity 
produced by cold-work. 

(c) The formation of clustered slip lines. 


Fig. 11 
To 
Th 
o 
T2 
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Work-hardening at temperatures such that 7,>7,>7,>T7. 


At a temperature 7’, more than twice as high, true recovery may occur 
during straining. 

Finally, following a suggestion of Brown, we can explain with this model 
the behaviour of the stress-strain curve of a crystal such as aluminium 
strained at room temperature and then lowered to liquid air temperature 
(fig. 12). On cooling to liquid air temperature there is a small jump in o, 
doubtless due to the change in G ; then hardening is very rapid until the 
curve characteristic of liquid air temperature is approached. Now in our 
model the quantity n/LZ—i.e. the strength of the groups of piled-up 
dislocations divided by the distance between them, will not depend on the 
temperature of deformation; moreover it is this quantity, according to 
formula (11), which determines the flow stress. The crystal at room 
temperature flows more easily because the material round any slip plane 
softens, owing to a rearrangement of the secondary groups. Thus when slip 
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continues from the point A in fig. 12, the annealed region is quickly 
used up ; thereafter the material differs from a crystal strained from the 
beginning at liquid air temperature only in that some of the slip bands are 
grouped in clusters, but not in the total number of dislocation groups or in 
the mean value of n/L. 

A further consequence of the formation of vacancies is probably, as has 
been suggested by various authors (Hall and Williamson 1951), the 
fragmentation or substructure that is formed in a cold-worked material 
(Gough and Wood 1936, 1938). The elements of the substructure are, 
we suggest, simply polygonized crystallites, and polygonization is enabled. 
to take place through the movement of vacancies formed by the cold work. 
This suggests that fragmentation would not take place at a low enough 
temperature. 


Fig. 12 


180°C 


Stress strain curve initially at room temperature and 
then at liquid air temperature. 


Finally, we give a very rough estimate of the critical temperatures 7’, 
below which the vacancies cannot move. Unfortunately we have only 
one definite measurement of the activation energy U for movement of a 
vacancy, that of Dugdale (1952) for platinum, and this is less than one half 
the probable value of W for this element, instead of being equal to it. 
The following table gives U-+- W, the activation energy for self-diffusion for 
the elements for which it is known, and for Al, Pt, Ni deduces it on the 
assumption that these quantities are proportional to the melting points. 
U is deduced on the assumption that the ratio U to W is the same as for 
platinum. Finally 7', is taken, rather arbitrarily, as the temperature at 
which a vacancy makes one jump in ten seconds. 

Commenting on these results, we may note that Seitz (1952, p. 45) reports 
unpublished experiments by Blewitt which record that the resistance of 
copper strained at liquid air temperature shows a partial recovery on 
annealing at —80° © ; with the above value of U, we may estimate that 
at this temperature a vacancy would change places about once in 1 000 sec, 
soananneal at this temperature is just possible. It is difficult on this basis 
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to understand the results of Bowen, Eggleston and Kropschott (1952) who 
after cold work at room temperature found an activation energy of 
28-3 kcal for further recovery of the residual resistance of copper. Perhaps 
this is the activation energy for the motion of interstitial atoms, 


Table 1 
Metal W+U U T', (degrees K) 
Cu 46-5 16 230 
Ag 45-9 16 230 
Au 53:0 18 260 
Pb 27:5 9-5 140 
y-Fe 67:9 23 350 
Al 39 13 200 
Pt 15 26 430 


Ni 64 23 390 


W and U are in keal/g atom. 
Experimentally determined quantities in bold-face type. 
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SUMMARY 


In a previous paper (Budden 1951 a, referred to as Paper I), it was 
shown that the propagation of an atmospheric could be treated either in 
terms of the successive reflection of rays at the earth and the ionosphere, 
or by considering the space between the earth and the ionosphere as a 
wave-guide, and discussing the properties of the various wave-guide 
modes. An outline of the theory was given for the case when the top 
surface of the wave-guide was the boundary of a homogeneous ionized 
medium, and the earth’s magnetic field was neglected. 

The present paper gives the full mathematical theory, and shows how 
the characteristics of the wave-guide modes may be determined in the 
most general case, provided that the reflecting properties of the earth 
and the ionosphere are known as functions cf the angle of incidence 
which is, in general, complex. It is found that the wave-guide modes 
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are of two types, which are described as ‘ quasi-transverse magnetic ” 
and ‘ quasi-transverse electric’. Allowance can be made for the curva- 
ture of the earth. 

A particular case is then discussed, in which the surface of the earth 
is assumed to be perfectly conducting, and the ionosphere is assumed to 
be a homogeneous ionized medium, with the steady magnetic field of 
the earth superimposed. The results of numerical calculations are 
given, for a few special cases, in the form of curves. These show (i) the 
attenuation in the various modes as a function of frequency, (ii) the 
polarization characteristics of the wave in typical modes, and (iii) the 
amplitudes of the waves which are excited in typical modes by a vertical 
electric dipole source. 


$1. INTRODUCTION 


In a previous paper (Budden 1951 a, hereinafter referred to as Paper I), 
it was shown that the propagation of an atmospheric could be treated 
either in terms of the successive reflection of rays at the earth and the 
ionosphere, or by considering the space between the earth and the ionosphere 
as a wave-guide, and discussing the properties of the various wave-guide 
modes. The method of wave-guide modes was shown to be the simpler 
for propagation to great distances, since there are then only one or two 
modes of large amplitude, and the modes of higher order are heavily 
attenuated. The characteristics of the various modes were evaluated 
by a method which was essentially equivalent to Eckersley’s phase 
integral method (Eckersley 1932), and a physical argument was advanced 
to justify the use of this method. A mathematical justification of the 
method was given only in outline in Appendix A of Paper I. The effects 
of the earth’s magnetic field and of the earth’s curvature were neglected. 

The object of the present paper is to give a fuller mathematical treat- 
ment of the same problem, and to show how the theory can be extended 
to include the effects of both the earth’s curvature and the earth’s magnetic 
field. As in Paper I, the ionosphere is still treated as homogeneous and 
sharply bounded, but is no longer regarded as isotropic. 


§ 2. SomrE GENERAL THEORETICAL IDEAS 
2.1. The Reflection of Very Long Radio Waves at the Ionosphere 


Recent experiments on the propagation of very long radio waves from 
commercial radio transmitters have shown that reflection of these waves 
at the ionosphere is accompanied by a change of polarization. (See for 
example, Bracewell ef al. 1951.) It is convenient to describe the incident 
and reflected waves in terms of component electric fields parallel to, and 


perpendicular to, the plane of propagation,* and to consider four co- 
Ee ee ee eee 


* The component in the plane of propagation is sometimes called the ‘ normal 


component ’, and that perpendicular to the plane of propagation the ‘ abnormal 
component ’. 


: 
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efficients denoted by it, ,2,, ,R,, ,R,, to indicate the complex ratio 
of a specified electric field in the wave after reflection to a specified 
electric field in the wave before reflection. The first subscript denotes 
whether the electric field specified in the incident wave is parallel (||) or 
perpendicular (L) to the plane of incidence, and the second subscript 
refers in the same way to the electric field in the reflected wave. 

The set of four conversion coefficients can be written in matrix form 


thus: 
R R 
R= |e oh eo a i er ieaiter 
we,  R ) 


1 


In this paper, as in Paper I, the ground will be treated as a perfect 
conductor. For reflection at the ground, the reflection coefficient matrix 


therefore becomes 
it 0 z 
R= ie os RRA eS) 


_ If we consider two successive reflections, the first at the ground, and the 
second at the ionosphere, the pair of reflections is together represented 


by the matrix 
R,, — Rk 
R= |" a eet eh. eee eo) 
; h — ft, 


[sgl 


It is readily shown that any number of successive reflections is represented 
by a matrix which is the product of the matrices appropriate to individual 
reflections, multiplied in the correct order by the ordinary rule for 
multiplication of matrices. 


2.2. Sources and Image Sources 


In Paper I the source of an atmospheric was assumed to be a vertical 
electric dipole, whose dipole moment varied with time in a known way. 
The dipole moment was Fourier analysed into sinusoidally varying 
components, and the propagation of the component frequencies was 
considered separately. In finding the impulse response of the wave- 
guide system, a Fourier integral was evaluated which included a factor 
representing the propagation characteristic of the system as a function 
of frequency. In the present paper we confine our attention to one 
component frequency of the source. The source is therefore considered 
to be an oscillating vertical electric dipole of angular frequency p. If the 
impulse response of a system is required, it can be found by evaluating 
the Fourier integral which contains the transmission characteristic as 
a function of p. 

The signal at a point distant from the source can be considered as 
made up of rays travelling direct from the source, together with rays 
reflected one or more times at the ionosphere and at the ground. It is 
well known that these rays could be produced by an array of point 
sources placed at the positions of the images of the source as seen by 


4L2 
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reflection in the two surfaces. If the source is at the ground, and the 
ionosphere is at height h, this array of images acts like a diffraction 
grating to give spectra at angle 7 with the vertical given by: 


2h cosi=NA, ©. 5 eee 


where N is an integer, and A is the wave-length. 

The zero order spectrum forms the zero order mode in the wave-guide. 
The two spectra for N=+1, and N=—1, may be regarded as the 
component waves of the mode of order 1. If now the ionosphere is 
taken to be an imperfect reflector, the image dipoles are no longer equal, 
but become weaker as we go further away from the central one. The 
images correspond to two diffraction gratings, one above and one below 
the ground, in each of which the sources are not all of equal intensity. 
A series of spectra is formed as before, but in each spectrum the energy 
is now spread over a small range of angles.* 

The following way of considering the spectrum of one order, produced 
by the grating of unequal sources, is also useful. It is possible to find 
a complex angle, iy, for each spectrum such that the energy in the 
spectrum may be considered to be confined entirely to a wave whose 
normal is inclined at this complex angle to the vertical. Consider the 
contributions to the Nth order spectrum of two adjacent points in one of 
the gratings. These contributions consist of two plane waves whose 
(complex) phase difference is 2h costiy. The ratio of the two contri- 
butions is given by: 

R, exp (—2jkh cos ty). errs ee ly 


R, is a complex number representing the effect of two reflections, one — 


at the ground, and one at the ionosphere. If the ionosphere is isotropic, 
&, is a pure number, but we shall later discuss the case where R, is a 
matrix. By analogy with an ordinary grating, having sources of equal 
amplitude, we may infer that the components from adjacent points of a 
grating must agree in both amplitude and phase, and the expression (5) 
must therefore be unity. We then have: 


exp (27kh cos ty)=Raexp(27N). . . . . . (6) 
‘This is the same as eqn. (3) of Paper I. 

Now consider the field which results from combining the spectrum of 
order +N from the top grating and that of order —N from the bottom 
grating. From symmetry it is clear that the resolved parts of the 
propagation constants in the horizontal direction for these two spectra 
are the same. In fact these spectra are the upgoing and downgoing 
component waves of the wave guide mode of order +n. It is clear on 
this picture why the properties of modes of order +-n and —n are different. 


The situation is similar to that which occurs with an optical Lummer-— 


Ghercke plate. The finite angular spread in each spectrum corresponds to 
the finite resolving power of the plate. 


: 
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Modes of one sign relate to spectra of the gratings for which the wave 
normal is inclined towards the direction of increasing amplitude of the 
grating points, and those of the other sign to spectra with the wave 
normal inclined towards the direction of decreasing amplitude. It may 
at first be difficult to see how an array of images above the top boundary 
of the guide can give rise to an upgoing wave (real part of cos 7 y positive) 
below the boundary. Such an upgoing wave must, however, exist, to 
satisfy the boundary conditions. A fuller mathematical treatment of 
the problem is given in § § 3.1 to 3.3. 


2.3. The Hertz Vectors 


The vertical dipole source will be assumed to have a dipole moment 
P,, exp (jpt). The field of such a dipole is conveniently represented by 
the electric Hertz vector II,. This is everywhere parallel to the dipole 
source, and its magnitude is a scalar whose value at any point in a medium 
of dielectric constant « and unit permeability is 


IT, =P, exp (jpt) exp (—jkr)/r, tn LIT see, 
where r is the distance from the point to the source, and 27/k is the 


wavelength in the medium. The electric and magnetic fields of the 
dipole are then given by 


E,=(1/e9e) grad div II,—p, eM, jo, . . . . . (8) 
H,,—9/ot (curl II), Sf So py eat ee oe ER ee eras () 


where «<, and py are the electric and magnetic permittivities of free space. 
If the waves reflected at the earth and the ionosphere have no abnormal 
component, then the image dipoles are all vertical electric dipoles, and 
the fields due to them can also be described by electric Hertz vectors. 
The full expression for J, is then a series of terms representing the 
contributions from the various image dipoles. Such a series was given 
in Appendix A of Paper I. 

If the first reflected wave from the ionosphere includes some abnormal 
component, we may consider this as coming from the same image point 
as that which gives rise to the normal component, but the image source 
must now include a part which gives rise to a horizontal electric field. 
This part must have the same polar diagram as the source dipole. An 
image with these properties is a vertical magnetic dipole. Hach image 
therefore consists of an electric and a magnetic dipole superimposed. 
The field of a vertical magnetic dipole P, exp(jpt) at the origin is 
conveniently described by the magnetic Hertz vector II, which is always. 
parallel to its source. Its magnitude is given by 


IT, =ceyP, exp (jpt) exp (—jkr)/r. 
The factor ce, is included so that J7, and JT, have the same dimensions. 
The associated electric and magnetic fields are given by 
E,=—+,c0/ot( curl II, ), PE he ne ee ne) oe sta LO) 
H,=c grad div II, —(«/c)0°II,/0#*, Pry Le ae | 
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where c is the velocity of light in free space. The complete field at any 
point is therefore described by the single column matrix 


ae Joy eee eile 


The field in each component wave of any mode can also be represented 
by a matrix of this form, and the field after reflection of such a component 
wave is given by the matrix product: RII. 

The contribution to II from a vertical electric dipole is 


P{ exp (jpi—jly) Hr, ) 


Peale oni “Sokvat te cae 


2.4. The Wave-Guide Modes 


The space between the earth and the ionosphere may be regarded as 
a wave-guide of infinite horizontal width, but of finite vertical depth h. 
If the ionosphere is isotropic, then in such a guide modes of two types 
may be propagated, usually described as TM and TE modes. In the TM 
modes the magnetic field is everywhere horizontal, and the electric vector 
is entirely in the plane of propagation. In the TE modes, the electric 
field is horizontal, and the magnetic field is entirely in the plane of 
propagation. Only TM modes can be excited by a vertical electric 
dipole and hence in Paper I only TM modes were discussed. Any TM 
mode has an order n, which may be a positive or negative integer or zero. 
The zero order TM mode played an important part in the theory of 
propagation of atmospherics. If the earth and the ionosphere are perfect 
reflectors, the zero order TM mode is unattenuated at all frequencies, 
and becomes what is sometimes called a TEM mode. 

Now suppose that one surface of the wave-guide is doubly refracting. 
Then there are no purely TM or TE modes. The modes are, however, of 
two types. Those of one type are very similar to TM modes, but the 
magnetic field has a small component in the direction of propagation, 
and modes of the other type are similar to TE modes, but the electric 
field has a small component in the direction of propagation. We shall 
call these ‘quasi-TM’ and ‘ quasi-TE’ modes, respectively. It is 
now possible for modes of both these types to be excited by a vertical 
electric dipole, although we shall show (§ 4.6) that the amplitudes of the 
quasi-TE modes are small. 

Since quasi-TE modes may play some part in the propagation when 
the ionosphere is not isotropic, it is useful to discuss the propagation of 
‘TE-modes in the simple case where the earth and ionosphere are perfect 
reflectors. In such a guide the TE-modes could be excited by a vertical 
magnetic dipole, and we can discuss the excitation of the modes, using a 
point diffraction grating formed by the source and its images. Consider, 


where 
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first, the image of the source formed by the ground. It is a vertical 
magnetic dipole of equal amplitude but opposite sign to the source. 
The source and its image may be considered to be a single element of the 
grating. A succession of images of this pair of dipoles is formed by the 
two surfaces of the wave-guide, and these form the remaining elements 
of the grating. Now a single element cannot radiate in a direction 
perpendicular to the grating, because the radiations from the two 
component dipoles in each element cancel. Hence the grating can have 
no spectrum with its wave normal parallel to the surfaces of the wave- 
guide. The spectra, or wave-guide modes, of lowest order occur where 
2h cosi=+X, and in general modes of higher order occur where 
2h cos i=-+-nd, where n is an integer other than zero. 

The values of cos for the various modes of both kinds are given by 
the positions of the poles of an integrand (eqn. (19), § 3.2, eqn. (30), § 3.3). 


§ 3. DERIVATION OF THE MATHEMATICAL EXPRESSIONS 


3.1. The Total Field of a Vertical Electric Dipole 


We now derive an expression for the field in the wave-guide due to a 
vertical electric dipole source. Cylindrical polar co-ordinates (p, 4, 2) 
are used, and the source is at (0,0, H). The component of the field 
contributed by direct radiation from the dipole is given by (13). The 
space part of this expression may be written : 


{exp (—jkr)}/r—=4 jk | Hi,(&Sp) exp {FjkO(2—H)}dO, . . . (18) 
—sign if 2>H, 
+sign if -<H, 


where S=sini, C=cosi. (See, for example, Ott (1942), Sommerfeld 
(1935).) The integrand has branch points at C=+1. The integral is a 
contour integral, and the contour must be indented above the real axis 
at C=-+1, and below it at C=—1. The integrand is the equivalent of 
a set of plane waves distributed symmetrically in azimuth, and with 
their wave normals all at an angle 7 to the vertical. The field due to the 
original dipole is therefore given by 


T=1jk | Hig(hSp) exp peel Hedc >. +. (18) 


where P is given by (14). 

As the evaluation of the total field for the series of image dipoles is 
complicated, we deal first with the case where the upper medium is 
isotropic, so that the reflection coefficient matrix R, reduces to a pure 
number R, (a function of C), and R, is unity. It may be helpful to 
visualize R, as a Fresnel reflection coefficient for a sharp boundary, 
although the analysis is not applied specifically to this case until § 4. 
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3.2. The Residue Series—Isotropic Medium 


There are four series of images at heights Imh+H, 2mh—H, —2mh+H, 
and —2mh—H, where m is a positive integer. The fields of these are 


given by 
Th ayk J Hij(kSp)Pf, m(C) dC, (r=1,2,3,4) . . (17) 
where 
fi, m=[exP {jkC(z—H—2mh)}|R.”, (m=OAl eae) (18a) 
fo, m=[exP {jkO(z+H—2mh)}|R.”, (m=1, 2,...) (185) 
fs, m=leXP {—jkC(z—H-+2mh) }|R,”, (m=1, 25.8 (18c¢) 
fa m=lexp {—jhC(e+H+2mh)}1R.". (m=0,1,2,...) (18d) 
The total field is given by 
n=1jk | Hi, (6p). Pf(C) dG. ee ee 
where 
AAO tN hd m+ (fe, moe as ~ : < “4 - (20) 


The series (20) are geometrical progressions, but before they can be summed 
we must examine their convergence. With the sign convention used in 
Sommerfeld’s expression (15), the reflection coefficient R, reduces to the 
familiar reflection coefficient when C' is real and positive. When C is 
negative, R, is in general greater than unity. For convergence of the 
series it is necessary that 


| Ry exp (—2ja) |<1) 2 


where « stands for hkC. This holds on the positive real axis, but not 
on the negative real axis. We may, however, deform the contour in 
the integral (19) so that (21) holds everywhere on it. The original 
contour is shown in fig. | a. Let QQ’ be the line whose equation is 

|.R, oxp (—27a) [== ee 
It can be shown that in the third quadrant the branch points of R, lie 
below the line QQ’, and there are no poles. The integrand vanishes at 
infinity, in the first and third quadrants, due to the term Hi,(kSp). 
Hence the contour may be deformed to the position A, fig. 1 6, without 
crossing any singularity of the integrand, and (21) will then be valid. This 
applies for all the series (18 a to d). When this deformation is made, we 
may sum the series as geometrical progressions, and we have from (20) 


J(C)=2 exp (—jkCH) . cos (kOz)[1—R, exp (—2ja )}-} 
+2 exp (jkCH) . cos (kCz)R, exp (—2ja) [1—R, exp (—2ja)]. (23) 
Suppose that in the integral (19) we replace C and 8 by their complex 


conjugates, c* and S*; this is not the same as replacing the whole 
integrand by its complex conjugate, since we do not replace 9g by —), 
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except where it occurs in the variables C and 8. The integrand is there- 
fore replaced by a new function, and in particular it has different 
singularities. The values of C and S on the real axis are, however, 
unaltered by the change, and provided that we now use the contour of 
fig. 1c, the value of the integral is the same as before. To make the 
series (18 a to d) convergent, we now deform the contour to the position 
B as shown in fig. 1 d, and the series can then be summed as geometrical 
progressions as before. 


Fig. | 


Contours for integration in the complex C-plane. 


We can therefore write (19) thus: 
=1jk { Hi,(bSp) . P.f(C) .dD+4jk | _Hig(hS*p) . P .f(C%) . d0* 
: t 
(24) 


where f(C) is given by (23). f(C) and f(C*) vanish everywhere at infinity. 
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Let the values of S, C, « at the poles of f(C) be S,, Cy, %,. Then f(C) 


has poles where 


2a,,+-§j,log, RalC,) = 27) (eee eae 


(n is an integer), and there is a similar equation for the poles of f(C*). 

Now suppose that the contours A and B are closed as in figs. 1 e and 1 be 
into contours [’, and I’y. When the radius of the curved parts of these 
contours becomes indefinitely great, the contributions to the integrals 
from these curved parts vanish because of the term Hi (kSp). The 
contributions to the two integrals from the imaginary axis are equal in 
magnitude and opposite in sign since there S and S* are purely real, 
and equal. Branch points of R, (and therefore of f(C)) are excluded 
from the contours by cuts as shown (figs. 1 e, andf). (There may be more 
than the two branch points shown.) The contributions to the integrals 
from these cuts are associated with waves of a type discussed by Sommer- 
feld (1909) (the surface wave), and by Ott (1942) (the ‘head’ wave). 
For the particular cases discussed in § 4, the branch points occur where 
the imaginary part of S is large. The associated waves would therefore 
be heavily attenuated, and would make a negligible contribution to the 
signal received at a distance. They are not wave-guide modes, and we 
shall not discuss them. We shall evaluate only the contribution to the 
integral from poles, and we use the symbol J7,, to denote this contribution. 
Previous authors (Bremmer 1949, Rydbeck 1944) have called this 
contribution the ‘ residue series ’. 

The poles of f(C) in I’, are on the line OQ’, and occur for positive 
(or zero) values of n in (25). Poles of f(C*) in I’, are on OP’ (figs. 1 d, 
1 f) and occur for negative values of n. The residue of f(C) at the pole 
els 


n 


4{cos (kz0',,) cos (KHC,)}/{jhk—Ro'2Ro}, . . . . (26) 


with a corresponding expression for poles of f(C*). A dash, ’, denotes 
differentiation with respect to C. Hence 


11,=(2n3|h)) =a Hig(kS,,,p) cos (k20',) cos (KHC,){1+jRy'/(2hkR,)}-} 


+ 2) Hild," p) cos (kzC,,*) cos (HC, *) {1-4 ]y [QMKR,)}—* | . 
Ss seg MO 


The term Ry'/(2hkR,) is small except at very low frequencies. It was 
neglected in Paper I, and the resulting equation was quoted in Paper I, 
eqn. (8). The terms of the series (27) refer to the various wave-guide 
modes. The characteristics of the modes are determined by eqn. (25), 
which gives the positions of the poles. This is equivalent to Kekersley’s 
phase integral relation, and was discussed in Paper I. 


Propagation of a Radio Atmospheric—II 1189 


3.3. The Residue Series—Non-Isotropic Medium 


Analysis similar to that of the last section can be used where eqns. (17) 
to (20) are in matrix form. We shall now assume that the transmitter 
is on the ground, so that H=0. Instead of (23), we then have 


S(C)= { exp (+jkOz)+R, exp (—jkCz)} exp (2jx)fexp (2ja)—R,}-1 
+ {exp (—jkCz)-+-R 5-1 exp (jkOz)}{ exp (2ja)—Ry}-1Rg. (28) 
This includes the reciprocal matrix 
OXP ede ae Rate me a 0 drs My (29) 
whose terms can be evaluated. (28) then becomes 


={ exp (+jkCz)+R, exp (—jkOz)}(1/4)M exp (2jx), . (30) 


where 
es ae (29%) —,R,,, ik, ; 31 
a yeas exp (2ja)-+ Rh, ve 
exp (27«a)+ ,R,, — Rk ) 

M=— et MN ee et Gy 
(Ce. exp (2ja)—,R, ae 

(30) has poles where 
en) Mee eae sence! 9b. .0'(33) 


There are an infinite number of these, and as before, the position of any 
one pole depends upon an integer. The residue at the pole at C=C, is 

{exp (—jkC,,z)+R, exp (gkC,2)}{1/A'(C,,) }M(C,,)Ro(C,,). - (84) 
To integrate eqn. (20) we close the contours as before (figs. le and 
f). The algebra involved in evaluating the matrix products is tedious, 
and we give below only the final result for the pole of order n. If m is 
negative, S,,, C,,, must be replaced by S,,*, C,,*. 


TTyn= —ThHiy(kS,,p) cos (kC, # [exp (2j«,,)+,f,] 
x exp (2ja,)P,/4'(C UA ace sy cae (SO) 
TT n=jrkHiy(kS,,p) sin (kC',z te (2jo%n) ty, Py/A(C eae oO) 
where 
A'(C,)={ exp (2ja,)-+,R, }{2jhk exp (2ja,)—,R,'} 
+ { exp (2jo,)—yR, }{2jhk exp (2ja,)+ Ry phy f+ yhy Ay (37) 


It is easily checked that (35) reduces to one term of (27) for the isotropic 
case. 

From eqns. (35) and (36) we may deduce the ratio of the horizontal 
components of the magnetic field in, and perpendicular to the direction 
of propagation, for the wave-guide mode of order n. From eqns. (9) 
and (11) this ratio is equal to 


—c{oll joz}{/{ol, jet}, -= . - 22+ as « (38) 
and using (35) and (36) it becomes 
{Oty Aye exP (2jRRO a) jn le ee (39) 
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3.4. The Curvature of the Earth 


So far we have assumed that the earth is flat, and that there is free 
space between the walls of the wave-guide. The curvature of the earth 
can be allowed for (e.g., Booker and Walkinshaw 1946) by considering 
the wave-guide to be filled with a fictitious medium having a vertical 
gradient of refractive index p, given by 


w=14+Ke Sg St 2h re 


where ¢ is measured vertically upwards from the base of the top wall 
of the guide, and is therefore negative within the guide, and K is the 
curvature of the earth. If Kh<1, we may incorporate this modification 
into the theory by using the phase integral method. We shall neglect 
the effect of the curvature on the ionized medium which forms the upper 
boundary of the guide, so that we assume that it remains homogeneous. 
Since the fictitious medium in the guide is not homogeneous, the com- 
ponent waves of any one mode now have curved wave-fronts. Let @ 
be the inclination to the vertical of the wave-normal of a component 
wave at a height z=h+¢. Let g=pcosé. If 06=7i at the surface of 
the ionosphere, then for any height in the guide 


g=p*—sin* 1 2KC- cost, = ee 
since K is small. 
.g== cosi+KC/ cost. se phe a> hs ae 


The phase integral relation is 
0 
exp (2) kq dC }=R5:. exp.(20n)) eee 
—h 


If K=0, the left hand side becomes exp (2jhkC) and the equation reduces 
to (22). If K-0, we have for the left hand side of (43), using (41): 


exp {2jhkC—jKkh?/C}, » at Dae = eae 

where C'=cos i. 
Hence the curvature of the earth may be allowed for by replacing 
exp (2jhkC) by (44) in the formulae. This is an approximate method, 


but is accurate as long as the second term of (44) is very much less than 
the first, that is, as long as 


| O | >\/(KA/2) 0 ee ee 


This was true for the numerical results discussed in §§ 4.4, 4.5 and 4.6. 


§4. APPLICATION TO THE IONOSPHERE 
4.1. The Reflection Coefficients of the Ionosphere 


The preceding theory is applicable to any model of the ionosphere, 
for which the reflection coefficients are known analytic functions of the 
angle of incidence, regarded as a complex variable. 
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Before we can proceed further, it is necessary to adopt some specific 
model for the ionosphere. Theories which take account of the rate at 
which ionospheric characteristics change with height (e.g. Wilkes 1947, 
Heading and Whipple 1952), do not give sufficiently simple analytic 
functions for this purpose. 

In Paper I the ionosphere was assumed to be a homogeneous sharply 
bounded medium containing N electrons per ¢.c., subject to collisions at 
an average rate v per second, but the earth’s magnetic field was neglected. 
A discussion has been given elsewhere (Budden 1951 b) of the reflection 
of very low frequency radio waves when the ionosphere is assumed to be 
homogeneous and sharply bounded, with the effect of the earth’s magnetic 
field included. In a full treatment of this problem the refractive index wu 
is given by the Appleton—Hartree magneto-ionic formula (Appleton 
1932). It was shown, however, that fair agreement with the experimental 
results at very low frequencies could be obtained using the quasi-longi- 
tudinal approximation to the magneto-ionic formula (Booker 1935) and 
that this approximation becomes more accurate as the frequency decreases. 
We are therefore encouraged to use it in discussing the propagation of 
radio atmospherics. We shall now apply it to the results of the preceding 
sections. 

Using this theory, we find (Budden, loc. cit.) that 


p2=1—(jp,[p) exp (biz), . 2 ee. (46) 


where p, and 7 are functions of the characteristics of the ionosphere, and 
do not include the wave frequency. It can be shown that approximately 


De AVe ym ees ae (47) 
UDTORGES Ue re ee ie pee ee eee Co) 
where 
por=4rNe*/eqm. 
v=collision frequency. 
Pr= oem. 
e, m=charge and mass of the electron. 
N=number of electrons per unit volume. 
H,=component of earth’s magnetic field in direction of propaga- 
tion. 
y=electric and magnetic permittivities of free space. (The 
symbol , is used later to mean the refractive index of the 
ordinary wave, and must not be confused with the py used 
here.) 

A plane wave incident on the ionosphere at angle of incidence i gives 
rise to a reflected wave below the ionosphere, and to both ordinary and 
extraordinary waves in the ionosphere. For each wave E, denotes the 
component of E in the plane of incidence and perpendicular to the wave 


€o> 
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normal, and H, denotes the component perpendicular to the plane of 
incidence. Then for upgoing waves in the northern hemisphere we 


have for the ordinary wave 


io? =1—(Jp,lp) exp (jz), » + ss ~ (49) 
E, of By o= — Js 2 ea 
and for the extraordinary wave 
u2=1—(9p,/p) exp (—Jr), = ses 
E ol Ewe td. © ee ee ee 


The wave normals of these two waves are assumed to make angles 4, 
and 6, respectively to the vertical. The angles i, 6, and 6, are related 
by Snell’s law. 

sin =p, sin jp, sin 0.) ee 
The expressions (49) to (53) were used, together with the boundary 
conditions which apply at the interface between the ionized medium and 
the free space below it (Budden, Joc. cit.), to find the components of R. 
The results apply both to #, and H, and to the components JJ, IT, of 
the Hertz vectors, and are as follows: 


y= {(Hol4e— 1)(cos 8,4 Cos 8,) cos 1+ (4,+p,)( cos? i— cos 8, cos 8,)}/D, 


(54) 
2 =2j COs 1(4, cos 0,—p, Cos 8,)/D, ° OS ToS Se ee ee 
2, =2)j cos a(n, cos O,—He cos 6,)/D, - oNcl S| Be er 
2, ={—(HoMe— 1)(cos 6,4 cos 4,) cos 1+(u,+p2,)( cos” i—cos 8, cos 8,) }/D, 

(57) 
where 


D=(p,p,+1)(cos 6, cos 8,) cos 1+ (u,+p,)( cos? 7-+cos 6, cos 8,). (58) 


4.2. Method of Obtaining Numerical Solutions 


In the following sections we give some results of calculations made 
using the preceding formulae. To find the characteristics of any one 
mode, it is necessary to find the positions of the poles of the integrand 
in eqn. (20) by solving eqns. (31), (33) using the expressions (46) to (58). 
Let 

hk =," 8 Te Sw gs ee 


Then eqn. (33) reduces to the following quadratic equation for tan « : 
Molt, tan® a—J tan «(14,4 4,)(0?+cos 8, cos 0,)/{C(cos 0, cos 6,)}—1=0. 
(60) 


It can be shown that one root of this quadratic gives the quasi-TM 
modes, and the other gives the quasi-TE modes. Suitable values of 
the parameters p,, 7, and h were selected as described in the next section. 
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The type and order of the mode for which solutions were required was 
also chosen. Then the equation was solved for each of a series of values 
of the angular frequency p. The method of solution was as follows. 
A trial value of C was assumed, and from it the values of cos 6, and 
cos #, were calculated, using (53). The quadratic was then solved for 
tan « from which « was calculated. Then (59) gave a revised value of 
C which was used to suggest a new trial value of C. The process was 
repeated until a self-consistent solution was obtained. The actual 
calculations were made using the EDSAC calculating machine in the 
Mathematical Laboratory of the University of Cambridge. 

In a few cases the curvature of the earth was allowed for, by using 
the equation 


a=hhO—KbhoC 8... (BL) 


instead of eqn. (59). Equation (61) is derived from (44). The effect 
of the earth’s curvature was found to be very small, and in most of the 
calculations the earth was assumed to be flat. 

It is also of interest to calculate, for any one mode, the ratio of the 
two horizontal components of the magnetic field of the wave, namely, 
the component in the direction of propagation, and the component at right 
angles to it. For an isotropic reflecting medium, this ratio is zero for 
a TM mode, and infinite for a TE mode. In general its value is given 
by (38), for which an algebraic expression was obtained. This was very 
complicated and will not be given here. We shall call this ratio the 
‘ polarization factor’. The results of the calculations are given in § 4.6. 

Finally we wish to known how much of a given mode is excited by a 
vertical electric dipole. For this purpose, a quantity was calculated 
which is proportional to the transverse horizontal component of the 
magnetic field of the wave, when the attenuation term is omitted. The 
expression used was derived from eqn. (35) using eqn. (9). It will be 
called the ‘excitation factor’. The results of the calculations are 
discussed in § 4.7. 


4.3. The Values of the Parameters 


The parameters p, p,, h, and 7 occur only in the combinations p/p,, 
h.p,, and 7. Solutions were obtained for fixed values of h . p, and 7, 
and for a series of values of p/p, (equivalent to a series of different 
frequencies, p). The frequency scales in the figures in this paper are 
chosen on the assumption that h=80 km. 

In a previous paper (Budden 1951 b), the results were given of some 
calculations of the reflection coefficients (54) and (55) for real angles of 
incidence, and these were compared with experimental measurements for 
a frequency of 16 ke/s (Bracewell et al. 1951). The values of the para- 
meters p, and + which gave best agreement with the experimental results 
were p,=2 10% sec!, r=60°. Most of the following results there- 


fore apply to these values. 
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In Paper I, some results were given of calculations in which the earth’s 
magnetic field was neglected, that is t=0, for several values of p,. All 
the modes in this case were pure TM modes. Some results for 7=0, 
p,=2% 10° sec"! are also given in this paper for comparison with the 
results when r=60°. They refer to TM modes of order —1, 0, and +1, 
and to TE modes of order —1 and +1. These TE modes cannot be 
excited by a vertical electric dipole, but it is interesting to compare them 
with the quasi-TE modes which occur if the top reflector is doubly 
refracting, and which can be excited by a vertical electric dipole. 


4.4. Numerical Results for Attenuation in Wave-Guide Modes 


In fig. 2, curves of attenuation versus frequency are shown for 
p,=2X 10° sec1, r=0°, h=80km. The full curves are for TM modes, 


Fig. 2 


Amplitude Reduction in 1000 km path. 


Frequency in ke/s.——> 


Fractional reduction in received e.m.f. for 1000 km path in modes of various 
order. h=80 km, p,=2 x 105 sec-!, r=0°. 


and the dotted curves for TE modes. In fig. 3, the corresponding curves 
are shown for p,=2x 10° sec-1, r=60°, h=80km: the full curves are 
for quasi-TM modes, and the dotted curves for quasi-TE modes. 

It is clear that the main features of the curves are the same both for 
7=0", and 7=60°. In particular, the quasi-TM mode of order zero is 
attenuated most heavily in a range of frequencies from about 1 ke/s to 
8 ke/s. This frequency range is determined mainly by the height of 
reflection, h, and only to a smaller extent by the characteristics of the top 
reflector. This Is illustrated in fig. 4, which shows the attenuation in the 
zero order mode for various values of p, when h=80 km and r=60°, 
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It was shown in Paper I, and is apparent from fig. 2 of this paper, that 
when 7=0, the TM modes of orders 0, and —1, have very similar character- 
istics. This similarity disappears when r=60°. Then the quasi-TM mode 

Fig. 3 
ie) 
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Amplitude Reduction in 1000 km path. 


@) 


Fractional reduction in received e.m.f. for 1000 km path in modes of various 
order. —f=80-km, p,-=2 < 10° sec—*, r==60°. 


Amplitude Reduction in 1000 km path. 


O 5 10 
Frequency kce/s. ——> 


Fractional reduction in received e.m.f. for 1000 km path in the zero order 
mode, as a function of frequency. 7=60°, h=80 km. 
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of order —1 is much more heavily attenuated, and its characteristics 
resemble more closely those of the quasi-TM mode of order +1. 

The attenuation was also calculated for the quasi-TM mode of order zero 
when the earth’s curvature was taken into account. It was found to be 
only very slightly greater (1:3% for a frequency of 5 ke/s), so that the 
effect of the curvature is very small. We should expect a larger attenu- 
ation because a wave which leaves the ground with a given angle of elevation 
is reflected from the ionosphere at an angle which is more nearly vertical 
if the earth is curved than if it is assumed to be flat. At nearly grazing 
incidence, the reflection coefficient decreases as the angle of incidence 
becomes more vertical (this is clearly shown in the curves given by Budden 
1951 b). Thus we should expect heavier attenuation when the earth is 


A Quasi us pers 


e) 
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Fig. 5 
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Magnetic Field. 
Scale for quasi-TM mode. 
Scale for quasi-TE mode. 


O 


Ratio: Longitudinal Magnetic Field/Transverse 


Polarization factor as a function of frequency, for quasi-TM mode of order 0 
and quasi-TE mode of order —l. p,=2X10® sec-1, r=60°, h=80 km. 


curved. The effect would become smaller if the angle of incidence were 
markedly different from grazing, and hence we should expect the earth’s 
curvature to have a still smaller effect on modes of order other than zero. 


4.5. Numerical Results for Polarization Factor 


Polarization factors were calculated as described in § 4.2 for the quasi- 
TM mode of order 0, and the quasi-TE mode of order —1. The polariz- 
ation factor is the ratio of the complex amplitudes of the horizontal 
components of the magnetic field in, and perpendicular to, the plane of 
incidence. The modulus of the polarization factor is plotted as a function 
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of frequency in fig. 5 for p,=2 x 105 sec“, r= 60°, h=80km. It is seen 
to be extremely small for the quasi-TM mode of order 0. It might be 
expected that the signal received from a distant low frequency transmitting 
station would consist predominantly of the quasi-TM mode of order zero, 
and hence we should expect that, if a rotatable loop receiving aerial were 
used, the bearing of the station would be sharp. This is known to be the 
case for signals from distant low frequency stations (see, for example, 
Round et al. 1925, Smith-Rose 1929). 

_ For the quasi-TE mode, the polarization factor is seen to be much larger. 
Hence, if a signal consisting predominantly of a quasi-TE mode were 
received on a rotatable loop receiving aerial, the orientation of the loop for 
maximum or minimum signal would bear no simple relation to the direction 
of the source. We shall see, in the next section, however, that the 
amplitude of a quasi-TE mode excited by a vertical electric dipole, is very 
small compared with the amplitude of the quasi-TM mode of order zero. 
It is therefore probable that quasi-TE modes play little part in the 
propagation to great distances of signals from very low frequency 
transmitters. 

Many experiments are known in which radio atmospherics are received 
on a pair of crossed loop receiving aerials. For example, the atmospheric 
direction finding chain of the British Meteorological Office uses a number of 
receiving stations of this type (see Adcock and Clarke 1947), in which the 
configuration of the horizontal magnetic field in the received signal is 
displayed on the screen of a cathode ray oscillograph. It is well known 
that most atmospherics received on a system of this type give a response 
which is practically a straight line on the screen. Hence the magnetic 
field of the wave is almost entirely in one direction which is assumed to be 
perpendicular to the direction of propagation, and a good bearing deter- 
mination is then possible on the source of the atmospheric. This suggests 
that quasi-TE modes, if present, have a very small amplitude. 
Occasionally, however, an atmospheric is received which gives a very 
elliptical or irregular response on the screen, and no accurate bearing 
measurement can be made. It is possible that this type of atmospheric 
contains an abnormally large proportion of quasi-TE modes. Since these 
would not normally be excited with large amplitudes by a vertical 
lightning flash, we must assume that the currents in the flash where they 
originate have a large horizontal component. 


4.6. Numerical Results for Excitation Factor 


Excitation factors were calculated as described in §4.2 for the 
quasi-TM mode of order 0, and the quasi-TE mode of order —1. The 
excitation factor is a complex quantity whose modulus is proportional to 
the amplitude of the transverse horizontal component of the magnetic 
field in the wave at a given distance from the source, when the attenuation 
term is omitted. The modulus of the excitation factor is plotted as a 
function of frequency in fig. 6, for p,=2 10° sec~!, 7=60°, h=80 km. 


4M2 
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The amplitude scale used for the quasi-TM mode of order 0 is 20 times 
the scale for the quasi-TE mode of order —1. From fig. 6 we see, for 
example, that for a frequency of 10 ke/s, the excited amplitude for the 
quasi-TM mode is 72 times the amplitude for the quasi-TE mode. For all 
frequencies above about 5 ke/s, the amplitude of the quasi-TE mode is 
very small compared with the amplitude of the quasi-TM mode. For lower 
frequencies the excitation factors become more nearly comparable, but 
in this range of frequency, the quasi-TE mode is heavily attenuated 
(see fig. 3). Hence, at a distance from a vertical dipole source, We should 
not expect quasi-TE modes to appear with appreciable amplitude. 


Fig. 6 
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Excitation factor as a function of frequency, for quasi-TM mode of order 0, 
and quasi-TE mode of order —1. p,=2> 105 see~!, t=60°, h=80 km. 


4.7. The Signal Received from an Atmospheric 


From the curves of excitation factor and attenuation as functions of 
frequency, it would be possible to calculate the signal received at a distance 
from a single atmospheric. This was done in Paper I, for the case when 
the earth’s magnetic field was neglected. The method involves the 
evaluation of a Fourier integral. The impulse responses of this type have 
not been calculated when allowance is made for the earth’s magnetic field, 
but since the attenuation curves are similar in form, whether or not the 
effect of the field is included, we should expect the results in the two cases 
to be very similar. 
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§ 5. CoNCLUSIONS 


The propagation of low frequency radio waves may be treated either in 
terms of successive reflections at the ionosphere, or in terms of wave guide 
modes. When the separation of the transmitter and receiver is large 
compared with the height of the ionosphere, it is simpler to use wave 
guide modes. The treatment given in Paper I, in which the earth’s 
magnetic field was neglected, led to a qualitative explanation of some of 
the phenomena associated with the propagation of radio atmospherics to 
great distances. 

The earth’s magnetic field is important in the theory of the ionosphere, 
and it was therefore necessary to investigate its effect on the wave-guide 
modes. This has been done in the present paper, and for the model 
discussed it has been shown that there is little change in the general form 
of the results, so that the treatment given in Paper I is justified. When the 
earth’s magnetic field is included, the modes are somewhat more compli- 
cated, and a method of calculating their polarization characteristics and 
degree of excitation has been given. 

The modes may also be affected by the curvature of the earth, which 
was neglected in Paper I. The present paper shows that this was probably 
justified, for frequencies below about 10ke/s. At higher frequencies, 
however, it is possible that the earth’s curvature may have an appreciable 
effect on the modes. 

In some recent experiments made by Dr. Weekes, the signal from a low 
frequency radio transmitter (frequency 16 ke/s), was measured for a range 
of distances extending up to 3650 km, from the transmitter. It is believed 
that the results of these measurements may be simply interpreted in terms. 
of wave-guide modes. 
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ABSTRACT 


Data for the determination of the mean lifetime of charged V-particles 
are discussed, following the procedure outlined by Wilson and Butler 
(1952). The tracks of twenty-seven V+-particles have now been photo- 
graphed on the Pic-du-Midi (2 867 m) by Armenteros et al. (1952) and by 
the authors in a circular cloud chamber, 28 cm in diameter. It is shown 
_ that, using these data, the chamber is too small for the lifetime to be 
determined, but a lower limit of 1:0 10-19 sec is obtained. Furthermore, 
the relation between the frequency of production of energetic V+-particles 
in penetrating showers of average energy (10-30) Bev, and various assumed 
lifetimes is discussed. Using recent photographic emulsion data on the 
frequency of particles of mass about 1200 m, in the same type of pene- 
trating shower an upper limit for the lifetime of 10~§ sec is obtained. 


§ 1. INTRODUCTION 


In a recent letter Wilson and Butler (1952) discussed several methods for 
measurement of the lifetime of unstable particles, such as charged 
V-particles. They concluded that the most direct procedure is to deduce 
a value from cloud-chamber photographs showing recognizable decays of 
V+-particles. The decays of twenty-seven V+-particles have been photo- 
graphed on the Pic-du-Midi (2867 m), using a circular cloud chamber 
28 cm in diameter. The momenta and angles of apparent deflection for 
twenty-two of these decays were given by Armenteros et al. (1952) in 
table 1 of their paper. A further five examples have now been obtained 
by the authors ; this additional data is summarized in table 1. 


Table 1. Measurements on Five Charged V-Particlest 
(1) (2) (3) (4) (9) (6) (7) 


; ; Sign and Transverse 

Event SE Me E/Tmin Angle momentum of momentum 
No of momentum of i) secondary ee 

: V (100 Mev/c) primary (deg) (100 mev/c) (100 mev/c) 

23 + — 1 13-5 +11-442-3 2-60-+0-50 
24 g — 1 15-20 — — 

25 — — 1 18 —6:0+1:0 1-86-40-31 

26 — -— 1 10 —16-+5 2:70 +0-90 

27 oe — 1 19 +8-0+1-6 2-64-L0-53 


* Communicated by Professor P. M. 8. Blackett, F.R.S. 
+ This table has the same form as table 1 of Armenteros ¢¢ al. (1952). 
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§ 2. Tue ExperIMEnTAL Data 


In order to determine the lifetime of V*-particles it is necessary to 
know two quantities for each decay in addition to those given in table fh 
Firstly, the length of the track produced by the unstable particle in the 
illuminated region of the chamber, and secondly, the distance the particle 
would have travelled in the chamber if it had not decayed must be 
measured. Wilson and Butler called these two lengths the ‘ path length ’ 
and the ‘ potential path length ’ respectively. Thus decay takes place after 
a path length, a, which is always shorter than the potential path length, 
a). Twenty-three pairs of values of x and a, are given in columns (2) 
and (3) of table 2. The tracks of the remaining four V*-particles are 
either too short for any measurements to be made or they are situated 
near the edge of the chamber in regions where stereoscopic measurements 
cannot be made. 

The values of a and a, have been determined by reprojecting the stereo-— 
scopic photographs. In most cases a is easily measured, the estimated — 
errors are given in column (2) of table 2; only approximate values can be 
given for a few of the events. The charged V-particles 3 and 10 were 
produced in a lead plate across the centre of the chamber ; 2 cannot be 
measured accurately for 3 because the decay occurred among a great 
many penetrating shower particles, the tracks of which prevent accurate 
reprojection of the decay. The lead plate was also in use when events 1, 
2 and 8 were obtained. When either of the lengths x or a, crossed the 
plate its thickness has been subtracted from the measured values. 

In order to measure a the track of each unstable particle has to be 
extrapolated until it leaves the illuminated region of the chamber. There 
is little difficulty in finding where the charged V-particle entered the fully 
illuminated volume of the chamber, but there is often some uncertainty 
in the determination of the exit point, supposing that the particle had not 
decayed. The illuminated depth of the chamber is about 7 cm; it is 
important, however, to determined this depth for each photograph because 
the illumination conditions may vary during the operation of the cloud 
chamber. Most of the decays of V*-particles are associated with showers 
of several particles, and the illuminated depth can be determined using these 
particles. If this procedure is not convenient, penetrating showers on the 
same piece of film are used. The estimated errors in the values of a are 
given in column (3) of table 2; in several cases it is impossible to give 
reasonably accurate values of 2». 

Very few of the observed V+-particle decays took place very close to 
the top of the cloud chamber. This is probably due to the difficulty of 
identifying decays which take place near the top of the chamber (or 
immediately underneath the lead plate) because of the confusion produced 
by the tracks of other particles in the penetrating showers. There is, 
therefore, a bias against selecting photographs showing decays with small 
«. Wilson and Butler pointed out this difficulty and suggested that 
the values of x and 2) should be measured from fiducial surfaces, draw 
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a predetermined distance within the chamber. For the present data 
it is convenient to measure x and a, from the surfaces of the well- 
illuminated volume and then to reduce the values of a and 2» (%) must 
also be reduced since it is measured from the beginning of the path 
length) by a constant amount. This length is sufficient for a decay to 
be identified and also for the direction of the charged V-particle to be 
found with reasonable precision. The values of « have all been reduced 
by 3cm; this figure is somewhat arbitrary and is discussed again in $3. 
Four decays with values of x less than 3 cm are rejected from the analysis. 
The corrected values of x, namely x’, are given in column (3) of table 2. 

Similarly, it is very difficult to detect the decay of a charged V-particle 
very close to the end of its potential path length. To avoid the possibility 
of a bias against observing this type of event all the values of x» are 
reduced by lem. This distance should be long enough for the apparent 
deflection of a track to be recognized, and thus the event identified. 
The corrected values of x), given as x» in column (5) of table 2, are 
obtained by reducing the values of x) by 4 cm. 

To measure the lifetime of the unstable particles it is necessary to 
convert the distances x’ and x,’ into the proper times of flight. The 
momenta and masses of the unstable particles must be known before this 
transformation can be made. Provisionally it is assumed that the 
charged V-particles are all of the same kind and have a mass of approxi- 
mately 1200 m,; this assumption is discussed in § 3. 

The momenta, P, of nine of the V+-particles have been measured ; 
the values are given in column (7) of table 1. In addition, three 
approximate lower limits are given for three high-energy events. 
Photographs 3 and 8 show heavily-ionizing primary particles, and the 
momenta given for these have been determined from the estimated 
ionizations and the assumed mass. Rough estimates of the momenta 
of the remaining nine events, which all have short primary tracks, can 
be made because the average value of the ratio of the primary to the 
secondary momentum, for the events with measured primaries, is very 
approximately 2-5. This rule was used when making the estimates 
labelled with an asterisk in column (7) of table 2. 

The values of the dilatation factor, P/M, are given in column (9) of 
table 2. P is the momentum and M is the rest mass, in energy units, 
of each charged V-particle. The proper times of flight, ¢, and the 
potential traversal times, 7’, are given in columns (10) and (11) of table 2. 


These values are obtained by dividing a’ and a,’ by Pc/M, where c is the 
velocity of light. 


§ 3. THe Mnan Liretime oF CHargep V-PAaRTICLES 


In principle, the mean lifetime of the charged V-particles can now be 
calculated from the pairs of times ¢ and 7’. These times have been 
measured assuming that all the V+-particles have a mass of 1200 mM,» 
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Armenteros et al. (1952) concluded that V+-particles are probably similar 
to the S-particles found by Bridge and Annis (1951) and Annis et al. 
(1952) and also to the «-particles found by O’Ceallaigh (1951). The S- 
and «-particles have a mass of about 1200 m, and decay at rest into a 
single charged meson. Recent work at Bristol, reported by Powell (1952), 
shows that the 1200 m, mass mesons observed in photographic emulsions 
are not all the same. There now appear to be at least two types, with 
masses of about 1100 m, and 1500m,. It is probable that the V+- 
particles are also not homogeneous, and therefore it is desirable to classify 
the particles responsible for the observed decays before proceeding to 
the lifetime determinations. Unfortunately, at present, the decays-in- 
flight observed in the cloud chamber cannot be classified unambiguously 
and, therefore, the data can only be analysed on the assumption that all 
the V+-particles have an average mass of about 1200m,. Owing to the 
small number of decays available, no evidence for the existence of two 
types of V+-particle can be obtained from the lifetime measurements. 
Unfortunately there is no possibility of measuring the separate lifetimes 
of different types if they occur in about equal numbers, even if they have 
appreciably different lifetimes. 

A preliminary examination of the usefulness of the data given in 
table 2 can be made by using the ratios w’/x)’ given in column (6). 
It can easily be shown that if the distribution of a large number of values 
of xv’/x,' is almost uniform from 0 to 1, with an average value close to 0-5, 
then the lifetime could be much longer than the mean traversal time, and 
no significant upper limit can be obtained. In fact the average value of 
x’ /x,' for eighteen decays is 0-38 and the mean value of 7’ is 2-6 x 10~1° sec. 
There are two sources of error in this determination of the average 
x'/x,'. Firstly, the error due to the individual errors of measurement 
of x and a,’ is 0-03; this approximate value is calculated on the 
assumption that the errors in columns (2) and (3) of table 2 are independent. 
Secondly, there is a larger statistical error due to the fact that only a small 
number of ratios is available; an estimate of this error is 0:07. The 
average ratio, therefore, is not significantly smaller than 0-5, consequently 
the lifetime could be very much longer than 2-6 x 10-1°sec, and no upper 
limit to the value can be calculated. In order to obtain a lifetime using 
a sample of about 20 events, the average ratio x’/x,’ needs to be less 
than about 0-25. 

Owing to the observed form of the momentum spectrum of the 
V+-particles and to the geometry of the cloud chamber the values of 7 
are spread from 0-6 10-1 sec to 8-010"! sec. It is obvious that the 
most important decays are those with the longer potential traversal 
times. Photographs 8, 11 and 22 have an average potential traversal 
time of 7< 10-19 sec, and the average of their ratios w/a,’ is 0-54. Thus 
there is no possibility of making a lifetime determination by considering 
the few important events alone. In future work it is clearly necessary 
to improve the selection of decays with long values of 7’. 
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The average value of «x’/x,’ is influenced by the correction procedure 
used when measuring a’ and 2)’. That the average of 2'/ay’ depends on 
the choice of the correction length is shown in table 3, where the correction 
procedure is indicated in column (1), and the corresponding average ratio 
is given in column (2). The amount of data diminishes as the correction 
lengths are increased, as shown in column (3). 


Table 3. The Influence of the Correction Procedure 


(1) (2) (3) 
Correction Average No. of 
procedure ratio decays 

Pull #375 0-46 19 
x —3 cm, x, —4 cm 0-38 18 
«x —5 em, % —6 cm 0-38 14 


The average ratios in column (2) become approximately constant after 
3cm have been subtracted from each value of x and 4cm from each 2p, 
and thus the bias effects are probably eliminated by these reductions. 

A statistical treatment of lifetime data in the form of values of ¢, within 
a potental traversal time 7’, was given by Peierls (1936). He showed 
that a minimum value of the lifetime, independent of 7’, is given by the 
arithmetic mean of the values of t. The average of t, calculated from the 
data given in column (9) of table 1, is 1:0 10-1!®sec. The statistical 
theory given by Peierls was discussed by Hole (1947) and a maximum 
likelihood procedure for finding the most probable lifetime from a set 
of values of t and 7' was given by Bartlett (1936). The available data on 
V+-particles is not sufficiently extensive to warrant a calculation of © 
this type. 


$4. Tor FREQUENCY OF PRODUCTION OF CHARGED V-PARTICLES 


Assuming that the V+-particles observed in the cloud chamber are 
predominantly of the same type it is possible to find an approximate 
relationship between the frequency of production and various assumed 
lifetimes. An idealized version of the Pic-du-Midi apparatus, shown in 
fig. 1, is used and the following simplifying assumptions are made. The 
production of V+-particles is assumed to occur uniformly throughout 
the lead absorbers, and the fraction of the shower particles which consists 
of V*-particles is regarded as independent of the energy of the primary 
particle. Attention is restricted to shower particles with momenta greater 
than 10° ev/c and with inclinations to the vertical of less than 30°. The 
inclinations of the shower particles in the producing layer are neglected, 
and it is assumed that they have a path length in the cloud chamber 
equal to the average value of 15 em. Furthermore, in this high momentum 
range, it can be assumed that the momentum spectrum follows an 
inverse power law. A negative exponent of three is used for the differ- 


ential momentum spectrum of all the shower particles (see Barker and 
. Butler 1951), 
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The ratio of the number, Ny, of V+-particles produced to the total 
number of shower particles, N,, is F=N,/N,, which can be written 
(N/N,)/(N/Ny), where N is the number of V-particles decaying in the 
chamber. There were 8000 (+25%) tracks of shower particles and | 
14 V-particle decays, satisfying the requirements laid down in the initial 
assumptions. Thus N/N, is equal to 1-75 10-3. The ratio N/Ny can 
be calculated as a function of 75, the assumed lifetime of the V-particles 
viz.: 


N fe [ [exp—a(x-+-b)][1—exp (—ac)|P-3 dP da 


Ny Sc) 


{ | PeaPax 
Py’ 0 


Fig. 1 


+—x—en 


— Nn —— >} OC -2q——_ 0 ———o 


Simplified arrangement of cloud chamber and lead absorber 


where «=M/P7,, M and P being the mass and momentum of a V-particle. 
The integrations are performed over the thickness of the lead absorbers, 
a, and for particles with momenta greater than Py. The relationship 
between F and 7, is shown in fig. 2 for the following conditions : 
P,=109 ev/c, a=15 cm, c=15 cm, b=5 or 10 cm, and M=1 200 m,. 

The various simplifying assumptions that have been made in these 
calculations are not very critical; for example, increasing the gap 
between the absorbers and the chamber from 5 cm to 10 cm does not make 
a very large difference to the frequency-lifetime curve. 

The estimated average energy of the primary particles incident on the 
absorbing layer which trigger the counter detection system is 20-40 Bev. 
The average energy, however, of all interactions contributing to the 
flux of shower particles with momenta greater than 10° ev/c is probably 
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somewhat lower than this value since a considerable number of secondary 
interactions occur in the absorber. Thus the curve in fig. 2 refers to 
V+-particles produced in showers of average energy about (10-30) Bev. 
If the V-particles were to interact in the lead above the chamber with 
the geometric cross section the values of Ny/J, in fig. 2 would be increased 
by about a factor of two. 

Data obtained in photographic emulsion work (Perkins 1952) make it 
likely that particles of mass about 1200 m, constitute less than 10°%,°ot mm 
the total number of shower particles in penetrating showers of average 
energy about 15 Bev. Using this value, an upper limit of about 10~* sec 
is found for the mean lifetime of the charged V-particles which are 
observed to decay in the cloud chamber. The minimum frequency of 
production of about 1° corresponds to a mean lifetime in the range 
10—1°_]0-® sec. 


Fig. 2 


' =O -9 =8 at 


lO 1O lO 
To (sec) 


The relations between the lifetime and frequency of production of V+-particles. 
(i) Curves A and A’ show Ny/Ns as a function of 7, for distances between 
___ the absorber and the chamber of 5 and 10 em respectively. 
(ii) Curve B shows N/Ny as a function of 7, for a distance of 5 em between 
the absorber and the chamber. 


§ 5. CONCLUSIONS 


Under the assumption that the charged V-particles decaying in flight 
in the cloud chamber are predominantly of one type and similar to the 
«-particles observed in photographic emulsions (Armenteros et al. 1952) 
limits can be placed on the mean lifetime of the particles. 

Following the suggestions of Wilson and Butler (1952) data concerning 
the distribution across the chamber of the decays in flight of eighteen 
V+-particles, observed in a chamber 28cm in diameter, have been 


? 
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presented. A minimum value of 10-19 sec hag been found for the mean 
life. Considerably more data are required before a precise estimate can 
be made. 

The frequency of production of V*-particles has been considered as a 
function of various assumed lifetimes. From the recent photographic 
plate experiments, which show that «-particles of mass about 1200 Me, 
constitute less than 10% of the total number of shower particles, an 
upper limit of 10~* sec is found for the mean lifetime. 

These results may be compared with those obtained by other workers. 
Daniel et al. (1952) suggest that the lifetime of «-particles is about 
10-°sec. However, estimates obtained from cloud-chamber investigations 
of slow unstable particles yield a lower limit of 10~® sec (Astbury et al. 
1952, Annis e¢ al. 1952). The results obtained in the present work, 
concerning the energetic V+-particles produced in interactions of average 
energy 10-30 Bev, are consistent with either of these values. 
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[Received June 10, revised July 30, 1952] 


ABSTRACT 


Etch figures similar in appearance to those observed in the mineralogical 
technique are shown to result from the rapid chemical dissolution of 
growth spirals on silicon carbide. Each spiral resulting from an inde- 
pendent screw dislocation gives an etch figure. The etch figure produced 
by rapid dissolution, therefore, becomes a convenient means for indexing 
screw dislocations whether growth spirals are evidenced or not. Prolonged 
rapid dissolution of a crystal with a screw dislocation extending between 
opposite crystal faces results in a hole through the crystal at the site of 
the screw dislocation. 


SMALL round holes in crystals of AlB, formed in an aluminium melt and 
recovered by dissolving away the melt with strong acid have been 
described (Horn, Kasper and Fullam 1952). Spiral growth was observed 
on the AIB, crystals prior to the formation of the hole. Treatment of 
the crystals in strong acid was found a necessary part of producing the 
hole. The site of the hole is the centre of the spiral. Since the time 
F. C. Frank originally proposed that crystal growth may be initiated by 
screw dislocations (Frank 1949), A. R. Verma (1951) among others has 
successfully demonstrated that spiral growth patterns are evidence for 
growth by a screw dislocation mechanism. The holes in AIB, crystals 
have, therefore, been attributed to the preferential chemical attack at a 
screw dislocation. 

Silicon carbide has been sevestipated more thoroughly for spiral growth 
than other substances. It seemed desirable to extend the study of the | 
rapid dissolution of growth spirals to include those on SiC. Consideration 
of the process of dissolution of growth spirals leads to the conclusion that 
in order to attack the region of the screw dislocation preferentially, the 
process must be rapid. Dissolution near equilibrium conditions will 
cause a general disintegration of the crystal. The difference is readily 
demonstrated in the case of the AIB, crystals depending on the concen- 
tration of acid used. The emphasis of these studies is on the results from 
rapid dissolution. 

Silicon carbide is attacked by alkali carbonate fusion (Scott 1939). 
If the reaction is carried out in air at about 1000° c, there appears to be 
a uniform attack of silicon and carbon atoms. At lower temperatures, 
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a carbon residue may adhere to the crystals. The details of the results 
obtained depend somewhat on the carbonates used. However, the 
conclusions drawn from the present studies are sufficiently qualitative 
that they are not influenced by the detail of the etching process beyond 
knowing that the etching has been rapid. 

After selecting silicon carbide crystals with the desired spiral growth 
patterns, the procedure for following the dissolution of the growth patterns 
consisted of making photographs of replicas prepared from the crystal 
surface after each stage of alkali carbonate fusion. 

In Plate LX XITI there are shown photographs of replicas of the surface of 
a SiC crystal selected because the entire surface (a) was marked by terraces 
resulting from the interaction of eight spirals, (b) of the same sense. 
Using the criterion that the step width per revolution for each spiral is 
constant, each spiral may be considered to be derived from an independent 
dislocation. In the photographs (c) through (f) there are shown replicas 
of the same area as in (b) after successive treatments of the crystal in a 
fused alkali carbonate mixture at approximately 1000°c for the total 
times indicated. After two minutes, the spirals are obliterated. The 
centre of each spiral becomes indexed by a very fine pit, which on 
successive treatments enlarges and deepens. The underlying hexagonal 
symmetry is eventually apparent from the rays within the etched figures. 

Although one must keep in mind that the detailed structure of etched 
figures depends on the etching conditions, it is seen from the above 
results that the centre of each isolated spiral arising from an independent 
screw dislocation is through rapid etching replaced by a pit or etch figure. 
Thus, the etch figure may, under such conditions, be used to identify 
the centre of a spiral or more particularly, to identify the position of a 
serew dislocation. 

From the results obtained with independent dislocation spirals, one 
may turn to a crystal known to be made up from many spirals, or 
multiple interacting dislocations. Plate LX XIV (a) shows a photograph 
of a replica of the surface of a SiC crystal seen to exhibit several areas and 
kinds of growth spirals. In (c) is shown enlarged the spiral to the right 
in (a). It appears to begin from a single dislocation and after several 
turns, merges with many small spirals to one side. In (b) is shown 
enlarged the area to the left in (a). This area contains two spirals compos- 
ed of the interaction of a number of apparently closely spaced dislocations. 
Photographs in Plate LX XV show the crystal area in Plate LX XIV (a) 
after treatment in fused Na,CO, for the times indicated. 

Several features may be noted from these results. As previously, 
the centre of the independent dislocation spiral is replaced by an etch pit. 
The centre of each complex spiral is also replaced by an etch pit. In 
addition, however, there are many new etch pits—etch pits for which 
there was not a visible corresponding surface spiral. Another striking 
feature of the etch pits is their arrangement. These are, in general, 
arranged in rows between which are the 30-, 60-, 120- and 90-degree angles 
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characteristic of the underlying hexagonal symmetry of the crystal in 
the plane of the face. One may only speculate concerning the significance 
of the arrangement of the etch figures. These linear arrays may be 
evidence of dislocations caused by deformation of the crystal during 
growth. Such deformation has been proposed by F. C. Frank (1951) — 
and is used by A. J. Forty (1952) to account for the large step heights 
encountered in spiral growth on CdI,. Alternatively, the rows of figures 
may be connected with adsorbed impurity or dislocations formed at the 
edges of the crystal during growth. Another feature of the etch figures, 
seen particularly in Plate LX XV (a), (b), and (c) is that they are arranged 
approximately evenly spaced in the rows, rather as though a repulsive 
force existed between them. 

One may now turn to see what happens upon etching a crystal with 
no visible growth spirals on the surface, but which is obtained from the 
same group of crystals among which there are those showing spiral growth. 
The photographs in Plate LX XVI are of the surface of SiC crystals etched 
as indicated. Although no spirals were detected on the surfaces of these 
crystals, etching produces a large number of etch figures similar in 
appearance to those observed from crystals with spirals. One derives 
the impression from looking at many of these crystals that there are, 
in fact, many more etch figures obtained from crystals with no observable 
spirals than from crystals on which spiral growth is observed. This is 
consistent with the idea that such crystals may have smooth surfaces 
because of the interaction of the growth steps from a large number of 
screw dislocations. Similarly to the previous two situations, the etch 
figures are arranged in rows which display the hexagonal symmetry in 
the face plane and being approximately evenly spaced in the rows suggest 
a repulsive force between dislocations. 

It is of interest now to see what happens upon prolonged etching. 
Though more complex situations are geometrically possible, it is expected 
that frequently a screw dislocation extends perpendicularly through from 
face to face of a thin crystal. Then, upon prolonged etching, a hole may 
be produced through a crystal. This process was proposed for the holes 
in AIB,. That the same is true for SiC is seen in the photographs in 
Plate LX XVII which show holes produced in SiC by alkali-carbonate fusion. 
For preparing photograph (c), gold was evaporated on the surface of the 
crystal in order to improve the contrast. When small, the holes are 
round ; when etched to a large size, the hexagonal symmetry becomes 
apparent. 

The results reported here appear to have much in common with the 
technique of producing etch figures for interpreting the symmetry of 
mineralogical specimens (Honess 1927), thus prompting some speculation. 
In common is the fact that the results depend on the etchant and the 
rate of etching, that in both cases the crystal symmetry is revealed by the 
figure developed, and that the figures themselves generally are arranged 
according to the crystal symmetry. The deeply etched rows of screw 
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dislocation. sites have the appearance of ‘solution channels’ often 
encountered in the mineralogical technique. Thus, one is led to think 
that the origin of many etch figures on natural minerals may be due to 
the presence of screw dislocations in the mineral. A corollary to this 
point of view is that natural minerals should be a source of materials 
showing spiral growth. It was on natural beryl that spiral growth was 
first recognized (Griffin 1951). Spiral growth has also been found on 
natural crystals of barite (Horn, unpublished), graphite (Horn 1952), 
and micas (Amelinckx 1952). A reinvestigation of minerals from the 
viewpoint of screw dislocations may be rewarding in understanding the 
properties of screw dislocations and, also, in providing additional evidence 
for the growth conditions in various deposits. 

Although the effect of the results reported here is to raise rather than 
to answer problems, the study is reported now because there is implied in 
the results a technique for revealing screw dislocations in crystals. 
The technique depends on the differential rates of chemical attack of a 
screw dislocation and good crystal when the etching is rapid. The 
technique appears to be applicable whether there is evidence of screw 
dislocations from spiral growth patterns or not and may, therefore, be 
presumed to apply to revealing screw dislocations emerging from a surface 
whether the dislocations are introduced subsequent to or during crystal 
growth. 


The author is pleased to acknowledge the benefits from discussion with 
and the interest of Dr. J. 8S. Kasper, Dr. M. H. Hebb, and Dr. L. Apker 
of this laboratory. The criticisms of Dr. F. C. Frank during the prepar- 
ation of the manuscript are also gratefully acknowledged. 
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Angular Correlations in the Reaction '°N(p, ay)?2C 


By J. Srep and A. P. FRENcH 
Cavendish Laboratory, Cambridge* 


[Received August 25, 1952] 


Tum bombardment of 5N by protons has been studied by Schardt eé al. 
(1952), who found several resonant states of 1*O for which the chief mode 
of decay is to the 4-47 mev level of #C by emission of low energy alpha- 
particles. 

We have studied the angular correlation between these alpha-particles 
and the subsequent gamma-rays at the 898 kev resonance, which is 
sufficiently strong and isolated to be regarded as a pure state. 

Targets of 1=N were obtained from the Atomic Energy Research 
Establishment, Harwell; they were formed by bombarding thin 
aluminium sheets in the electromagnetic separator. An aluminium 
backing was chosen in order to reduce the number of elastically scattered 
protons from the target. 

In our experiments, alpha-particles emitted at 90° to the incident 
proton beam were separated from the scattered protons by magnetic 
analysis, and were detected by a thin ZnS screen and E.M.I. photo- 
multiplier. The gamma-rays emitted in the reaction were detected by 
a 1l}in.x14in. mosaic of Nal crystals and a photomultiplier. The 
crystals were situated some 2 in. from the target, giving a solid angle 
of about % sphere, and their detection efficiency was about 30% for 
the 4:47 Mev gamma-rays. By rotation of the gamma-phosphor about 
the target, the («, y) correlation could be studied over the angular range 
90° to 180°. The arrangement was such that an («, y) correlation could 
be studied in a plane either containing or perpendicular to the incident 
proton direction. 

It is intended to give a fuller description of the apparatus elsewhere. 

Both the alpha and gamma single counting rates were continuously 
monitored, and the number of coincidences was obtained as a fraction 
of the number of alpha-particles counted at a given angle. The 
coincidence resolving time was 0-3 psec. 

The experimental points are shown in fig. 1, and it is seen that a high 
order of anisotropy is present in the («, y) correlations in both planes. 
The results can be fitted if one assumes that the incoming protons have 
orbital momentum 2 and that the compound nucleus 160* is (2—) (ie. 


has spin 2 and odd parity). There is evidence to show that the first 
ee ee eee 
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excited state of ?C is (2-+-) (Lewis 1952, Hubbard et al. 1952, Kraus and 
French 1952); it could thus be formed from the 160* by emission of 
alpha-particles with J=1. 
On the above assumption, the theoretical correlation functions are 
7 6 
9) en eee ee 2 
de (0.5) 1 ie 6 

for the («, y) correlation in a plane perpendicular to the incident proton 
direction and 
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, y) angular correlations in 4®N(p, ay) at 898 kev: (a) in plane containing” 
: eas beam ; (b) in plane perpendicular to beam. The alpha particles. 
were observed at 90° to the beam in each case. Errors shown are standard 


deviations. 
for the correlation in a plane containing the incident proton direction. 
In these expressions, f is the relative probability of forming the ao 
compound nucleus from the parallel and anti-parallel spin configurations. 
respectively of the colliding particles. 
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The full curves of fig. 1 are obtained by putting f=0-72 and correcting 
for the finite geometry of the apparatus. The value of f is deduced from 
the observed angular distribution of the gamma-rays relative to the 
proton beam at 898 kev (Kraus and French 1952), and it can be seen 
that an adequate fit to the experimental points is obtained. One should 
note that alpha-particles with 1=3, as well as /=1, might be emitted 
from a (2—) state of 1*O*; it was not, however, found necessary to 
consider these in fitting either the (p,y) distribution or the («, y) 
correlations. 


We are grateful to the staff of the Atomic Energy Research Establish- 
ment, Harwell, for providing separated nitrogen isotope targets. One of 
us (J. 8.) is indebted to the Department of Scientific and Industrial 
Research for a grant. 
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An Effect of Electron Bombardment on Order in Cu,Au Alloy 


By J. Apam, A. GREEN and R. A. DUGDALE 
Atomic Energy Research Establishment, Harwell, Berks.* 


[Received September 8, 1952] 


‘THE experiment described in this letter was designed in order to study 
the part played by lattice defects, introduced by high energy electron 
bombardment, in the ordering of a Cu,Au alloy. 

The effect of a pile irradiation on the electrical resistance and unit-cell 
dimensions of an ordered Cu,Au alloy containing 25-8 atomic per cent Au 
was studied recently by Adam and Dugdale (1951). Their main con- 
clusion was that, under suitable conditions, neutron bombardment is 
capable of producing simultaneous ordering and disordering effects. 
Using an initially disordered Cu,Au resistor, Blewitt and Coltman (1952) 
found that the relaxation rate for ordering at 200° c was considerably 
increased by neutron bombardment in a pile. They attributed this 
effect to the accelerated growth of ordered domains by the diffusion of 

Frenkel effects formed by the bombardment. 
In the present experiment a resistor of the 25:8°% alloy was prepared 
in the form of a strip 0-0015 in. thick. After annealing to remove the 
effects of cold work a high degree of order was induced in it by heating 


a 
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for six days at 381°c (which is about 7°c below the measured critical 
temperature), followed by slow cooling to room temperature at 20° c 
per hour. After this treatment the resistivity at 0°c was about 
45 ohms cm™1. X-ray diffraction photographs of similarly treated 
specimens show strong and sharp superlattice lines which indicate a 
high degree of order. The specimen was then bombarded for 2 hours with 
1 ev electrons, the integrated flux being about 3x 1017 electrons/sq. cm, 
the specimen temperature varying between 40°c and 45°c. The 
bombardment caused no measurable change in resistance. The specimen 
was then heated at 100°c for various periods of time; its resistance 
{measured at 0° c) decreased as shown in the graph. Subsequent heating 
at 130° c resulted in a further decrease in resistance. Further electron 
bombardment of the specimen followed by heating at 100° 0, again 
produced a marked decrease in resistance, and a third irradiation resulted 
in a similar effect. Preliminary experiments with 3 Mev electron 
bombardments gave similar results. 


(6) 100 200 
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Fractional change in resistance v. annealing time, at 100° c and 130° ¢, for an 
ordered Cu, Au resistor after electron bombardment (1 Mev). 


On the other hand, heating the specimen for 3 hours at 100° c before 
the first bombardment produced no measurable change in resistance. 

1 mev electrons: are capable of giving up to 68 ev recoil energy to a 
Cu atom, and up to 22 ev to a Au atom in elastic collisions with the 
nuclei. Seitz (1949) has suggested a nominal figure of 25 ev for the 
energy required to form an interstitial atom by bombardment. Assuming 
this to be the case in Cu,Au, the electrons are capable of knocking at 
least the Cu atoms out of their normal sites into interstitial positions, 
forming the so-called Frenkel defects. The cross section for 25 ev or 
more to be given to a Cu atom is of the order of 20 barns (estimated from 
the formula given by Mott and Massey 1949). Thus, the electron 
bombardment may be expected to have created Cu Frenkel defects with 
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a concentration of about 6 per million. If one or other of the components 
of these defects—viz. interstitial atoms or atomic vacancies—acquired 
sufficient mobility at 100° © to diffuse some distance through the lattice 
before being captured, movement of some Cu and Au atoms from 
‘wrong’ to ‘right’ positions could take place. The increase in order 
thus produced would explain the observed decrease in resistance. 

According to the theory of Sykes and Jones (1936) the specimen here 
used would initially contain large domains within which the equilibrium 
degree of order would correspond to a temperature in the neighbourhood 
of 250°c. When large domains are present the electrical resistance is 
not very dependent on domain size and it would seem most probable 
therefore that the diffusing lattice defects increased the order within 
the domains. The experiment described here suggests that by introduction 
of artificially produced lattice defects it is possible to reach an equilibrium 
degree of order corresponding to a temperature lower than that obtainable 
by purely thermal treatment in a reasonable time. 


The authors’ thanks are due to Messrs. T. M. Fry and J. H. W. Simmons. 
for stimulating discussions on this subject, to Messrs. E. R. Wiblin and 
J. D. Milne and members of the Linear Accelerator and Van der Graaf 
sections for electron irradiations, and to the Director, Atomic Energy 
Research Establishment, for permission to publish this letter. 
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Critical Shear Stress and Temperature 
By E. N. pa C, ANDRADE, F.R.S.* 


| Received September 8, 1952] 


Ir is generally accepted that, as regards the mechanical behaviour of 
single crystals, there is a contrast between the metals of hexagonal 
structure which, in general, possess one set of glide planes, and those, 
like the cubic metals, with a multiplicity of glide planes. This view fs 
supported by consideration of the curves showing the variation of 
resolved shear stress with glide over a wide range of glide. The diagram 
originally prepared by E. Schmid (1934) clearly separates the two classes 
the cubic metals showing very much greater hardening with glide. It % 
———— $a eee 
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to be noted that the group of cubic metals (Al, Ni, Cu, Ag, Au) is also 
the group of relatively high melting point, but even when allowance is 
made for this the distinction still exists. 

However, Dr. Henderson and I (1951) have recently found that the 
pure f.c. cubic metals, gold and silver, show the hexagonal type of glide 
(‘ easy glide ’) at low temperature when the glide is not too large. Rohm 
and Kochendorfer (1950) have shown that a single crystal of aluminium 
(f.c. cubic), subjected to approximately simple shear, gives the hexagonal 
type of stress-strain curve as against the type normally found with 
cubic metals in tension and, in particular, that found by other workers 
with aluminium under the ordinary test conditions. The difference 
between the mechanical behaviour of the two classes of metals appears, 
then, to be less fundamental than has been suggested. With suitably 
chosen conditions glide in cubic metals can apparently be restricted to 
one set of glide planes, without even local interference by other sets, 
and then the behaviour is of a hexagonal nature. 
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The critical shear stress shows no distinction between the two types 
of metals. Schmid, in the paper quoted, gives = diagram of critical 
shear stress against temperature for four metals (Cd, Zn, Mg, Bi), but 
no systematic comparison with cubic metals appears to gee been made. 
It has been suggested (Andrade 1938, Andrade and Chow 1940) that 
in comparing the mechanical properties of different metals it is advisable — 
to take as the measure of the temperature the ratio 6=T/T,,, where T 
and T’,, are respectively the temperature in question and the temperature 
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of melting, both absolute. Accordingly I have prepared a diagram 
(fig. 1) showing the variation of critical shear stress with 6 for gold and 
silver (Andrade and Henderson), cadmium (Boas and Schmid), magnesium 
(Schmid 1931), zine (Fahrenhorst and Schmid) and bismuth (Georgieff 
and Schmid) and have inserted single values for lead (Makin, unpublished), 
copper (Sachs and Weerts) and mercury (Greenland). The critical shear 
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stress is notoriously a quantity which is very structure-sensitive, but 
the values represented in the diagram are for pure metals and seem 
consistent and reliable. For aluminium the values are so variable that 
I have not included the metal. ; 

It will be seen that the curves for the different metals show a very 
similar course, which they do not do if the critical shear stress o is plotted 
against 7’. The values for the different crystal types show no suggestion 
of grouping according to crystal structure. With magnesium, which is 
hexagonal and normally glides on the basal planes, pyramidal glide on 
(1011) faces is established at about 225° c (@=0-54), which accounts for 
the peculiar form of the curve for this metal. 

According to Orowan (1934) the variation of critical shear stress with 7’ 
is given by 

o=0,—BT}?, 
In fig. 2, o is plotted against 61/*, the points being taken from the smooth 
curves in fig. 1 (which were drawn before this plot was made). For gold, 
silver, cadmium and zine Orowan’s law holds well within the range 
4=0-05 to 0-4, although for higher temperature the experimental values 
are higher than it would indicate. . For bismuth values are not available 
for low 6’s. 
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The Sintering of Powders and Diffusion 


By R. 8S. Barnes 
Atomic Energy Research Establishment, Harwell* 


[Received August 15, 1952] 


Iv is the purpose of this communication to apply the result of some 
recent interdiffusion experiments (Barnes 1952), using solid sandwiches 
of metals which form a complete series of solid solutions, to the density 
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behaviour of mixed powders which also form solid solutions, and to the 
more fundamental problem of the mechanism of the sintering process. 

It has been observed that the density of mixed copper and nickel 
powders decreases initially (by sometimes as much as 20%) before it 
starts to increase on heating at temperatures where interdiffusion can 
occur (Rhines and Meussner 1943, Alexander 1951, Butler and Hoar 
1952). Changes in lattice parameter cannot account for this large 
decrease, nor can the expansion of gas entrapped during compacting, as 
compacts consisting of copper powder only do not behave similarly 
(Butler and Hoar 1952). 

Density changes similar to these have been observed when some 
mixed spinel powders are heated together (Rigby 1952). 

Diffusion experiments have shown that voids form on the copper-rich 
side of the marked join in copper-nickel couples (Seith and Kottmann 
1952, Barnes 1952). Similar voids resulting from the interdiffusion of 
alpha—brass and copper, and aluminium—bronze and copper have been 
found by Biickle and Blin (1952). The formation of these voids in 
material which before diffusion had been quite solid reduced the density 
of the sandwiches by as much as 20% for almost complete diffusion 
(Barnes 1952). These experiments were able to show that the voids 
did not appear as a result of gas, either in the metals themselves or 
entrapped during the welding process. Thus in sintered mized powders. 
an extra volume increase, over and above any caused by included gas, 
must result from the formation of such voids, and it is suggested that this 
is the cause of the initial density decrease, and also the initial decrease in 
hardness and rapid rise in the electrical resistance observed by Rhines 
and Meussner (1943). 

The voids resulting from the interdiffusion process have crystallo- 
graphic shapes, therefore they should be distinguishable from voids. 
formed by the closing of the pores between the original particles, which 
initially would have less regular shapes. A careful mechanical polish or 
microradiographic examination should show this. 

The interdiffusion experiments can best be explained on the basis of a 
vacancy mechanism of diffusion where there is a resultant flow of 
vacancies from the nickel, where they are generated, to the copper, 
where they are lost. The vacancies can be lost in two ways, either by 
annihilation at lattice discontinuities, e.g. grain boundaries and edge 
dislocations, which results in a shift, relative to both copper and nickel 
surfaces, of the markers placed at the original interface ((b) in figure), or by 
coalescing to form voids, which results in a volume increase equal to the 
ae eae aed ve @ spill of the peat relative to the nickel 
to both the copper and_ nickel eS (d) ; peprpeeten ett 
nae, es (d), both these processes must: 
ee ee eee ee parece to form voids is proportional to 
(Barnes 1962) Tae inulin and to +/¢ (where is the annealing time) 

: et number of vacancies which have crossed the 
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marked interface is proportional to 8/ (the shift of the markers relative 
to the nickel reference surface), and this is also proportional to \/t. Thus 
the number of vacancies which have been annihilated at lattice dis- 
continuities must also be proportional to 1/t. The relative shifts suggest 
that approximately equal numbers of vacancies must be lost by each of 
the two processes, and as they are both proportional to \/t it seems 
highly probable that a given proportion of the vacancies meet lattice 
discontinuities and are annihilated, the rest forming voids. Little 
sintering of the voids can occur during interdiffusion as this would be 
expected to give a dependence other than /f. 

Thus, when interdiffusion causes the material to be swpersaturated with 
vacancies, they are lost both by void formation and at lattice dis- 
continuities, the energy needed to increase the area of the internal and 
external surfaces being derived from the process of mixing the two types 
of atom. However, when the material is of homogeneous composition 
{either because the original compact consisted of particles of only one 
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Representation of the cross-section of a couple (a) before the diffusion anneal, 
(6) after diffusion, all the vacancies haveng been annihilated, (c) all the 
vacancies going to form voids, and (d) what happens in practice. 


material or because of complete diffusion) there will be a tendency to 
retain only the equilibrium number of vacancies, and these will not 
enlarge existing voids as there is no source of energy to increase the 
surface area. But vacancies will still be annihilated at lattice dis- 
continuities without increasing the surface energy, locally reducing the 
number of vacancies in the lattice. This loss will be made up by a steady 
generation of vacancies at voids which will diminish in size as a result. 
Small voids would decrease in size more rapidly than large voids, because 
the gain in energy per vacancy generated is larger at small voids. Large 
voids might be expected even to grow at the expense of smaller voids 
because of the resultant decrease in surface energy, but the total volume 
of the voids would steadily decrease, and this was observed by Rhines, 
Birchenall and Hughes (1950). The steady annihilation of vacancies at 
lattice discontinuities will eliminate whole planes of atoms, causing a 
reduction in the total volume. Thus the sintering process probably 
consists of a special form of plastic deformation where the change in 
shape results from a flow of vacancies from the voids to lattice 
discontinuities. 
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In the Fermi—Dirac distribution function, namely 
1 
¢, as is well known, is the thermodynamic potential, or the free energy 
at constant pressure, per electron, and is given by 
C= —T8+ 00, 5 si saa hh ee 
where the letters have their usual significance. Obviously ¢ is a function 
of the temperature, and is determined by the equation 


=2| Ns (E) dE, Oe PA ey 


where » is the number of electrons per unit volume, and N (Z) is the 


density of energy states. For an almost completely degenerate gas, i.e. 
at temperatures 7'’<7'), where 


Ty=blk= sa (5) Supine Seeman 


is the degeneracy-temperature, the integration of (3) yields to a first 
approximation 


SOE aR) clue ee oe Bes = ri 
where 


york? /2C. s aelaetey RW shy fei « NM en toe GS 
* Communicated by the Authors. 
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Since the second term in (5) is merely a correction term, ¢ appearing in 
the denominator of (6) may be replaced by ¢,, or by ¢*, which we shall 
define presently, which are both of nearly the same magnitude as ¢. 

In expression (5) ¢, is frequently referred to as the value of £ at T= 0, 
and for that reason is sometimes regarded as independent of 7’, and 
hence the temperature variation of ¢ is taken to be determined completely 
by the second term, which involves 7’ explicitly. It is the main purpose 
of this note to emphasize that this will be the case only if the density 
of electrons n is kept constant, whereas at constant pressure the density n, 
and hence also @) as will be seen from (4), will vary with the temperature. 
In other words , will not represent the free energy at the same pressure 
at T=—0. The latter energy will be given by 


Co*=Co(1+ 7/3), Aue a a ener ea 


and the expression for the temperature variation of ¢ at constant pressure 
will therefore be given by 
Cee (lays eee eto ea.) «(S) 


Physically ¢)*—€ has the following significance. Let 47 be the thermionic 
work function, as usually defined, of a metal at temperature 7’. If the 
electrons in the metal, i.e. in the condensed phase, can be regarded as an 
almost completely degenerate assemblage, having a finite latent heat of 
evaporation, then the temperature variation of 4 due to the temperature 
variation of € will be given by 


nap n Us eee se ot. (9) 


Incidentally it may be mentioned that thermodynamically 
Tes e 
it—f=T | Bat — | iawn 4. (10) 
wee Cae 


where ¢,, is the specific heat of the electrons in the condensed phase at 
constant pressure. In view of (8) and (10) one obtains 


Ue Ve eee 2!) ss). 9 (LT) 


which is the same expression as for c,. 

It should be mentioned here that if the electronic structure of the 
metal, instead of corresponding to a nearly empty parabolic band, as 
we have taken it to be till now, corresponds to a nearly full parabolic 
band, the expressions for c, and c, will remain the same as before, but 
will now refer to the holes, and hence ¢,*—{ in eqns. (8), (9) and (10) 
should be replaced by ¢—¢)*. In other words ¢ will now increase with 
the increase of the temperature, though the specific heats remain positive 
as they should. 

That in a highly degenerate electron gas, the specific heat at constant 
pressure should be the same as that at constant volume can also be 
demonstrated otherwise. At higher temperatures naturally c, will 
increase more rapidly than c,, and they will tend to the values 3k and 


3k respectively per electron. 
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CXXII. Notices of New Books and Periodicals received 


Descartes and the Modern Mind. By Autpert G. A. Bauz. (New Haven: 
Yale University Press, 1952.) [Pp.492.] Price $ 10.00. 


Tuis is an impressive essay in the history of ideas. The author shows how 
Descartes’ philosophical position arises from the scholastic tradition as embodied 
in and modified by Thomas Aquinas. In developing his thesis Professor Balz 
generously acknowledges his indebtedness to Gilson and Koyré and gives good 
reasons for any disagreement with these scholars. 

His mode of presentation will not be to everybody’s taste. Thus, I believe 
that little help will be derived from the proposal to use “ the name Cartesius to 
signify Reason itself ’’ and to regard René Descartes ‘‘ as one who proposes to 
report what he has learned from Cartesius, the very voice of Reason ”’ (p. 67). 
It would be a pity if this and similar passages were to deter anybody from 
reading a very learned, well documented, and interesting contribution to the 
history of philosophy. See 


ERRATUM 


On the Theory of Beta-Radioactivity, by A. M. Smrru, 1952, Phil. Mag. 
43, 915. 


Figures 3 (a) and 3(b) should be intérchanged. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 
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